捕获创建仅在R中具有一个非空行或列的表的数据

Question

我正在对一堆变量执行chisq.test ，但它意外失败：

Error in chisq.test(var1, var2) : 'x' and 'y' must have at least 2 levels

快速检查得出：

> table(var1,var2)
    var2
var1    1    2    3    4    5    6
   0   34   48  153  710 1667 2177
   1    0    0    0    0    0    0

因此，列联表中有一行为空，只剩下一行。 如何测试此类变量并在循环中跳过它们？

Answer 1

您可以使用nlevels来获取一个因子中的级别数。 例如，您可以执行以下操作：

if (nlevels(var1) > 1 & nlevels(var2) > 1) {
    chisq.test(var1, var2)
} else {
    warning('Not enough levels to perform chi-sq test.')
}