# Read the data for Northern and Southern European cerambycids
data(ceram.s)
data(ceram.n)
# Resample 100 times the multiple-site dissimilarities
# for 10 countries.
beta.ceram.s<-beta.sample(ceram.s, index.family="sor", sites=10, samples=100)
beta.ceram.n<-beta.sample(ceram.n, index.family="sor", sites=10, samples=100)
# Plot the distributions of beta.SIM in Southern Europe (red)
# and Northern Europe (blue)
hist(beta.ceram.s$sampled.values$beta.SIM, breaks=10, col="red", xlim=c(0,1))
hist(beta.ceram.n$sampled.values$beta.SIM, breaks=10, col="blue", add=TRUE)
Run the code above in your browser using DataLab