R中的Wilcox.test - 没有足够的x观察答案

【问题标题】：Wilcox.test in R - not enough x observationsR中的Wilcox.test - 没有足够的x观察
【发布时间】：2019-03-12 17:03:27
【问题描述】：

我是 R 新手，正在尝试对水质数据执行 wilcox.test。我的数据很长，我对数据进行了子集化以创建“preWWTF”和“postWWTF”，以便在上游废水处理设施升级之前或之后对数据进行分组。我使用的代码是：

wilcox.test(x=preWWTF$Result [preWWTF$Loc_Analyte=="BarkTop_DP"],
        y=postWWTF$Result [postWWTF$Loc_Analyte=="BarkTop_DP"],
        paired = FALSE)

我收到错误“没有足够的（有限）'x' 观察”。没有 NA 或空白值。然而，preWWTF 的观测值少于 postWWTF。是否有一些代码语言可以用来“截断”发布的 WWTF 数据，以便在前后有相同数量的观察？我认为这是导致问题的原因。谢谢。

【问题讨论】：

这个问题是不可重现的，所以我们都只是猜测。您是否检查以确保您的两列都是数字（而不仅仅是编码为字符或因子的数字）？
我认为<- 是一个错误，它正在做你不打算做的奇怪事情。如果您想更改 preWWTF$Loc_anlayte 的内容，请在不同的行上执行此操作，以便检查函数中的内容。
尝试使用[preWWTF$Loc_Analyte == "BarkTop_DP"]
更好的做法是在函数调用上方的行中创建输入，然后您可以在调用函数之前检查它们是否是您认为的。
Seth 和 Dave2e - 使用 == 而不是

标签： r statistics

【解决方案1】：

我只是在这里写一些可能有助于调试您的问题的东西。

将您的数据分配给名为 x 的新对象

x=preWWTF$Result [preWWTF$Loc_Analyte=="BarkTop_DP"]

检查 x

summary(x)
View(x)

将您在此处看到的内容与您在源 data.frame preWWTF 中看到的内容进行比较，确保您提取了您期望的值。

将你的 y 分配给一个对象

y= postWWTF$Result [postWWTF$Loc_Analyte=="BarkTop_DP"]

检查你

summary(y)
View(y)

如果这些摘要中有任何您不理解的内容，请在此处发布。只要一切都是数字并且没有任何内容是 NA 或 INF，然后继续运行您的模型。

wilcox.test(x, y,
         paired = FALSE)

【讨论】：

谢谢！这发现有我看不到的空数据行（本质上是 NA）。感谢您花时间帮助我。