Ошибка "Нечисловой аргумент математической функции" при запуске ScheffeTest с объектом aov

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 (формула = x ~ g, drop.unused.levels = TRUE):
неверный тип (список) для переменной 'x'

Как правильно бегать ScheffeTest с aov объект?