Question

我在nparcomp :: mctp（）中遇到以下问题。我键入以下命令并产生以下错误消息，我找不到解决方案。

数据：

data <- data.frame(
  Values = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 9, 3, 3, 1, 1, 2, 3, 0, 3, 2, 1, 2, 2, 0, 1, 4, 1, 0, 2, 2, 8, 4, 2, 2, 3, 2, 1, 1, 0, 0, 0, 3, 8, 7, 3, 2, 6, 3, 2, 3, 3, 4, 4, 10, 0, 0, 5, 4, 2, 6, 21, 5, 14, 3, 24, 5, 3, 3, 5, 3, 4, 3, 3, 3, 0, 9, 3, 2, 0, 4, 7, 3, 20, 2, 3, 1, 7, 5, 5, 22, 40, 23, 7, 45, 38, 11, 9, 7, 3, 32, 5, 1, 4, 6, 5, 0, 5, 4, 2),
  Weeks  = sort(rep(1:12, 10)),
  Rep    = rep(1:10, 12)
  )

使用过的功能命令：

npar.1 <- nparcomp::mctp(Values ~ Weeks, 
           data = data, 
           type = "Tukey", alternative = "two.sided", 
           asy.method = "normal")

错误讯息：

Error in checkmvArgs(lower = lower, upper = upper, mean = delta, corr = corr,  : ‘lower’ not specified or contains NA 
In addition: Warning message: 
In cov2cor(CV) : diag(.) had 0 or NA entries; non-finite result is doubtful

我也试过了其他的渐近逼近方法，但是没有成功。

所以，我找到了这条指令（dataset limitation in R package "nparcomp"）并遵循修改源代码的指令。我安装了修改后的源代码如下：

install.package("path", repo = NULL, type = "source")

不幸的是，没有任何改变的结果。

请问，有没有人知道如何解决这个问题或有类似的困难？

非常感谢。

Answer 1

如果我正确理解你的数据，我相信你会看到这一点，因为有一些周，无论他们的代表是什么，只有零值。因此，当 Tukey 检验试图找到一个用于进行比较的大均值时，它会遇到困难。

在不知道您的数据的情况下，我不知道删除第 1 周和第 2 周是否合适，但似乎解决了错误。

nparcomp :: mctp中的错误：“checkmvArgs中的错误...未指定或包含NA”，尽管数据不包含NA

1 个答案: