我正在使用Logistf库使用R使用Firth Logistic回归基于AnAttr对Neilsen评级进行建模。
通常,我在SAS中使用Pearson's Hosmer和Lemeshow测试模型拟合。我仅在查看logistf软件包的内容和文档的输出中看到一个似然比测试,找不到等效的测试。
可能性比检验是否足以验证模型(假设重要的预测变量)?有没有办法从模型输出预测值?这是由“数据输出”处理的吗?如果是,数据将发送到哪里?
代码
fit1 <-logistf(data=DiffScores, RatingMedStan ~ AnAttr, dataout=TRUE)
summary(fit1)
Output
Model fitted by Penalized ML
Confidence intervals and p-values by Profile Likelihood Profile Likelihood
coef se(coef) lower 0.95 upper 0.95 Chisq p
(Intercept) 0.2125129 0.6871306 -1.08186402 1.638132 0.1080064 0.74242534
Age 0.5086712 0.3088479 0.02466206 1.277481 4.2991303 0.03813187
Likelihood ratio test=4.29913 on 1 df, p=0.03813187, n=12
Wald test = 2.712595 on 1 df, p = 0.09955889