varEvol: Plot mean variance versus the number of time points.

In hallucigenia-sparsa/seqtime: Time Series Analysis of Sequencing Data

Plot mean variance versus the number of time points.

Description

Plot the mean variance across all taxa for an increasing number of time points. If groups are given, compute the mean of beta diversity (assessed with Bray Curtis) across samples for each time point. Groups are interpreted as replicates, such that each group is assumed to cover the same time points in the same order.

Usage

varEvol(x, groups = c(), plot = TRUE, sd = FALSE)

Arguments

x	a community matrix with rows representing taxa and columns time points
groups	a vector with group assignments and as many entries as there are samples
plot	do the plot
sd	use the standard deviation instead of the variance

Value

a list with variances (or standard deviations if sd is true) and the intersection, slope, pval and adjR2 of the linear regression of variances against time

Examples

N=20
ts=glv(N=N,generateA(N))
out=varEvol(ts)

hallucigenia-sparsa/seqtime documentation built on Jan. 9, 2023, 11:53 p.m.