Description Usage Arguments Details Value Author(s) References See Also
This function performs a standard normal homogeneity test, but allows for unequally spaced observations in time.
1 2 |
data |
The data to be analyzed for changepoints. |
period |
The SNHT works by calculating the mean of the data on the previous period observations and the following period observations. Thus, this argument controls the window size for the test statistics. |
time |
Numeric vector specifying times for the observations. If not supplied, it is assumed that each observation occurs on one time period. If supplied, then the algorithm will create a new dataset with the same number of observations for each time unit by adding missing values. |
estimator |
See ?robustSNHT |
scaled |
See ?snht. |
rmSeasonalPeriod |
See ?snht. |
The SNHT works by calculating the mean of the data on the previous period and on the following period. The test statistic at each observation is then computed as described in Haimberger (2007). Essentially, though, it just compares the means of these two periods and normalizes by the standard deviation.
Note: if there are not enough observations both before and after the current observation, no test is performed.
Large values of the test statistic suggests the presence of a changepoint. Haimberger (see references) suggests values larger than 100 should be considered changepoints. However, this does not apply if scaled = TRUE.
Returns a data.frame, with columns score, leftMean, and rightMean, and time. Statistic is the SNHT test statistic described above, and leftMean (rightMean) are the means to the left (right) of the current observation.
Note that new (missing) observations were introduced to the dataset to ensure the same number of observations occur per day.
Josh Browning (jbrownin@mines.edu)
L. Haimberger. Homogenization of radiosonde temperature time series using innovation statistics. Journal of Climate, 20(7): 1377-1403, 2007.
Other snht functions: robustSNHT
,
snht
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.