我正在尝试用R执行惩罚多项逻辑回归。我使用pmlr
包,这非常好。不幸的是,此代码不估计伪R2。为了自己估算它,我需要使用plmr
执行“仅拦截”模型。我使用以下命令执行此操作:
Res <- pmlr(y~1, data=x, weights=NULL,
alpha=0.05, penalized = TRUE,
method = c("likelihood","wald")[1], joint=FALSE)
使用函数y~1
只应合并拦截,对吧?
不幸的是,我收到以下错误消息:
错误1:m:长度为0的参数
有人可以告诉我如何只能通过拦截执行回归分析吗?在此先感谢!!!
亲切的回归 S Mau