lc<-seq(0.9,0.99,0.005)
alpha<-1-lc

x<-1-alpha/2
quantnorm<-qnorm(x)
tquant14<-qt(x,14)
tquant28<-qt(x,28)
t<-cbind(lc,alpha,x,tquant14,tquant28,quantnorm)
colnames(t)<-c("lc","alpha","p","t14","t24","normstd")

print(t)