pool.var: Compute a pooled variance

pool.varR Documentation

Compute a pooled variance


Computes a pooled variance from samples in a paleontological time-series


pool.var(y, nn = NULL, minN = NULL, ret.paleoTS = FALSE)



either a paleoTS object, or a vector of sample variances


a vector of sample sizes


sample size below which variances are replaced with pooled variances. See Details.


if TRUE, a paleoTS object is returned. If FALSE, the value of the pooled variance is returned.


A pooled variance of a set of populations is the weighted average of the individual variances of the populations, with the weight for each population equal to its sample size minus one.

For many kinds of traits, variation levels tend to be similar among closely related populations. When this is true and sample sizes are low, much of the observed differences in variance among samples will be due to the high noise of estimated the variances. Replacing the observed variances of all populations (or only those with nn < minN) with the estimated pooled variance can reduce this noise.


if ret.paleoTS = TRUE a paleoTS object with all (or some) variances replaced with the pooled variance; otherwise the pooled variance


cant_all <- pool.var(cantius_L, ret.paleoTS = TRUE)   # replace all variances with pooled variance
cant_n5  <- pool.var(cantius_L, minN = 5, ret.paleoTS = TRUE)  # replace only pops with n < 5

paleoTS documentation built on Aug. 9, 2022, 1:06 a.m.