#Styracaceae number of spp per genus vectorized


Styracaceae<-c(3,2,1,2,3,1,1,2,2,5,8, 87) #numbers of species in same order as genus list
Styracaceae<-sort(Styracaceae)
num_genera<-length(Styracaceae)
nsp<-sum(Styracaceae)
list2<-c(1:(nsp+1))

mem<-NULL

for(i in 1:100){
breaks<-sample(list2,(num_genera-1)) #creating break points that are integers
breaks<-sort(breaks)
nl<-c(breaks,(nsp+1))
n2<-c(1,breaks)
r<-sort(nl-n2)
mem<-rbind(mem,r)
}

plot(1:num_genera,Styracaceae, xlab="Rank genus based on no. spp.",ylab="Number of species",ylim=c(0,nsp),col="blue", pch=7)
for (i in 1:100) points(1:num_genera,mem[i,],col="red",pch=4)


random_is_bigger<-NULL
for (i in 1:num_genera){
random_is_bigger<-c(random_is_bigger, length(which(mem[,i]> Styracaceae[i]))/100)
}
random_is_bigger


list3<-c(1:(nsp-1))
mem<-NULL
for(i in 1:100){
breaks<-sample(list3,(num_genera-1)) 
breaks<-sort(breaks)
d<-c(min(list3)-1,breaks)
e<-c(breaks,max(list3)+1)
r<-sort(e-d)
mem<-rbind(mem,r)
}
#mem
#rowSums(mem)