<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />こんにちは、みんな。 Mann-Kendall を実行して、一時的な傾向が有意かどうかをテストしようとしています。 JMPのマニュアルやオンラインで見つけることができません。 誰にもアイデアはありますか？
ありがとう！
MSとデイブに感謝します！ Kendall の Tau と Mann-Kendall 検定が同じ (および/または関連) であるかどうかもわかりませんでした。 多変量テストに関する情報と、スクリプトについてもありがとうございます。
定義 Mann-Kendall test statistic と Kendall Tau b はしばしば異なる用語で記述されますが (それぞれ時系列と相関の観点から)、私が知る限り、それらは本質的に同じです。
以下は、私がかなり前に書いたスクリプトです。 Theil の勾配を計算して表示し、Mann-Kendall 検定を実行します。 比較のために、パラメトリック線形回帰直線がプロットに表示されます。 「By」関数と、季節的なパターンを説明する可能性を含めることで、スクリプトを一般化することを計画しました。 しかし、私はまだそれらの計画を実現する時間がありませんでした. これがあなたにいくつかのアイデアを与えることを願っています。
dt =現在のデータ テーブル() ;
dlg =列ダイアログ(
「データ テーブルから Theil 勾配と Kendall Tau を推定する」 ||
シャア( dt <<名前を取得) 、
イース=列リスト( 「Yを選択」 、
マックスコル( 1 ) 、
ミンコル( 1 ) 、
データ・タイプ(数値)
) 、
Xs =列リスト( 「×を選択」 、
マックスコル( 1 ) 、
ミンコル( 1 ) 、
データ・タイプ(数値)
) 、
グループ=列リスト( 「（動作していない）によって」 、マックスコル( 1 ) )
) ;
列=桁( dt 、 dlg [ 「イース」 ] ) ;
列X =桁( dt 、 dlg [ 「×」 ] ) ;
colGroup =桁( dt 、 dlg [ "グループ" ] ) ;
// Theil 勾配の計算
n = N 行( dt ) ;
斜面= J ( n 、 n 、 . ) ;
為に(私= 1 、私< n 、私++、
為に( j =私+ 1 、 j <= n 、 j ++、
斜面[私、 j ] = (列[ j ] -列[私])) / (列X [ j ] -列X [私]))
)
) ;
彼ら=分位点( 0.5 、斜面) ;
//Kendall の Tau b を使用して、勾配がゼロと異なるかどうかをテストします
mv =多変量(
よ(列、列X ) 、
散布図行列( 0 ) 、
相関多変量( 0 ) 、
ケンダルのタウ( 1 )
) ;
// データをプロットし、線形回帰と Theil の勾配を追加します
ビブ=二変量(
よ(列) 、
バツ(列X ) 、
適合線( {線の色( "赤" )} ) 、
フィットスペシャル(スロープ(彼ら) 、 {線の色( "緑" )} ) 、
見えない
) ;
rbiv =ビブ<<報告;
rbiv [ 「リニアフィット」 ] <<このようにすべて閉じます;
rbiv [アウトラインボックス( 7 )] <<タイトルを設定( 「テイルズスロープ」 ) ;
デシベル=報告( mv )[ 「多変量」 ] <<親;
デシベル<<先頭に追加( rbiv ) ;
デシベル[テキストボックス( 4 )] <<テキストを設定( 「テイルズスロープ」 ) ;
rbiv <<ウィンドウを閉じる;
Mann-Kendall 検定が、ノンパラメトリック相関の Kendall Tau 統計にどの程度関連しているか (もしあったとしても) はわかりません。後者は、多変量プラットフォームでアクセスできます: [分析] > [多変量メソッド] から、多変量レポートから、ホットスポットから [ノンパラメトリック相関] を選択します。
テストが完全に異なる場合は、この方法に関する参考文献に興味があります。水文学の分野向けに開発された多くの方法 (Hurst Index など) を実装するアドインの構築を検討しています。
よろしく
デイブ
