包含零向量的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) 

和產生的情節，如