更新时间:2023-12-02 13:51:40
library(elrm)
data(drugDat)
drug.elrm <- elrm(formula=recovered/n~sex+treatment,interest=~sex+treatment,r=4,
iter=100000,burnIn=1000,dataset=drugDat)
> summary.elrm(drug.elrm)
Call:
[[1]]
elrm(formula = recovered/n ~ sex + treatment, interest = ~sex +
treatment, r = 4, iter = 1e+05, dataset = drugDat, burnIn = 1000)
Results:
estimate p-value p-value_se mc_size
joint NA 0.14886 0.00173 99000
sex 0.27092 0.69385 0.01204 2649
treatment 0.76739 0.07226 0.00314 13160
95% Confidence Intervals for Parameters
lower upper
sex -0.6217756 1.212499
treatment -0.1216884 1.852346
# If you look at the summary function, it simply outputs formatted results to
# the screen. So instead, we can just work with the original drug.elrm object
names(drug.elrm)
# shows you everything in this object
# to see the p-values
drug.elrm$p.values.se
joint sex treatment
0.001734482 0.012039701 0.003143006
# to get the p-value for joint
drug.elrm$p.values.se[[1]]
# now for the CI
drug.elrm$coeffs.ci
lower upper
sex -0.6217756 1.212499
treatment -0.1216884 1.852346
> drug.elrm$coeffs.ci[[1]]
[1] -0.6217756 -0.1216884
> drug.elrm$coeffs.ci[[1]][1]
[1] -0.6217756
>