我有30个研究的数据集,其中我正在使用几个非参数检验来获得p值,以作为研究之间距离的度量,以将研究聚类在一起。
ks测试具有nboots
,Cramer-von具有replicates
,AD具有Nsim
,我相信它们都是引导p值的一种形式。
nestedRanksTest
是一个引导Wilcox测试的程序包,但我似乎无法使其与我的数据配合使用。我使用Wilcox.test
计算研究之间的p值的函数是:
mat <- outer(1:30, 1:30, Vectorize(function(i,j)
{wilcox.test(as.numeric(rep(seq(0,14,1),as.vector(data[i,]))),
as.numeric(rep(seq(0,14,1),as.vector(data[j,]))))$p.value}))
有没有可以做到这一点的包裹?