post2 dens(inv_logit(post2$a[,1,1]),lwd=3,col=4,xli
m=c(0,0.5),xlab="probability admission")
dens(inv_logit(post2$a[,2,1]),lwd=3,col=4,lty
=2,add=TRUE)
dens(inv_logit(post2$a[,1,2]),lwd=3,col=2,add
=TRUE)
dens(inv_logit(post2$a[,2,2]),lwd=3,col=2,lty
=2,add=TRUE)
> precis(m1,depth=3)
mean sd 5.5% 94.5% n_eff Rhat4
a[1] -1.82 0.09 -1.96 -1.67 1483 1
a[2] -0.89 0.07 -1.00 -0.78 1562 1
total e ect shows disadvantage
direct e ect confounded
D2,G2
D2,G1
0.0 0.1 0.2 0.3 0.4 0.5
0 10 20 30 40
probability admission
Density
D1,G1
D1,G2