| modelFrame | R Documentation |
This function develops recursive evaluation of functions for one-level modeling (FOLM) and LME detrending of dendroclimatic chronologies.
modelFrame(rd, fn = list("rtimes",
"scacum", "amod"),
lv = list(2, 1, 1),
form = "tdForm",
...)
rd |
|
fn |
|
lv |
|
form |
|
... |
Further arguments in |
Defaults model fluctuations in
tree-ring width chronologies via recursive
implementation of four FOLM:
rtimes, scacum,
amod, and
frametoLme. Nevertheless,
other FOLM can be implemented to model
aridity-index fluctuations(see example with
climatic data). Processed chronologies are
detrended with lme function
and other nlme methods
. Internal algorithm uses
shiftFrame
arguSelect and
ringApply
functions. Consequently, arguments that are
not iterated over factor-level labels in the
processed data are specified in 'MoreArgs'
lists (see examples). Arguments in
modelFrame objects can be updated
with update function.
Threefold list with fluctuations in fluc,
groupedData object in model, and model call in
call.
Wilson Lara <wilarhen@gmail.com> [aut, cre] (<https://orcid.org/0000-0003-3527-1380>), Felipe Bravo <fbravo@pvs.uva.es> [aut] (<https://orcid.org/0000-0001-7348-6695>)
Lara W., F. Bravo, D. Maguire. 2013. Modeling patterns between drought and tree biomass growth from dendrochronological data: A multilevel approach. Agric. For. Meteorol., 178-179:140-151.
##TRW chronology (mm) and inside-bark radii
data(Pchron,envir = environment())
## Parameters of allometric model to compute Diameter at Breast
## Height over bark (DBH, cm) from diameter inside bark (dib, cm)
## and Total Tree Biomass (TTB, kg tree -1 ) from DBH (Lara
## et. al. 2013):
biom_param <- c(2.87, 0.85, 0.05, 2.5)
## Modeling tree-biomass fluctuations while accounting for
## within-plot source variability (see defaults in "modelFrame"
## function)
## \donttest{
## trwf <- modelFrame(Pchron,
## to = 'cm',
## MoreArgs = list(mp = c(2,1, biom_param)),
## log.t = FALSE,
## on.time = FALSE)
## }
## Climatic records:
data(Temp,envir = environment())
data(Prec,envir = environment())
## Aridity-index fluctuations:
## \donttest{
## aif <- modelFrame(rd = list(Prec, Temp),
## fn = list('moveYr','wlai'),
## lv = list('year','year'),
## form = 'lmeForm')
## summary(aif$'model')
## }
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.