我正在尝试计算从我的数据集中的几个不同列(8个单独的字段/列)中的R中的wilcox.test函数计算的伪中值(Hodges-Lehmann估计)。
> head(fantasyproj2)
V1 V2 V3 V4 V5 V6 V7 V8
1 25 25.87 20.35 26.65 27.20 27.0 17.970 27.70
2 24 27.48 19.81 19.57 22.20 27.5 16.350 20.04
3 22 19.89 17.62 21.99 19.12 26.0 22.484 23.70
4 21 17.72 16.09 15.55 18.60 18.5 17.450 14.59
5 21 21.56 17.90 18.46 20.80 23.0 16.540 19.76
6 20 NA 17.73 15.84 19.00 20.5 19.080 15.05
它们都是数字,并且包含一些缺失值(缺失值中没有数据,并且在R中已更改为NA)。当我运行代码时:
fantasyproj$hodges.lehmann <- apply(fantasyproj2[,c(1,2,3,4,5,6,7,8)],1, function(x) wilcox.test(x, conf.int=TRUE, na.action=na.exclude)$estimate)
我收到错误:
Error in uniroot(wdiff, c(mumin, mumax), tol = 1e-04, zq = qnorm(alpha/2, : f() values at end points not of opposite sign
关于这个错误的文献不多except,它说将exact = TRUE
添加到语句中,但这并没有帮助。任何帮助将不胜感激!