使用aov对象运行ScheffeTest时出现错误'数学函数的非数字参数'

Question

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

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.default中的错误（formula = x~g，drop.unused.levels = TRUE）：
        变量'x'

的类型（列表）无效

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