Description Usage Arguments Details Value Author(s) References See Also Examples
This function calculates moving response and correlation functions from treering chronologies and monthly climatic data. The calculation is performed repeatedly for consecutive time windows. Function parameters may be bootstrapped to calculate their significance and confidence intervals.
1 2 3 
chrono 

clim 

method 
string specifying the calculation method. Possible values are “response” and “correlation”. Partial strings are ok. 
start 
integer value to determine the first month to be used as a predictor in the response or correlation function. A negative value denotes a start month from previous year, a positive value denotes a start month from current year. 
end 
integer value to determine the last month to be used as a predictor in the response or correlation function. A negative value denotes an end month from previous year, a positive value denotes an end month from current year. 
timespan 
integer vector of length 2 specifying the time interval (in years) to be considered for analysis. Defaults to the maximum possible interval. 
vnames 
character vector with variable names. defaults to corresponding column names of 
sb 
logical flag indicating whether textual status bar should be suppressed. Suppression is recommended for e.g. Sweave files. 
win.size 
integer giving the window size for each recalculation. 
win.offset 
integer giving the number of years between each window start. 
startlast 
logical flag indicating whether the first window should start at the rear end (youngest part of the series) or not. 
boot 
logical flag indicating whether bootstrap resampling is to be performed. 
ci 
numerical value to set the test level for significance test (values 0.01, 0.05 and 0.1 are allowed); the confidence intervals are adapted accordingly. 
The functions dcc
and mdcc
clone the
functionality of programme DENDROCLIM2002 (Biondi and Waikul, 2004), and
will calculate bootstrapped (and nonbootstrapped) moving (mdcc
and static (dcc
) response and correlation functions in a similar
fashion as described in the above mentioned paper.
In case of response function analysis 1000 bootstrap samples are taken from the original distribution and an eigen decomposition of the standardized predictor matrix is performed. Nonrelevant eigenvectors are removed using the PVP criterion (Guiot, 1990), principal component scores are then calculated from the matrices of reduced eigenvectors and standardized climatic predictors. Response coefficients are found via singular value decomposition, and tested for significance using the 95% percentile range method (Dixon, 2001). In case of correlation function analysis, the coefficients are Pearson's correlation coefficients. The same method for significance testing is applied.
Input chronology data can be a data.frame
such as produced by
function chron
of package dplR. It has to be a data.frame
with at least one column containing the treering indices, and the
corresponding years as rownames
.
For climatic input data, there are three possibilities: Firstly, input
climatic data can be a data.frame
or matrix
consisting of
at least 3 rows for years, months and at least one climate parameter in
the given order. Secondly, input climatic data can be a single
data.frame
or matrix
in the style of the original
DENDROCLIM2002 input data, i.e. one parameter with 12 months in one row,
where the first column represents the year. Or thirdly, input climatic
data can be a list of several of the latter described data.frame
or matrices
. As an internal format dispatcher checks the format
automatically, it is absolutely necessary that in all three cases, only
complete years (months 112) are provided. It is not possible to mix
different formats in one go.
The window for response/correlation function analysis is specified via
start
and end
, where e.g. 4 means previous April etc.
The window size for moving response and correlation functions is set via
win.size
, and the distance from one window start to the next is
set with the parameter win.offset
. Parameter startlast
indicates, wether the first window is started from the rear (youngest
part) of the series or not.
Bootstrapping (boot
) is by default disabled to get the results
faster.
A list
containing data.frame
s for coefficients and
confidence intervals for each parameter and time window used for the
moving functions.
Christian Zang
Biondi, F. & Waikul, K. (2004) DENDROCLIM2002: A C++ program for statistical calibration of climate signals in treering chronologies. Computers & Geosciences 30:303311
Dixon, P.M. (2001) Bootstrap resampling. In: ElShaarawi, A.H., Piegorsch, W.W. (Eds.), The Encyclopedia of Environmetrics. Wiley, New York.
Guiot, J. (1991) The boostrapped response function. TreeRing Bulletin 51:3941
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19  data(muc.clim) # climatic data
data(muc.spruce) # spruce data
# calculate and plot moving response function
dc.mov1 < mdcc(muc.spruce, muc.clim)
mdcplot(dc.mov1)
# calculate and plot moving correlation function with different window
# parameters
data(rt.spruce)
data(rt.temp)
data(rt.prec)
dc.mov2 < mdcc(rt.spruce, list(rt.temp, rt.prec),
vnames = c("temp", "prec"),
method = "corr", win.size = 20,
win.offset = 5)
mdcplot(dc.mov2)


  0%

==  3%

====  6%

=======  9%

=========  12%

===========  16%

=============  19%

===============  22%

==================  25%

====================  28%

======================  31%

========================  34%

==========================  38%

============================  41%

===============================  44%

=================================  47%

===================================  50%

=====================================  53%

=======================================  56%

==========================================  59%

============================================  62%

==============================================  66%

================================================  69%

==================================================  72%

====================================================  75%

=======================================================  78%

=========================================================  81%

===========================================================  84%

=============================================================  88%

===============================================================  91%

==================================================================  94%

====================================================================  97%

====================================================================== 100%

  0%

============  17%

=======================  33%

===================================  50%

===============================================  67%

==========================================================  83%

====================================================================== 100%
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.