From e05283245623be01d9cc9c96e6c551c0308f9826 Mon Sep 17 00:00:00 2001 From: AndyBunn Date: Mon, 24 Feb 2025 13:53:36 -0800 Subject: [PATCH] New deve version 1.7.9 and some small changes to help files for sss and rwi.stats --- ChangeLog | 8 +++++++- DESCRIPTION | 2 +- man/rwi.stats.running.Rd | 6 +++++- man/sss.Rd | 37 +++++++++++++++++++++++++++++++------ 4 files changed, 44 insertions(+), 9 deletions(-) diff --git a/ChangeLog b/ChangeLog index 5225938..b0c8919 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,4 +1,10 @@ -* CHANGES IN dplR VERSION 1.7.8 +* CHANGES IN dplR VERSION 1.7.9 + + +File: rwi.stats ans sss help +---------------- + +* Changed verbiage after emails with Stefan Klesse and Allan Buras File: powt diff --git a/DESCRIPTION b/DESCRIPTION index 35f5e93..3e46335 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -2,7 +2,7 @@ Encoding: UTF-8 Package: dplR Type: Package Title: Dendrochronology Program Library in R -Version: 1.7.8 +Version: 1.7.9 Authors@R: c(person("Andy", "Bunn", role = c("aut", "cph","cre", "trl"), email = "bunna@wwu.edu"), person("Mikko","Korpela", role = c("aut", "cph", "trl")), person("Franco", "Biondi", role = c("aut", "cph")), diff --git a/man/rwi.stats.running.Rd b/man/rwi.stats.running.Rd index 57f9890..75394b7 100644 --- a/man/rwi.stats.running.Rd +++ b/man/rwi.stats.running.Rd @@ -133,7 +133,11 @@ rwi.stats.legacy(rwi, ids=NULL, period=c("max", "common")) output. The signal-to-noise ratio is calculated following Cook and Pederson (2011). - Note that Buras (2017) cautions against using the expressed population signal as a statistic to determine the whether a chronology correctly represents the population signal of a data set. He reccomends the use of subsample signal strength (\code{\link{sss}}) over EPS. + Please note that Buras (2017) cautions against the use of expressed + population signal (EPS) for evaluating the climate-reconstruction potential + of a tree-ring sample. Instead, he recommends the use of subsample signal + strength (\code{\link{sss}}) for evaluating the loss of predictive power + back in time when sample replication drops. If desired, the \code{\var{rwi}} can be filtered in the same manner as the family of cross-dating functions using \code{\var{prewhiten}} and diff --git a/man/sss.Rd b/man/sss.Rd index 64c5b20..7fb9450 100644 --- a/man/sss.Rd +++ b/man/sss.Rd @@ -22,21 +22,46 @@ sss(rwi, ids = NULL) \code{data.frame(tree=1:ncol(\var{rwi}), core=rep(1, ncol(\var{rwi})))}. } } \details{ - This calculates subsample signal strength (sss) following equation 3.50 in Cook and Kairiukstis (1990) but using notation from Buras (2017) because writing the prime unicode symbol seems too difficult. The function calls \code{\link{rwi.stats}} and passes it the arguments \code{ids} and \code{prewhiten}. + This calculates subsample signal strength (sss) following equation 3.50 in + Cook and Kairiukstis (1990) but using notation from Buras (2017) because + writing the prime unicode symbol seems too difficult. The function + calls \code{\link{rwi.stats}} and passes it the arguments \code{ids} + and \code{prewhiten}. - To make better use of variation in growth within and between series, an appropriate mask (parameter \code{\var{ids}}) should be provided that identifies each series with a tree as it is common for dendrochronologists to take more than one core per tree. The function \code{\link{read.ids}} is helpful for creating a mask based on the series \acronym{ID}. + To make better use of variation in growth within and between series, an + appropriate mask (parameter \code{\var{ids}}) should be provided that + identifies each series with a tree as it is common for dendrochronologists + to take more than one core per tree. The function \code{\link{read.ids}} is + helpful for creating a mask based on the series \acronym{ID}. - Subsample signal strength is calculated as \eqn{\frac{n[1+(N-1)\bar{r}]}{N[1+(n-1)\bar{r}]}}{n*(1+(N-1)*rbar) / N*(1+(n-1)*rbar)} where \code{n} and \code{N} are the number of cores or trees in the subsample and sample respectively and \code{rbar} is mean interseries correlation. If there is only one core per tree \code{n} is the sample depth in a given year (\code{rowSums(!is.na(rwi))}), \code{N} is the number of cores (\code{n.cores} as given by \code{\link{rwi.stats}}), and \code{rbar} is the mean interseries correlation between all series (\code{r.bt} as given by \code{\link{rwi.stats}}). If there are multiple cores per tree \code{n} is the number of trees present in a given year, \code{N} is the number of trees (\code{n.trees} as given by \code{\link{rwi.stats}}), and \code{rbar} is the effective mean interseries correlation (\code{r.eff} as given by \code{\link{rwi.stats}}). + Subsample signal strength is calculated as \eqn{\frac{n[1+(N-1)\bar{r}]}{N[1+(n-1)\bar{r}]}}{n*(1+(N-1)*rbar) / N*(1+(n-1)*rbar)} + where \code{n} and \code{N} are the number of cores or trees in the + subsample and sample respectively and \code{rbar} is mean interseries + correlation. If there is only one core per tree \code{n} is the sample + depth in a given year (\code{rowSums(!is.na(rwi))}), \code{N} is the + number of cores (\code{n.cores} as given by \code{\link{rwi.stats}}), + and \code{rbar} is the mean interseries correlation between all series + (\code{r.bt} as given by \code{\link{rwi.stats}}). If there are multiple + cores per tree \code{n} is the number of trees present in a given year, + \code{N} is the number of trees (\code{n.trees} as given by + \code{\link{rwi.stats}}), and \code{rbar} is the effective mean interseries + correlation (\code{r.eff} as given by \code{\link{rwi.stats}}). -Readers interested in the differences between subsample signal strength and the more commonly used (running) expressed population signal should look at Buras (2017) on the common misuse of the expressed population signal as well as Cook and Pederson (2011) for a more general approach to categorizing variability in tree-ring data. + Readers interested in the differences between subsample signal strength and + the more commonly used (running) expressed population signal should look at + Buras (2017) on the common misuse of the expressed population signal as well + as Cook and Pederson (2011) for a more general approach to categorizing + variability in tree-ring data. } -\value{ A \code{numeric} containing the subsample signal strength that is the same as number if rows of\code{rwi}. +\value{ A \code{numeric} containing the subsample signal strength that is +the same as number if rows of\code{rwi}. } \references{ - Buras, A. (2017) A comment on the Expressed Population Signal. Dendrochronologia 44:130-132. + Buras, A. (2017) A comment on the Expressed Population Signal. + Dendrochronologia 44:130-132. Cook, E. R. and Kairiukstis, L. A., editors (1990) \emph{Methods of Dendrochronology: Applications in the Environmental Sciences}.