curve(dnorm(x,0,1),from=-3,to=4, col=2, xlab='',ylab='',xaxt='n',yaxt='n') curve(dnorm(x,1,1),add=T,col=1) segments(0,dnorm(0,0,1),0,0,lty=2,col=2) segments(1,dnorm(1,1,1),1,0,lty=2,col=1) axis(1,at=1, label=expression(theta)) text(1.5,dnorm(1.5,1,1),expression(hat(theta)[1]),pos=4) text(-0.5,dnorm(-0.5,0,1),expression(hat(theta)[2]),pos=2, col=2) mtext('Density',side=2,line=0.75)