[Rd] variance of a scalar (PR#546)
Prof Brian Ripley
Prof Brian Ripley <ripley@stats.ox.ac.uk>
Fri, 19 May 2000 10:04:09 +0100 (BST)
> Date: Fri, 19 May 2000 09:57:10 +0100 (BST)
> From: Jonathan Rougier <J.C.Rougier@durham.ac.uk>
> To: Prof Brian Ripley <ripley@stats.ox.ac.uk>
> cc: r-devel@stat.math.ethz.ch, R-bugs@biostat.ku.dk
> Subject: Re: [Rd] variance of a scalar (PR#546)
>
> On Fri, 19 May 2000, Prof Brian Ripley wrote:
>
> > > I was surprised to find that the variance of a scalar, using
> > > var(), is NA. Surely this should be zero?
> > >
> >
> > Do you mean a vector of length 1? (S has no scalars.) No, it should be
> > NA, as sum((x-xbar)^2)/(n-1) is NaN, or, in statistical terms, one has
> > no idea at all of the variability from a single sample.
>
> I thought that might be the reason. There is no mention of using n-1
> rather than n in the denominator on the help page. Perhaps that might be
> added?
Um. What professional statistics package uses n, then? Are you
suggesting that there is serious room for doubt?
--
Brian D. Ripley, ripley@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 272860 (secr)
Oxford OX1 3TG, UK Fax: +44 1865 272595
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel 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-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._