Hi Eddie, I don't have the vegan package, but this may help: accum<-function(x,y) return(x+(y-x)/7) nspec<-90 for(i in 2:50) nspec[i]<-accum(nspec[i-1],220) plot(nspec/nspec[1]*100,type="l",xlab="Sites", ylab="Species accumulation (%)")
Jim On Wed, Jun 17, 2020 at 1:15 AM Eddie Tsyrlin <tsyrl...@student.unimelb.edu.au> wrote: > > I need to express species accumulation curve in percentage terms, i.e. the > vertical axis is from 1% to 100% > I have community data similar to BCI (from vegan package). > I can construct the 'usual' species curve (see below) but I need to convert > the species number to percentage. > > data("BCI") > bcispec<-specaccum(BCI, method = "random", > permutations = 100) > plot(bcispec) > > -- > Eddie Tsyrlin, PhD Candidate | PEARG | Building 404, Grd Floor, Rm G21 > > School of Biosciences | University of Melbourne > *Post:* Bio21, Flemington Rd, Parkville VIC 3052 > *m: *0413 995525 > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.