Learn R Programming

DendroSync (version 0.1.4)

bet.aSE: Between-group synchrony for a homoscedastic unstructured model

Description

The function calculates the between-group synchrony (a^) and standard error (SE) for a homoscedastic unstructured model or full model (mUN).

Usage

bet.aSE(model)

Arguments

model

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

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.

Details

The function calculates between-group synchrony for a homoscedastic unstructured model (mUN).

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

Run this code
# NOT RUN {
## Calculate between-group synchrony and SE
 # for an unstructured homocedastic model of conifersIP data:
 data(conifersIP)
 
 #Fit the homoscedastic set of varcov models (mBE, mNE, mCS, mUN)
 # using taxonomic grouping criteria (ie. Species)
 ModHm <- dendro.varcov(TRW ~ Code, varTime = "Year", varGroup = "Species", 
                        data = conifersIP, homoscedastic = TRUE)
 summary(ModHm)
   
 #Obtain the unstructured model between-group synchrony and SE
 # for each varGroup stratum.
 bet.aSE(ModHm$mUN)#Unstructure
   
# }

Run the code above in your browser using DataLab