raplot.paired {MLpreemption}

Rank-abundance plot for two samples

Description

Function raplot.paired creates a rank-abundance plot on screen, and can show decaying lines with uncertainty zones for two samples fitted by maximum likelihood.

Usage

raplot.paired(x, y, xlab = "Species rank", ylab = "log (Relative abundance)",
              main = "Rank-abundance", sym = c(1, 2), alpha = 0.05)

Arguments

x

Count vector of the first sample

y

Count vector of the second sample

xlab

Label x axis

ylab

Label y axis

main

Main title for the plot

sym

Symbols for first and second sample (c(1,2) by default)

alpha

Significance level (0.05 by default)

Value

NULL

Author(s)

Jan Graffelman ([e-mail])

References

Graffelman, J. (2021) Maximum likelihood estimation of the geometric niche preemption model

See Also

raplot

Examples

data("Mehrabi")
raplot.paired(Mehrabi[,1],Mehrabi[,2])