0

我有一个数据框，我对其进行了双向方差分析。

dput(m3)

结构（列表（Delta = c（-40，-40，-40，-40，-31.7，-29.3，-27.8，-26.7，-26.2，-25.4，-24.7，-23.1，-23，-22.9， -22.4，-22.2，-21.4，-21，-20.8，-15.1，-14.9，-14.1，-6.2，-6.2，-6，-5.3，-4.9），位置=结构（c（1L，1L， 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 3L, 2L, 3L, 3L, 3L) , .Label = c("int", "pen + int", "ter + pen"), class = "factor"), Between = c(0L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 0L, 2L, 1L, 0L, 1L, 0L, 2L, 0L, 2L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L ), 相对 = 结构 (c(5L, 6L, 6L , 7L, 8L, 3L, 3L, 4L, 5L, 4L, 3L, 5L, 3L, 5L, 7L, 5L, 4L, 6L, 3L, 3L, 6L, 2L, 1L, 2L, 1L, 1L, 1L), .Label = c("1&2", "2&3", "2&4", "2&5", "3&4",“3&5”、“3&6”、“4&6”）、class =“因子”））、class =“data.frame”、row.names = c（NA，-27L））

library(agricolae)
aov.2sum=aov(Delta.~Location*X.between, data=m3)

我想使用HSD.test分析数据，就像使用相同功能的另一个数据帧一样。我遵循包手册中的代码格式，如下所示。

tx <- with(m3, interaction(Location, X.between))
amod <-aov(Delta~tx, data=m3)
test=HSD.test(amod, "tx", group=TRUE)

然后我收到以下错误

(x, value = value)中的错误.rowNamesDF<-：不允许重复的“row.names”另外：警告消息：设置“row.names”时的非唯一值：“int.0”、“pen + int.1” , '笔 + int.2', 'te + int.0', 'te + int.1'

经过进一步分析，我发现我的重复行名错误与我的 X.between 功能有关。当我使用以下代码时，我得到相同的重复行名错误：

HSD.test(amod, "X.between", group=TRUE)
>> Error in data.frame(row.names = means[, 1], means[, 2:6]) : 
  duplicate row.names: 0, 1, 2

如何为 HSD.test 选择行名？那我怎样才能改变我的行名呢？或者只是避免这种重复错误？

感谢您的所有帮助。

4

0 回答 0