bet.het.aSE {DendroSync}

Between-group synchrony for a heteroscedastic unstructured model

Description

The function calculates the between-group synchrony (a^) and standard error (SE) for a heteroscedastic unstructured model (mHeUN).

Usage

bet.het.aSE(model)

Arguments

model

a class "lme" unstructured model (mHeUN) produced by dendro.varcov with homoscedastic equals FALSE.

Details

The function calculates between-group synchrony for a heteroscedastic unstructured model (mHeUN).

Value

The function returns a matrix containing between-group synchrony and SE for each combination of varGroup levels. This function is used internally in sync.

Author(s)

Josu G. Alday, Tatiana A. Shestakova, Victor Resco de Dios, Jordi Voltas

References

Shestakova, T.A., Aguilera, M., Ferrio, J.P., Gutierrez, E. & Voltas, J. (2014). Unravelling spatiotemporal tree-ring signals in Mediterranean oaks: a variance-covariance modelling approach of carbon and oxygen isotope ratios. Tree Physiology 34: 819-838.

Shestakova, T.A., Gutierrez, E., Kirdyanov, A.V., Camarero, J.J., Genova, M., Knorre, A.A., Linares, J.C., Resco de Dios, V., Sanchez-Salguero, R. & Voltas, J. (2016). Forests synchronize their growth in contrasting Eurasian regions in response to climate warming. Proceedings of the National Academy of Sciences of the United States of America 113: 662-667.

See Also

sync for a clear description of synchrony evaluation.

Examples

## 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