This function performs a standard normal homogeneity test, but allows for unequally spaced observations in time.
robustSNHTunequal(data, period, time, estimator = NULL, scaled=TRUE
,rmSeasonalPeriod = Inf)The data to be analyzed for changepoints.
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.
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.
See ?robustSNHT
See ?snht.
See ?snht.
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.
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.
L. Haimberger. Homogenization of radiosonde temperature time series using innovation statistics. Journal of Climate, 20(7): 1377-1403, 2007.
Other snht functions: robustSNHT,
snht