[R] problems with rexp ?

[Setzer.Woodrow@epamail.epa.gov]

instead of (1 - pexp(data)), try using pexp(data, lower.tail=FALSE), for
a more accurate calculation of upper tail probabilities.

R. Woodrow Setzer, Jr.                                            Phone:
(919) 541-0128
Experimental Toxicology Division                       Fax:  (919)
541-5394
Pharmacokinetics Branch
NHEERL MD-74; US EPA; RTP, NC 27711


|---------+------------------------------>
|         |           Juan Ramon Gonzalez|
|         |           <jrgonzalez at ico.scs|
|         |           .es>               |
|         |           Sent by:           |
|         |           owner-r-help at stat.m|
|         |           ath.ethz.ch        |
|         |                              |
|         |                              |
|         |           05/03/02 09:33 AM  |
|         |                              |
|---------+------------------------------>
  >---------------------------------------------------------------------------------------------------------------------|
  |                                                                                                                     |
  |        To:      r-help at stat.math.ethz.ch                                                                            |
  |        cc:                                                                                                          |
  |        Subject: [R] problems with rexp ?                                                                            |
  >---------------------------------------------------------------------------------------------------------------------|




Does anyone know if R have any problems with the exponential random
number
generation (function rexp)?

I comment it because I executed

data<-sort(rexp(100))
plot(data,dexp(data)/(1-pexp(data)),type="l")

and the graphic isn't constant.

(Note: exponential distribution have a constant hazard failure rate).

Thank you,

Juan

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.
-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._