> sum((x-mean(x))^2)/(n) [1] 0.4894708 > ((n-1)/n) * var(x) [1] 0.4894708
hth, Kingsford On Mon, Oct 19, 2009 at 9:30 AM, Peng Yu <pengyu...@gmail.com> wrote: > It seems that var() computes sample variance. It is straight forward > to compute population variance from sample variance. However, I feel > that it is still convenient to have a function that can compute > population variance. Is there a population variance function available > in R? > > $ Rscript var.R >> set.seed(0) >> n = 4 >> x = rnorm(n) >> var(x) > [1] 0.6526278 >> sum((x-mean(x))^2)/(n-1) > [1] 0.6526278 >> > > ______________________________________________ > R-help@r-project.org mailing list > 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. > ______________________________________________ R-help@r-project.org mailing list 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.