data(cliff.sp)
siteCD1<-data.frame(t(cliff.sp[1,]))
#Shannon-Weiner diversity
SW<-function(data){
d<-data[data!=0]
p<-d/sum(d)
-1*sum(p*log(p))
}
bootstrap(siteCD1[,1],SW,R=1000,matrix=FALSE)Run the code above in your browser using DataLab