如何解决解决在R中使用stargazer创建logit回归结果表的错误
我希望展示一个漂亮的Logistic回归表,该表使用观星仪测量高血压,其中包括系数,标准误和显着性(以星号表示)。当我尝试插入占星师的规格时,我看到以下错误消息:“%错误:无法识别的对象类型。”我提供了一些示例数据/下面运行的代码。如何解决?谢谢!
library(stargazer)
library(mfx)
structure(list(AGE = c(40L,23L,24L,18L,30L,33L,32L,63L,22L,24L),IMMIGRANT = c(0,1,1),FAMSIZE = c(2L,2L,3L,6L,1L,1L),HLTH_INS = c(1,0),HYPERTEN = c(0,SMOKE = c(0,PSU = c(2L,2L)),row.names = c(NA,-10L),class = "data.frame")
#The regression works without adjusting for clustered SE
logit<-logitmfx(HYPERTEN~AGE+IMMIGRANT+FAMSIZE+HLTH_INS+
SMOKE,data=sample,atmean=TRUE,robust=T)
logit_mfx_coef <- logit$mfxest[,1]
logit_mfx_se <- logit$mfxest[,2]
stargazer(logit,type="text",title = "Predicting Probability of Hypertension",intercept.bottom=FALSE,coef = logit_mfx_coef,se = logit_mfx_coef,column.labels="Logit mfx",digits=4,align=TRUE)
解决方法
将logit$fit
作为第一个参数传递给stargazer()
。
logitmfx()
操作返回了一堆东西,但是stargazer()
希望将拟合的模型对象(或数据框)作为其第一个参数。
stargazer(logit$fit,type="text",title = "Predicting Probability of Hypertension",intercept.bottom=FALSE,coef = logit_mfx_coef,se = logit_mfx_coef,column.labels="Logit mfx",digits=4,align=TRUE)
输出:
Predicting Probability of Hypertension
=============================================
Dependent variable:
---------------------------
HYPERTEN
Logit mfx
---------------------------------------------
Constant 0.0000
(0.0000)
AGE
IMMIGRANT
FAMSIZE
HLTH_INS
SMOKE
---------------------------------------------
Observations 10
Log Likelihood -0.0000
Akaike Inf. Crit. 10.0000
=============================================
Note: *p<0.1; **p<0.05; ***p<0.01
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。