[Rd] ppoints (PR#7538)
Prof Brian Ripley
ripley at stats.ox.ac.uk
Mon Jan 24 10:37:44 CET 2005
On Wed, 19 Jan 2005 tobias.verbeke at telenet.be wrote:
> Dear r-bugs,
>
> Whilst playing with ppoints I discovered
> that when one uses it directly, occasional
> NA's in a vector also become data fractions:
>
> ppoints(c(1,2,NA,4))
>
> Would it be a good idea to add a warning message
> as in:
>
> ppoints <- function (n, a = ifelse(n <= 10, 3/8, 1/2))
> {
> if(any(is.na(n))) warning("'n' contains NA's")
> if(length(n) > 1) n <- length(n)
> if(n > 0)
> (1:n - a)/(n + 1-2*a)
> else numeric(0)
> }
Why? There are 4 points in your vector, and the result is perfectly
valid as documented, even if they were all NAs.
> Another minor remark concerning ?ppoints. It says:
>
> n: either the number of points generate or a vector of
> observations. ^^^^^
As you see, that does not line up, but the typo has been fixed.
--
Brian D. Ripley, ripley at stats.ox.ac.uk
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UK Fax: +44 1865 272595
More information about the R-devel
mailing list