<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />大家好。 我正在嘗試運行 Mann-Kendall 來測試時間趨勢是否顯著。 在 JMP 手冊或在線任何內容中都找不到它。 有人有想法麼？
謝謝！
謝謝女士和戴夫！ 我不確定 Kendall 的 Tau 和 Mann-Kendall 測試是否相同（和/或相關）。 感謝您提供有關多變量測試的信息，然後也感謝您提供腳本！
即使 Mann-Kendall 檢驗統計量和 Kendall Tau b 的定義經常用不同的術語來描述（分別從時間序列和相關性的角度），據我所知，它們本質上是相同的。
下面是我很久以前寫的一個腳本。 它計算並顯示 Theil 斜率並執行 Mann-Kendall 檢驗。 圖中顯示了一條參數線性回歸線，僅用於比較。 我計劃通過包含“By”函數以及解釋季節性模式的可能性來概括腳本。 但我還沒有時間實現這些計劃。 我希望這會給你一些想法。
dt =當前數據表() ;
dlg =列對話框(
“從數據表中估計 Theil 斜率和 Kendall Tau” ||
字符( dt <<得到名字) ,
是=列列表( “選擇 Y” ,
最大科爾( 1 ) ,
明科爾( 1 ) ,
數據類型(數字)
) ,
Xs =列列表( “選擇 X” ,
最大科爾( 1 ) ,
明科爾( 1 ) ,
數據類型(數字)
) ,
團體=列列表( “通過（不工作）” ,最大科爾( 1 ) )
) ;
結腸=柱子( dt , dlg [ “是” ] ) ;
colX =柱子( dt , dlg [ “X” ] ) ;
科爾組=柱子( dt , dlg [ “團體” ] ) ;
//計算泰爾斜率
n = N行( dt ) ;
連續下坡= Ĵ ( n , n , . ) ;
為了(一世= 1 ,一世< n ,一世++,
為了( j =一世+ 1 , j <= n , j ++,
連續下坡[一世, j ] = (結腸[ j ] -結腸[一世]) / ( colX [ j ] - colX [一世])
)
) ;
泰爾=分位數( 0.5 ,連續下坡) ;
//Kendall 的 Tau b 測試斜率是否不為零
MV =多變量(
是的(結腸, colX ) ,
散點圖矩陣( 0 ) ,
相關性多元( 0 ) ,
肯德爾的 Tau ( 1 )
) ;
//繪製數據並添加線性回歸和泰爾斜率
比夫=雙變量(
是的(結腸) ,
X ( colX ) ,
擬合線( {線條顏色( “紅色的” )} ) ,
適合特別(坡(泰爾) , {線條顏色( “綠色的” )} ) ,
無形的
) ;
rbiv =比夫<<報告;
rbiv [ “線性擬合” ] <<像這樣關閉所有;
rbiv [大綱框( 7 )] <<設置標題( “泰爾的斜坡” ) ;
D b =報告( MV )[ “多元” ] <<父母;
D b <<前置( rbiv ) ;
D b [文本框( 4 )] <<設置文本( “泰爾的斜坡” ) ;
rbiv <<關閉窗口;
我不確定 Mann-Kendall 檢驗在多大程度上（如果有的話）與非參數相關的 Kendall Tau 統計量相關。後者可以在多變量平台下訪問：分析>多變量方法，然後從多變量報告中從熱點中選擇非參數相關性。
如果測試完全不同，我會對您對此方法的任何參考感興趣。我正在考慮構建一個插件來實現為水文學領域開發的許多方法（例如赫斯特指數）
問候
戴夫
