y=3.2 mu=c(2,2.5,3,3.5,4) sigma=1 z=(y-mu)/sigma z prior=rep(1/5,5) prior lik=dnorm(z) lik prior*lik posterior=lik*prior/sum(lik*prior) posterior lik=dnorm(y,mean=mu,sd=sigma) lik ?dnorm mu=seq(4,9.5,by=0.5) mu prior=rep(1/12,12) prior qnorm(0.025,2,1)