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)
> head(a)
[,1]
[1,] 93.73555
[2,] 101.83562
[3,] 91.64403
[4,] 115.95044
[5,] 103.29411
[6,] 91.79562
> tail(a,20)
[,1]
[99999981,] 94.13305
[99999982,] 99.05620
[99999983,] 92.38097
[99999984,] 101.97245
[99999985,] 90.67295
[99999986,] 95.95812
[99999987,] 75.85622
[99999988,] 89.77369
[99999989,] 101.50076
[99999990,] 119.73697
[99999991,] 90.86393
[99999992,] 92.32930
[99999993,] 102.45092
[99999994,] 124.03421
[99999995,] 86.92416
[99999996,] 100.61147
[99999997,] 93.50381
[99999998,] 101.24958
[99999999,] 95.95971
[100000000,] 96.43605
> m.a <- mean(a)
> sse.a <- sum((a-m.a)^2)
> n.a <- length(a)
> df.a <- n.a - 1
>
> sse.a / n.a
[1] 100
> sse.a/ df.a
[1] 100
>
variance/ro2.1773184821.txt.gz · Last modified: by hkimscil