如何获得MICE之后的赔率比率,然后是pool()

时间:2017-03-23 14:46:16

标签: r

我想进行逻辑回归分析,根据儿童创伤的类型预测结果的几率(即自杀 - 这是二元结果(0/1))(即性行为 - 这是报告了二进制变量(0/1)。

我的数据集中有一些缺失值(即data.noid),因此我使用MICE使用以下命令:

set.seed(123)
imp20 = mice(data.noid, m=20, printFlag=FALSE, maxit = 50, 
                                               seed=2525, control = list(maxiter = 500))

它对我的数据运作良好。所以我继续进行glm分析,然后进行pool()round()函数,以获得我的结果(即自杀)和童年创伤(即童年性虐待)之间的关联。我符合以下条件:

fit.mi20 = with(data=imp20, exp = glm(suicide~ sexabuse , family = binomial))
combFit = pool(fit.mi20) 
round(summary(combFit),2)

Results: 
                   est   se     t     df Pr(>|t|) lo 95 hi 95 nmis  fmi lambda
(Intercept)       -2.01 0.32 -6.24 130.08     0.00 -2.65 -1.37   NA 0.14   0.13
sexabuse2          0.50 0.53  0.95 153.93     0.35 -0.54  1.54   NA 0.07   0.06

我希望你能告诉我如何获得OR(而不是现在提出的'est')。我将不胜感激任何建议和指导。

0 个答案:

没有答案