## Calculate between-group synchrony and SE
##for a heteroscedastic unstructured model for conifersIP data:
data(conifersIP)
#Fit the heteroscedastic set of varcov models (mBE, mHeNE, mHeCS, mHeUN)
# using taxonomic grouping criteria (ie. Species)
ModHt <- dendro.varcov(TRW ~ Code, varTime = "Year", varGroup = "Species",
data = conifersIP, homoscedastic = FALSE)
#between-group synchrony and SE for each varGroup stratum combination
# in heteroscedastic unstructured models.
bet.het.aSE(ModHt$mHeUN)#Unstructured model
Run the code above in your browser using DataLab