data("SnedecorCochran")
if (FALSE) {
# 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