有关MK检验的问题

baonan0811

请问有人对魏凤英老师在《现代气候统计诊断与预测技术》一书中讲解MK检验时用到的上海1900-1990年91年平均温度的例子做过实习吗？

为什么我用书中的数据作出的检验结果和书中不一致，但是用同事已经做过的MK检验数据和结果对比是一致的。
出现这种情况要么就是我和同事关于MK检验的算法编写有错误，要么就是书中的数据和检验结果不一致。

baonan0811

第二章图中红色的为UB线黑色的为UF线

对比两张图中UB和UF交点位置不一样，书中为1921，而我检验的结果是1931

本以为可能是我的程序写的有问题，但是拿同事已经做过检验的数据经行对比，得出的UF\UB数据是一致的。

baonan0811

这里是有关MK检验的代码

1. class Mann_Kendall
   
2.     {        
   
3.         private float[] _myUF;
   
4.         private float[] _myUB;
   
5. 
   
6.         public float[] UF
   
7.         { get { return _myUF; } }
   
8. 
   
9.         public float[] UB
   
10.         { get { return _myUB; } }
    
11. 
    
12.         public Mann_Kendall(List<float> list)
    
13.         {
    
14.             float[] myArry = list.ToArray();
    
15.             list.Reverse();
    
16.             float[] myRevArry = list.ToArray();
    
17. 
    
18.             CalculateUF(myArry);
    
19.             CalculateUB(myRevArry);
    
20.         }
    
21. 
    
22.         public Mann_Kendall(float[] arry)
    
23.         {           
    
24.             float[] tempArry = new float[arry.Length];
    
25.             Array.Copy(arry, tempArry, arry.Length);
    
26.             Array.Reverse(tempArry);
    
27. 
    
28.             CalculateUF(arry);
    
29.             CalculateUB(tempArry);
    
30.         }
    
31. 
    
32.         private void CalculateUF(float[] x)
    
33.         {
    
34.             _myUF = GetResult(x);
    
35.         }
    
36. 
    
37.         private void CalculateUB(float[] x)
    
38.         {
    
39.             float[] result = GetResult(x);
    
40.             List<float> listUB = new List<float>();
    
41.             foreach (float item in result.Reverse())
    
42.             {
    
43.                 listUB.Add(-item);
    
44.             }
    
45.             _myUB = listUB.ToArray();
    
46.         }
    
47. 
    
48.         private float[] GetResult(float[] x)
    
49.         {
    
50.             List<int> listR = new List<int>();
    
51.             for (int i = 0; i < x.Length; i++)
    
52.             {
    
53.                 int sum = 0;
    
54.                 for (int j = 0; j <= i; j++)
    
55.                 {
    
56.                     sum = x[i] > x[j] ? sum + 1 : sum;
    
57.                 }
    
58.                 listR.Add(sum);
    
59.             }
    
60. 
    
61.             float[] U = new float[listR.Count];
    
62.             for (int i = 0; i < listR.Count; i++)
    
63.             {
    
64.                 if (i == 0)
    
65.                 { U[i] = 0; }
    
66.                 else
    
67.                 {
    
68.                     int s = 0, n = i + 1;
    
69.                     float e = 0, var = 0;
    
70.                     e = n * (n - 1) / 4f;
    
71.                     var = n * (n - 1) * (2 * n + 5) / 72f;
    
72.                     for (int j = 0; j <= i; j++)
    
73.                     {
    
74.                         s = s + listR[j];
    
75.                     }
    
76.                     U[i] = Convert.ToSingle(Math.Round((s - e) / Convert.ToSingle(Math.Sqrt(Convert.ToDouble(var))), 2));
    
77.                 }
    
78.             }
    
79. 
    
80.             return U;
    
81.         }  
    
82. 
    
83.     }

复制代码

baonan0811

同样的数据，做的滑动T检验，取滑动步长为10年，虚线为0.05水准线-2.89。
分析结果显示突变年份应该在1930年左右，
可以看出我的MK检验结果应该没问题

言深深

麻烦楼主把你的测试数据传上来吧，我做过mk，曾经也把mk的程序发在论坛了，刚刚，不，从下午开始就想来测试一下，但是却发现敲了几个小时，还是没敲对！！！！！！！！！！！

baonan0811

言深深 发表于 2011-12-8 22:03

                               
登录/注册后可看大图

麻烦楼主把你的测试数据传上来吧，我做过mk，曾经也把mk的程序发在论坛了，刚刚，不，从下午开始就想来测试 ...

数据在单位电脑，明天上班发。
麻烦版主帮忙检验下看看到底是我的问题还是书上的图有问题。

baonan0811

数据如下，麻烦版主帮忙测试一下

majianyong

baonan0811 发表于 2011-12-9 08:42

                               
登录/注册后可看大图

数据如下，麻烦版主帮忙测试一下

楼主，我给你做了一下MK检验，结果与魏凤英老师的相同，结果见附件j:\m-k.png

majianyong

累计距平结果见下图：大概1935年为突变明显年，所以各种统计方法结果不一致，多互相比较，采取一种有效的，另外亦可以用小波分析查看突变年，这里我就不做了j:\累计距平.png

majianyong

另外，按照你的滑动T检验，你取的是两边10年步长，如1900-1909--1910-1919，求得的第一个统计值赋给了1909年，既然是滑动，所以这里我赋给1901-1909里面任何一年都是可以的，所以你说通过滑动T检验的结果来验证你的MK结果是正确的是不合理的。还是那句话，多比较几种检验方法，找个合适的。