R version 3.3.2 (2016-10-31) -- "Sincere Pumpkin Patch" Copyright (C) 2016 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > par4 = '0.95' > par3 = '1' > par2 = '0' > par1 = '‎3.184' > par4 <- '0.95' > par3 <- '1' > par2 <- '0' > par1 <- '‎3.184' > #'GNU S' R Code compiled by R2WASP v. 1.2.327 () > #Author: root > #To cite this work: Wessa P., (2013), Normal CDF (v1.0.10), Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_cdfnorm.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > 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)) + } Warning message: NAs introduced by coercion > r [1] NA > ru [1] NA > rb [1] NA > 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 [1] 1.644854 > q1u [1] -1.644854 > q2 [1] 1.959964 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/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) + } Error in if (par1 != "") { : missing value where TRUE/FALSE needed Execution halted