# NOT RUN {
# Calculate a covariance matrix on a null community
data(knz_001d)
a1_cyclic <- cyclic_shift(subset(knz_001d, subplot == "A_1"),
time.var = "year",
species.var = "species",
abundance.var = "abundance",
FUN = cov,
bootnumber = 10)
# Return CI on a1_cyclic
confint(a1_cyclic)
# }
Run the code above in your browser using DataLab