我正在运行纵向序数逻辑回归,并试图从我的回归输出中提取特定系数到潜在图形。作为下面的示例,我的输出位于变量 olr1
中。
> summary(olr1)
Call:
mixor(formula = concern_infect_me1 ~ age + sex + rec_date2, data = arr_full,
id = id, link = "logit")
Deviance = 8172207
Log-likelihood = -4071296
RIDGEMAX = 0.5
AIC = -4086121
SBC = -4086182
Estimate Std. Error z value P(>|z|)
(Intercept) -6.3494e+06 1.5792e+01 -4.0206e+05 < 2.2e-16 ***
age -2.8328e+05 9.2873e-01 -3.0502e+05 < 2.2e-16 ***
sex1 1.7476e+07 2.9853e+01 5.8538e+05 < 2.2e-16 ***
sex2 9.4361e+06 1.3226e+02 7.1345e+04 < 2.2e-16 ***
rec_date24 -1.5589e+07 2.7552e+01 -5.6579e+05 < 2.2e-16 ***
rec_date25 1.7540e+06 1.0285e+01 1.7053e+05 < 2.2e-16 ***
rec_date26 9.4236e+06 1.8388e+01 5.1250e+05 < 2.2e-16 ***
rec_date27 3.8766e+07 1.9112e+02 2.0283e+05 < 2.2e-16 ***
rec_date28 5.4125e+06 1.0620e-07 5.0966e+13 < 2.2e-16 ***
rec_date29 -1.2950e+06 3.3328e+01 -3.8857e+04 < 2.2e-16 ***
rec_date210 2.6721e+05 5.3520e+00 4.9928e+04 < 2.2e-16 ***
rec_date211 -6.9163e+03 1.2743e+00 -5.4276e+03 < 2.2e-16 ***
rec_date212 1.0876e+06 3.3343e+01 3.2619e+04 < 2.2e-16 ***
rec_date213 3.8482e+06 6.7283e+00 5.7195e+05 < 2.2e-16 ***
Random.(Intercept) 1.6208e+16 5.4539e+10 2.9718e+05 < 2.2e-16 ***
Threshold2 5.1255e+06 8.4926e+01 6.0353e+04 < 2.2e-16 ***
Threshold3 5.1255e+06 7.8226e+01 6.5521e+04 < 2.2e-16 ***
Threshold4 5.1255e+06 5.1224e+01 1.0006e+05 < 2.2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
我试图将 rec_date24 到 rec_date213 的系数估计理想地作为一个列表,以便我可以将其设置为 y 变量的列表来制作表格,然后制作图表,以便某些变量 {{1} } 将保存值列表:-1.559e+07, 1.7540e+06,..., 3.8482e+06。
我是 R 的相对新手,所以任何简单的解释都将不胜感激。在此先感谢您的帮助!