# NOT RUN {
data("SnedecorCochran")
# }
# NOT RUN {
# analyze using uniform prior:
bma1 <- bayesmeta(y=SnedecorCochran[,"mean"],
sigma=sqrt(SnedecorCochran[,"var"]),
label=SnedecorCochran[,"no"],
tau.prior="uniform")
# analyze using Jeffreys prior:
bma2 <- bayesmeta(y=SnedecorCochran[,"mean"],
sigma=sqrt(SnedecorCochran[,"var"]),
label=SnedecorCochran[,"no"],
tau.prior="Jeffreys")
# compare results:
print(bma1)
print(bma2)
forestplot(bma1)
forestplot(bma2)
# }
Run the code above in your browser using DataLab