[R] problems with rexp ?

Liaw, Andy andy_liaw at merck.com
Fri May 3 16:08:28 CEST 2002


Have you looked at the range of the numbers?  I got

> range(dexp(data)/(1-pexp(data)))
[1] 1 1

and 

range(dexp(data)/(1-pexp(data)))-1.0
[1] -1.665335e-15  1.110223e-15

Andy


> -----Original Message-----
> From: Juan Ramon Gonzalez [mailto:jrgonzalez at ico.scs.es]
> Sent: Friday, May 03, 2002 9:34 AM
> To: r-help at stat.math.ethz.ch
> 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
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._.
> _._._._._._._._._
> 


------------------------------------------------------------------------------
Notice:  This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (Whitehouse Station, New Jersey, USA) that may be confidential, proprietary copyrighted and/or legally privileged, and is intended solely for the use of the individual or entity named in this message.  If you are not the intended recipient, and have received this message in error, please immediately return this by e-mail and then delete it.

==============================================================================

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



More information about the R-help mailing list