我正在使用包“betareg”(版本 3.0-3)在 R 中进行 beta 回归。对于蒙特卡洛模拟,我想将分散 - 参数“phi”建模为 logit 模型。不幸的是,该软件包不支持 phi 的链接功能“logit”。我已经看到可以提供一个自写的“glm - 链接”对象,但我不知道如何做到这一点。有谁知道,我如何解决我的问题并将“phi”建模为逻辑模型。
betareg(formula, data, subset, na.action, weights, offset,
link = c("logit", "probit","cloglog", "cauchit", "log", "loglog"),
link.phi = NULL, type = c("ML", "BC", "BR"),
control = betareg.control(...), model = TRUE,y = TRUE, x = FALSE, ...)
link.phi:
精度模型(phi)中链接函数的字符规范。目前支持“identity”、“log”、“sqrt”。默认值为“log”,除非公式的类型为 y ~ x,其中默认值为“identity”(为了向后兼容)。或者,可以提供类“link-glm”的对象。