# 2/16/09
set.seed(13715)
X=rnorm(50);qqnorm(X);qqline(X)
?shapiro.test
shapiro.test(X)
?rt
args(rt) # print out the argument list
X=rt(n=50,df=5);qqnorm(X);qqline(X)

theta=seq(-10,10,.01)
rho1=-theta/(1+theta^2)
plot(y=rho1,x=theta,type='l')
abline(h=0);abline(v=0)
abline(h=-1/2,col='blue')
abline(h=1/2,col='blue')
abline(v=1,col='green')
abline(v=-1,col='green')