cdf of the standard normal distribution
Never mind. I have figured out a way to do it.
Here it is:
if (fits$link == "probit") {
pct0 <- integrate(dnorm, -Inf, a)
pct0 <- pct0[[1]]}
THanks any way for your answers.
Charlie Liu
Roger Peng
<rpeng at stat.u To: Chunhua Liu/RTP/USEPA/US at EPA
cla.edu> cc: "r-help at stat.math.ethz.ch" <r-help at stat.math.ethz.ch>, Jeff
Gift/RTP/USEPA/US at EPA
02/19/02 Subject: Re: [R] cdf of the standard normal distribution
02:13 PM
Use the pnorm() function.
-roger
_______________________________
UCLA Department of Statistics
rpeng at stat.ucla.edu
http://www.stat.ucla.edu/~rpeng
On Tue, 19 Feb 2002, Liu.Chunhua at epamail.epa.gov wrote:
Dear Experts, I need to calculate the cdf of the standard normal distribution, i.e. H(x) = 1/sqrt(2*pi) integral(exp(-z^2/2) dz), where z is b/w -infi to infi. I know there should be a way to do it in R, but did not know to do it. I'd appreciate any help you could offer. Charlie Liu Graduate student intern at EPA/ECO -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.
-.-.-.-.-
r-help mailing list -- Read
http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To:
r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._.
_._._._._
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._