par(lwd=1.5,cex=1.5,mar=c(5,4,1,1)) mu <- 0; h0 <- 1/2; priorsd <- sqrt(1/h0) ybar <- 1.96; sigma <- 10 n <- 100; hs <- 100/(10^2); likesd <- sqrt(1/hs) hpost <- h0 + hs postsd <- sqrt(1/hpost) postmean <- 1.96/(3/2) # = 1.31 x<-seq(-5,5,length=101) plot(x,dnorm(x,postmean,postsd),type="l", xlab=as.expression(as.name('theta')),ylab="density",cex=1.5) lines(x,dnorm(x,mu,priorsd),lty=2) lines(x,dnorm(x,ybar,likesd),lty=3) legend("topleft", c("posterior","prior","likelihood"),lty=1:3) # legend(-4, .4, c("posterior","prior","likelihood"),lty=1:3) abline(v=postmean)