限制表因素

我有一個數據集像下面限制表因素

particip group device width length accep thresh rating d-rating 
1   RA  Dingo nom  nom  Y  5  8  3 
1   RA  Dingo nom  long Y  4  6  2 
1   RA  Dingo fat  nom  Y  4  6  2 
1   RA  Dingo fat  long N  6  4  -2

，我在其上運行的方差分析，像這樣

aov.AMIDS_d <- aov(d.rating ~ group*device*width*length + Error(particip/(device*width*length))+group,data.AMIDS_d)

該工程確定，直到我嘗試打印條件的手段像這樣

print(model.tables(aov.AMIDS_d,"means"),digits=3)

和它說

Error in model.tables.aovlist(aov.AMIDS_d, "means") : 
design is unbalanced so cannot proceed

根據設計，它應該是平衡的，所以我需要檢查我的數據結構。我試過

table(data.AMIDS_d)

但是這給了我很多的表我突破了max.print並且看不到東西。 如何限制某些列的表格或僅限於某些因素？它正在爲整型列的每個值做一個單獨的表，這是爲什麼有這麼多的一個重要原因。

來源

2013-05-03 Krysta

您可以確定哪些列的因素，那麼只有通過那些列table：

w <- sapply(data.AMIDS_, function(x) is.factor(x)) 
table(data.AMIDS_d[, w])

來源

2013-05-03 19:18:50

謝謝！這工作得很好。但它表示在所有單元中都有相同數量的觀察值，所以我必須找到另一種方法來找出不平衡來自哪裏。。。 – Krysta 2013-05-03 19:31:15

回答

相關問題