<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 统计量相关。后者可以在多变量平台下访问：分析>多变量方法，然后从多变量报告中从热点中选择非参数相关性。
如果测试完全不同，我会对您对此方法的任何参考感兴趣。我正在考虑构建一个插件来实现为水文学领域开发的许多方法（例如赫斯特指数）
问候
戴夫
