如何计算r中Williams趋势测试的p值（SAS中“PROBMC”的类似函数）

我想计算R中William trend test的p值，因为我已经知道了t统计量。 在SAS中，我可以使用函数PROBMC，如下所示

 PROBMC(distribution, q, prob, df, nparms<, parameters>)

---

R中是否有类似的函数来执行此操作？

我认为您可能运气不好

使用

库（“sos”）；查找DFN（“Williams趋势分布”）

并搜索结果，找到两个包，

PMCMCR

和

StatCharrms

，它们具有执行Williams测试的功能，但看起来这些只使用论文中的列表值来获得p=0.05的临界值，而不是直接计算分布/p值

获得完整分布/p值的计算看起来相当复杂，因此不太可能有人决定在R中实现它

如Williams（1971）所述（见参考文献），完整的计算非常冗长，分三个阶段进行

这将是一个很好的计算统计项目的人

       if parameters t=2.6,  k = 6, [nu] = 42, and t = 2.60 then probability is .9924467341.
       using (prob=probmc("williams",2.6,.,42,6);)