交叉验证clmm / clmm2模型

Question

有人有任何建议在R中交叉验证clmm / clmm2（序数包）模型吗？ 我是一个统计新手所以道歉... 我假设通过编写成本函数可以很容易地使用R中的cvTools包。

这是一些虚拟示例代码（我将k设置为低于此示例）。我的数据是序数，使用logit链接和灵活的阈值。在这个例子中，我使用了glm的成本函数，但这当然是错误的。我只是不知道成本函数是什么。注意预测功能只能与clmm2一起使用，而不能与clmm一起使用。

library(ordinal)
library(cvTools)

dat <- subset(soup, as.numeric(as.character(RESP)) <= 24)
dat$RESP <- dat$RESP[drop=TRUE]
head(dat)
m1 <- clmm2(SURENESS ~ PROD, random = RESP, data = dat, link="probit",
        Hess = TRUE, method="ucminf", threshold = "symmetric")

cvFit(m1, data = dat, y = dat$SURENESS,cost= function(r, pi = 0) mean(abs(r-pi) > 0.5), K = 2, R = 10)
warnings()