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')
|