# to calculate standard error of OTU and sequence counts for each host library(plotrix) count=read.csv("Desktop/Niwot/F_manuscript/edits/analyses/alpha/stats.csv",header=TRUE) map=read.csv("Desktop/Niwot/F_manuscript/edits/analyses/NiwotDGF_map3.csv", header=TRUE) count$host<-map$Host se1<-aggregate.data.frame(count$Count,by=list(count$host),std.error) se2<-aggregate.data.frame(count$otus,by=list(count$host),std.error)