File:TnormPDF.png

R source code: dtnorm <- function(x,mu=0,sd=1,a=-Inf,b=Inf) { dnorm(x,mu,sd) / (pnorm(b,mu,sd) - pnorm(a,mu,sd)) } x <- seq(-10,10,len=200) y1 <- dtnorm(x,mu=-8,sd=2 ,a=-10,b=10) y2 <- dtnorm(x,mu= 0,sd=2 ,a=-10,b=10) y3 <- dtnorm(x,mu= 9,sd=10,a=-10,b=10) y4 <- dtnorm(x,mu= 0,sd=10,a=-10,b=10)

bitmap(file="tnormPDF.png",type="png256",width=4,height=4,res=300,pointsize=12) par(mar=c(3,3,1,1)) par(mgp=c(1.8,0.7,0)) plot(range(x),range(c(y1,y2,y3,y4)),type="n",xlab="x",ylab="probability density") lines(x,y1,lwd=3,lty=1) lines(x,y2,lwd=3,lty=2,col="blue") lines(x,y3,lwd=3,lty=3,col="red") lines(x,y4,lwd=3,lty=4,col="orange") dev.off