为什么我得到的p值等于NA

时间:2018-08-06 08:12:36

标签: r

我正在使用energy软件包进行距离协方差检验和距离相关检验。

测试给出的结果是:

> test<-dcov.test(DATA_CATEG[1:nrow(DATA_CATEG),15],
+                 DATA_CATEG[1:nrow(DATA_CATEG),3])
> test

    Specify the number of replicates R (R > 0) to perform the test of independence

data:  index 1, replicates 0
nV^2 = 25.986, p-value = NA
sample estimates:
     dCov 
0.1248535 

是编程问题还是什么?假设我的数据集是干净的,并且没有任何NA值。

这是为了使其可复制:

    > dput(DATA_CATEG[1:20,c(3,15)])
structure(c(2L, 3L, 2L, 1L, 3L, 2L, 2L, 1L, 2L, 1L, 2L, 3L, 2L, 
2L, 1L, 3L, 2L, 2L, 2L, 2L, 3L, 4L, 7L, 3L, 3L, 6L, 3L, 9L, 26L, 
17L, 22L, 13L, 10L, 30L, 4L, 7L, 7L, 6L, 14L, 16L), .Dim = c(20L, 
2L), .Dimnames = list(NULL, c("IMPERFECTIONS", "OUTCOME")))

预先感谢您的帮助

1 个答案:

答案 0 :(得分:0)

您忘记指定R参数。此方法需要它。

您可以在文档中找到有关该主题的更多信息(有关该主题的建议值和一些文献): https://www.rdocumentation.org/packages/energy/versions/1.7-4/topics/dcov.test