1

有没有人有任何建议来交叉验证 R 中的 clmm/clmm2(序数包)模型?我是一个统计新手,所以很抱歉......我认为通过编写成本函数可以很容易地使用 R 中的 cvTools 包。

这是一些虚拟示例代码(我仅为此示例设置了低 k)。我的数据是有序的,使用logit链接,并且阈值灵活。在这个例子中,我使用了 glm 的成本函数,但这当然是错误的。我只是不知道成本函数是什么。注意 predict 函数只能用于 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()
4

0 回答 0