1

我运行了以下代码（取自此处）：

set.seed(123)
Njk   <- 10
P     <- 2
Q     <- 2
R     <- 3
DV_t1 <- rnorm(P*Q*Njk, -3, 2)
DV_t2 <- rnorm(P*Q*Njk,  1, 2)
DV_t3 <- rnorm(P*Q*Njk,  2, 2)
dfSPFpq.rL <- data.frame(id=factor(rep(1:(P*Q*Njk), times=R)),
                         IVbtw1=factor(rep(1:P, times=Q*R*Njk)),
                         IVbtw2=factor(rep(rep(1:Q, each=P*Njk), times=R)),
                         IVwth=factor(rep(1:R, each=P*Q*Njk)),
                         DV=c(DV_t1, DV_t2, DV_t3))
aovSPFpq.r <- aov(DV ~ IVbtw1*IVbtw2*IVwth + Error(id/IVwth), data=dfSPFpq.rL)
summary(aovSPFpq.r)

现在我想运行ScheffeTest：

我试过了：

 library(DescTools)
 ScheffeTest.aov(aovSPFpq.r)

pf(psi^2/(MSE * sscoeff * dfgrp), df1 = dfgrp, df2 = dferr, lower.tail = FALSE) 中的错误：
数学函数 ScheffeTest(aovSPFpq.r) 的非数字参数
model.frame 中的错误。默认（公式 = x ~ g，drop.unused.levels = TRUE）：
变量“x”的类型（列表）无效

ScheffeTest使用aov对象运行的正确方法是什么？

4

0 回答 0