计算许多列的伪中位数

时间:2015-12-15 16:26:43

标签: r

我正在尝试计算从我的数据集中的几个不同列(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添加到语句中,但这并没有帮助。任何帮助将不胜感激!

0 个答案:

没有答案