| par2 <- as.numeric(par2)par3 <- as.numeric(par3)
 if (par3 <= 0) par3 = 1
 if (par1 != '') {
 par1 <- as.numeric(par1)
 oldpar1 <- par1
 par1 <- (par1 - par2) / par3
 r <- pnorm(par1)
 ru <- pnorm(par1,lower.tail=F)
 rb <- pnorm(abs(par1)) - pnorm(-abs(par1))
 }
 r
 ru
 rb
 if (par4 != '') {
 par4 <- as.numeric(par4)
 q1 <- qnorm(par4,par2,par3)
 q1u <- qnorm(par4,par2,par3,lower.tail=F)
 q2 <- qnorm((1-par4)/2,par2,par3,lower.tail=F)
 }
 q1
 q1u
 q2
 load(file='createtable')
 a<-table.start()
 a<-table.row.start(a)
 a<-table.element(a,'Normal Density Function',2,TRUE)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'mean',header=TRUE)
 a<-table.element(a,par2)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'standard deviation',header=TRUE)
 a<-table.element(a,par3)
 a<-table.row.end(a)
 if (par1 != '') {
 a<-table.element(a,'Area under Normal Density',2,TRUE)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'q value',header=TRUE)
 a<-table.element(a,oldpar1)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'z value',header=TRUE)
 a<-table.element(a,par1)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'P[X <= z]',header=TRUE)
 a<-table.element(a,r)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'P[X > z]',header=TRUE)
 a<-table.element(a,ru)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'P[-|z| <= X <= |z|]',header=TRUE)
 a<-table.element(a,rb)
 a<-table.row.end(a)
 a<-table.row.start(a)
 a<-table.element(a,'P[-|z| <= X] + P[X > |z|]',header=TRUE)
 a<-table.element(a,1-rb)
 a<-table.row.end(a)
 }
 if (par4 != '') {
 a<-table.element(a,'Quantiles',2,TRUE)
 a<-table.row.end(a)
 a<-table.row.start(a)
 dum <- paste('q given P[X <= q] = ',par4,sep='')
 a<-table.element(a,dum,header=TRUE)
 a<-table.element(a,q1)
 a<-table.row.end(a)
 a<-table.row.start(a)
 dum <- paste('q given P[X > q] = ',par4,sep='')
 a<-table.element(a,dum,header=TRUE)
 a<-table.element(a,q1u)
 a<-table.row.end(a)
 a<-table.row.start(a)
 dum <- paste('|q| given P[-|q| <= X] + P[X > |q|] = ',1-par4,sep='')
 a<-table.element(a,dum,header=TRUE)
 a<-table.element(a,q2)
 a<-table.row.end(a)
 }
 a<-table.end(a)
 table.save(a,file='mytable.tab')
 
 |