Re: [R] species accumulation curve with percentages

Hi Eddie, Upon reading your initial request more carefully, the last command should be: plot(nspec/nspec[50]*100,type="l",xlab="Sites", ylab="Species accumulation (%)",ylim=c(1,100)) Jim On Wed, Jun 17, 2020 at 8:23 AM Jim Lemon wrote: > > Hi Eddie, > I don't have the vegan package, but this

Re: [R] species accumulation curve with percentages

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 wrote:

Re: [R] species accumulation curve with percentages

Probably: ?sum ?cumsum e.g. > x <- runif(100, 10,20) > cumsum(x)/sum(x) *100 Cheers, Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Jun 16, 2020 at

[R] species accumulation curve with percentages

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")