variance:ro2
This is an old revision of the document!
> rnorm2 <- function(n,mean,sd) { mean+sd*scale(rnorm(n)) }
> set.seed(1)
> a <- rnorm2(100000000, 100, 10)
> a.mean <- mean(a)
> ss <- sum((a-a.mean)^2)
> n <- length(a)
> df <- n-1
> ss/n
[1] 100
> ss/df
[1] 100
variance/ro2.1773184538.txt.gz · Last modified: by hkimscil