Hi All,
I recently discovered an interesting issue with the punif() function. Let
X~Uiform[a,b] then the CDF is defined by F(x)=(x-a)/(b-a) for (a<= x<= b).
The important fact here is the domain of the random variable X. Having said
that, R returns CDF for any value in the real domain.
I understand that one can justify this by extending the domain of X and
assigning zero probabilities to the values outside the domain. However,
theoretically, it is not true to return a value for the CDF outside the
domain. Then I propose a patch to R function punif() to return an error in
this situations.
Example:
> punif(10^10)
[1] 1
Regards,
Hamed.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
