r - 累積リンク混合モデルによる確率予測

Question

累積リンク混合モデルをordinalパッケージに当てはめようとしていますが、予測確率の取得について理解できないことがあります。ordinalパッケージの次の例を使用します。

   library(ordinal)
data(soup)
## More manageable data set:
dat <- subset(soup, as.numeric(as.character(RESP)) <=  24)
dat$RESP <- dat$RESP[drop=TRUE]
m1 <- clmm2(SURENESS ~ PROD, random = RESP, data = dat, link="logistic",  Hess = TRUE,doFit=T)
summary(m1)
str(dat)

現在、新しいデータセットの確率の予測を取得しようとしています

newdata1=data.frame(PROD=factor(c("Ref", "Ref")), SURENESS=factor(c("6","6")))

と

predict(m1, newdata=newdata1)

しかし、私は次のエラーが発生しています

Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : 
  contrasts can be applied only to factors with 2 or more levels

このエラーが発生するのはなぜですか? 間違った構文に何かありますpredict.clmm2()か？一般に、predict.clmm2() が出力する確率はどれですか? Pr(J<j)またはPr(J=j)? _ 誰かが、カテゴリ (順序) 序数混合モデルを特に R に適合させることに関する情報 (サイト、書籍) 資料を教えてもらえますか?