Pettitt - R 中的显着性水平
Pettitt - significance levels in R
我正在尝试在 R 上实现 Pettitt 测试。我自己编码，但参考了 R 的内置函数 pettitt.test。(https://www.rdocumentation.org/packages/trend/versions/1.1.4/topics/pettitt.test)
我可以成功实现它，但我的问题是关于 p 值的。我正在用这个等式计算 p 值：
p=2e^((-6U ̂^2)/(n^3+n^2 ))
U 是检验统计量的值，n 是我的样本长度。 （我已经检查过，我得到的结果与内置函数相同pettitt.test）
对于 5% 的置信区间，如果此 p 值低于 0.5，备择假设应该为真？或 0.05？我对 R 的文档（之前共享的网络链接）表示低于或等于 0.5 感到困惑，而 Pettitt 的论文（1979 年）第 5 页也表示低于或等于 0.5
对澄清这个问题有什么帮助吗？
论文说的是
... where the approximation holds good, accurate to two decimal places, for p(OA)<0.5
（强调）。
也就是说，这与要拒绝原假设的临界水平无关（通常 p=0.05），而是关于近似值适用的值范围。
大多数 R 函数和封装只是 return p 值，让用户决定他们想做什么来解释它（即拒绝或未能拒绝给定 alpha 水平的原假设，将其视为反对零假设等的连续（“Fisherian”）证据强度。
如果您希望，如果 p 的计算值 >0.5（例如“p 值的近似值可能不可靠”），您可以让您的函数发出警告
Pettitt, A. N.“变点问题的非参数方法。”皇家统计学会杂志。 C 系列（应用统计）28，没有。 2 (1979): 126–35。 https://doi.org/10.2307/2346729.
我正在尝试在 R 上实现 Pettitt 测试。我自己编码，但参考了 R 的内置函数 pettitt.test。(https://www.rdocumentation.org/packages/trend/versions/1.1.4/topics/pettitt.test)
我可以成功实现它，但我的问题是关于 p 值的。我正在用这个等式计算 p 值：
p=2e^((-6U ̂^2)/(n^3+n^2 ))
U 是检验统计量的值，n 是我的样本长度。 （我已经检查过，我得到的结果与内置函数相同pettitt.test）
对于 5% 的置信区间，如果此 p 值低于 0.5，备择假设应该为真？或 0.05？我对 R 的文档（之前共享的网络链接）表示低于或等于 0.5 感到困惑，而 Pettitt 的论文（1979 年）第 5 页也表示低于或等于 0.5
对澄清这个问题有什么帮助吗？
论文说的是
... where the approximation holds good, accurate to two decimal places, for p(OA)<0.5
（强调）。 也就是说，这与要拒绝原假设的临界水平无关（通常 p=0.05），而是关于近似值适用的值范围。
大多数 R 函数和封装只是 return p 值，让用户决定他们想做什么来解释它（即拒绝或未能拒绝给定 alpha 水平的原假设，将其视为反对零假设等的连续（“Fisherian”）证据强度。
如果您希望，如果 p 的计算值 >0.5（例如“p 值的近似值可能不可靠”），您可以让您的函数发出警告
Pettitt, A. N.“变点问题的非参数方法。”皇家统计学会杂志。 C 系列（应用统计）28，没有。 2 (1979): 126–35。 https://doi.org/10.2307/2346729.