包含零向量的R beanplot比较

Question

我希望能够使用beanplot来显示分类预测变量中4个增量的可视比较，如绘制在响应变量上。我的问题是，其中一个增量是一个充满零的向量，所以可以理解，beanplot例程的shapiro.test阶段将不起作用。

（这是我的错误）

Error in shapiro.test(x) : all 'x' values are identical 

有没有一种方法，我可以绕过这一点，并说，如果有零向量则该增量什么也不做，进入下一个让我在x轴上最多有4个增量，而只有向量包含值大于0的增量才有3个beanplots？

非常感谢您的帮助。

一切顺利，

贝克斯

Answer 1

我猜问题是，当一个类别的var是0。你可以测试这种情况下，并为该组NA设置的值。这应该与beanplot一起工作。首先，一些测试数据

dd<-data.frame(
    cat =sample(letters[1:5], 50, replace=T), 
    val = runif(50) 
) 
dd$val[dd$cat=="c"]<-0 

而现在的数据修复和绘图代码

library(beanplot) 

catvar <- with(dd, tapply(val, cat, var)) 
if(any(catvar==0)) { 
    dd$val[dd$cat %in% names(catvar[catvar==0])] <- NA 
} 

beanplot(val~cat, dd) 

和产生的情节，如