2

我有一个如下的数据集

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 并看不到任何东西。 如何将表格限制为某些列或仅针对因子？它还为整数类列的每个值创建了一个单独的表，这是有这么多列的一个重要原因。

4

1 回答 1

1

您可以确定哪些列是因素，然后只将这些列传递给table：

w <- sapply( data.AMIDS_, function(x) is.factor(x) )
table( data.AMIDS_d[, w] )
于 2013-05-03T19:18:50.280 回答