Differences

This shows you the differences between two versions of the page.

--- variance:ro2 [2026/03/10 23:20] – hkimscil
+++ variance:ro2 [2026/03/10 23:29] (current) – removed hkimscil
@@ Line 1: / Line 1: @@
-<code>
-> 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
->
-</code>