r - aov オブジェクトを使用して ScheffeTest を実行すると、「数学関数の引数が数値ではありません」というエラーが発生する

Question

次のコードを実行しました（ここから取得）：

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