Homogeneity test of Global Trends using Chi-Square on Mann-Kendall
Updated 19 Feb 2009
GlobalTrends- Homogeneity tests for multiple seasons and stations. This function will test for trends when seasonality is present and over multiple observation stations, all of which are Chi-square statistics. There are so many statistical tests being done, this function is more like a script or program than a function, but I prefer operating with functions.
This function relies heavily on Matlab's Statistical Toolbox for obtaining Chi-square values and ktaub.m function.
These tests will allow for ties, missing data, and multiple observations per time index, since it uses the enhanced ktaub.m function that was recently updated.
There are numerous narratives as part of the output to the screen providing supporting interpretation of the results.
This function is based on Chapter 17.5 in Gilbert.
[ChiOutput K M sig] = GlobalTrends( datain, alpha )
To stay consistent with the previous trend statistics posted (ktaub, sktt, and SenT), data are expected to be in the following structure:
datain(:,1) = Year
datain(:,2) = season
datain(:,3) = value
datain(:,4) = station
alpha = scaler (e.g. 0.05)
ChiOutput structure is like the following (labels are not included)
Total: Chi-square df p-value
Homogeneity: Chi-square df p-value
Season: Chi-square df p-value
Station: Chi-square df p-value
Station-Season: Chi-square df p-value
Trend: Chi-square df p-value
And depending on significances of Stations, Seasons, and StationSeasons, one of three other outpus may occur:
K: significance of Seasons per station
M: significance of Stations for seasons
And when seasonal trend tests should not be done
sig: individual station-season p-values are given by row
There is a lot of output to the screen as well, but using fprintf, one could easily redirect output to a file.
- Matlab Statistical Toolbox
Richard O. Gilbert, Pacific Northwest National Laboratories,
"Statistical methods for Environmental Pollution Monitoring", 1987, Van Nostrand Reinhold, New York Publishing, ISBN 0-442-23050-8.
One last note. I’ve posted enough functions that it makes sense to constitute forming a Matlab toolbox…so it’s in the works (slowly).
King County- Department of Natural Resources and Parks
Jeff Burkey (2023). Homogeneity test of Global Trends using Chi-Square on Mann-Kendall (https://www.mathworks.com/matlabcentral/fileexchange/22440-homogeneity-test-of-global-trends-using-chi-square-on-mann-kendall), MATLAB Central File Exchange. Retrieved .
MATLAB Release Compatibility
Platform CompatibilityWindows macOS Linux
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!
Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
I've gotten a couple of inquiries asking for an example dataset to use. I've added an example MAT file along with the m-file.
Updated Title to be more descriptive.
Updated description text. Data are not required to be linear per se, but formated like time-series. Up comming is a Log-person Type III function (aka USGS 17B) flood fequency analysis. Also, I've posted enough functions that a toolbox is coming.