命令:library("dosresmeta")
data("cc_ex")
fix("cc_ex")
modl<-dosresmeta(formula = logrr ~ dose, type ="cc", cases = case, n = n, lb=lb,ub=ub,data = cc_ex)
summary(modl)
predict(modl,delta=1)
predict(modl,delta=5)
结果:
> predict(modl,delta=1)
delta pred ci.lb ci.ub
1 0.04542885 0.004928432 0.08592926
> predict(modl,delta=5)
delta pred ci.lb ci.ub
5 0.2271442 0.02464216 0.4296463
书中结果:
> predict(modl,delta=1)
delta pred ci.lb ci.ub
1 1.0465 1.0050 1.0897
> predict(modl,delta=5)
delta pred ci.lb ci.ub
5 1.2550 1.0249 1.5367


雷达卡



京公网安备 11010802022788号







