# (C) Petteri Pajunen 2006 # X11(width=5,height=5) t <- seq(3,10,by=0.01) dfree <- 6 x1 <- dnorm(t,mean=0,sd=sqrt(dfree/(dfree-2))) # this makes the variances equal x2 <- dt(t,df=dfree) plot(t,x1,type="l",lwd=2,col="red",main="RED: NORMAL, BLUE: STUDENT-T") lines(t,x2,type="l",lwd=2,col="blue") scan(quiet=TRUE) t <- seq(3,10,by=0.1) dfree <- 6 x1 <- dnorm(t,mean=0,sd=sqrt(dfree/(dfree-2))) # this makes the variances equal x2 <- dt(t,df=dfree) plot(t,x1/x2,type="l",lwd=2,main="DENSITY RATIO: NORMAL/STUDENT-T") scan(quiet=TRUE) dev.off()