muleMan {BIOdry}

Multilevel dendroclimatic correlograms

Description

Multivariate correlograms between TRW fluctuations and climatic fluctuations.

Usage

muleMan(rd, cd, rd.var = NULL, cd.var = NULL, ...)

Arguments

rd

dataframe or groupedData. TRW fluctuations such as that produced by modelFrame.

cd

dataframe or groupedData. Aridity-index fluctuations such as that produced by modelFrame.

rd.var

character or NULL. Column name of the TRW fluctuations to be compared. If NULL then the first column is processed.

cd.var

character or NULL. Column name of the aridity-index fluctuations to be compared. If NULL then the first column is used.

...

Further arguments in mgram

Details

Function mgram in package ecodist is implemented to compare the dendroclimatic fluctuations. Models being compared should have common higher-level factors (see example).

Value

data.frame object of multivariate correlations.

Author(s)

Wilson Lara <[e-mail]>, Felipe Bravo <[e-mail]>

References

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.

Examples

##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)
 
 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:
 
     aif <- modelFrame(rd = list(Prec, Temp),
                       fn = list('moveYr','wlai'),
                       lv = list('year','year'),
                       form = 'lmeForm')
 
##Multivariate comparison:
 
 mcomp <- muleMan(trwf,
                     aif,
                  nperm = 10^3)
 str(mcomp)