###################################################
### code chunk number 3: Cs_mci_002
###################################################
covariates <- rbind(
Temp = fulldat[years, "Temp"],
TP = fulldat[years, "TP"]
)
# demean the covariates
the.mean <- apply(covariates, 1, mean, na.rm = TRUE)
covariates <- covariates - the.mean
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.