Differences

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

--- c:ps1-1:ci_note [2024/11/19 23:06] – hkimscil
+++ c:ps1-1:ci_note [2024/11/19 23:18] (current) – hkimscil
@@ Line 5: / Line 5: @@
 p1 <- rnorm2 (100000000, mean=70, sd=10)
-a1 <- qnorm(.95)
-a2 <- qnorm(.975)
-a3 <- qnorm(.995)
-a1
-a2
-a3
+set.seed(101)
 s.size <- 16
 s1 <- sample(p1, s.size, replace=T)
@@ Line 21: / Line 16: @@
 se <- sd.s1/sqrt(s.size)
-se2 <- se*2
+se2 <- 2*se
-se3 <- se*3
+se3 <- 2*se
 m.s1 - se2
@@ Line 28: / Line 24: @@
 # or
-se2 <- se*a2
+c2 <- qnorm(.975)
-se3 <- se*a3
+c3 <- qnorm(.995)
+c2
+c3
+se2 <- c2*se
+se3 <- c3*se
 m.s1 - se2
 m.s1 + se2
@@ Line 35: / Line 37: @@
 <code>
-> a1 <- qnorm(.95)
+> rnorm2 <- function(n,mean,sd) {
-> a2 <- qnorm(.975)
+>    mean + sd * scale(rnorm(n))
-> a3 <- qnorm(.995)
+>    }
-> a1
+>
-[1] 1.644854
+> p1 <- rnorm2 (100000000, mean=70, sd=10)
-> a2
+>
-[1] 1.959964
+> set.seed(101)
-> a3
-[1] 2.575829
->
 > s.size <- 16
 > s1 <- sample(p1, s.size, replace=T)
 > s1
- [1] 67.33346 80.38354 71.51355 69.21986 60.04729 65.60540 69.31297 52.70064 76.02403
+ [1] 64.87746 78.68130 50.61449 62.18639 68.08899 63.07744 54.48954 69.26928 55.92358
-[10] 79.09528 60.57269 62.36968 69.86095 74.76315 64.24179 72.98447
+[10] 70.06123 60.38630 87.91813 78.55023 69.78835 63.15719 64.35994
 > m.s1 <- mean(s1)
 > sd.s1 <- sd(s1)
 > m.s1
-[1] 68.5018
+[1] 66.33936
 > sd.s1
-[1] 7.442683
+[1] 9.59219
 >
 > se <- sd.s1/sqrt(s.size)
-> se2 <- se*2
+> se2 <- 2*se
-> se3 <- se*3
+> se3 <- 2*se
+>
 >
 > m.s1 - se2
-[1] 64.78046
+[1] 61.54327
 > m.s1 + se2
-[1] 72.22314
+[1] 71.13546
 >
 > # or
-> se2 <- se*a2
+> c2 <- qnorm(.975)
-> se3 <- se*a3
+> c3 <- qnorm(.995)
+> c2
+[1] 1.959964
+> c3
+[1] 2.575829
+>
+> se2 <- c2*se
+> se3 <- c3*se
+>
 > m.s1 - se2
-[1] 64.85495
+[1] 61.63928
 > m.s1 + se2
-[1] 72.14864
+[1] 71.03945
->
 >
 </code>