http://www.commres.net/ 2026-05-01T06:51:45+00:00 http://www.commres.net/ http://www.commres.net/_media/wiki/logo.png text/html 2025-11-30T22:51:24+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/summary_of_hypothesis_testing/code01?rev=1764543084&do=diff rm(list=ls()) rnorm2 <- function(n,mean,sd){ mean+sd*scale(rnorm(n)) } n.p <- 100000 m.p <- 100 sd.p <- 10 p1 <- rnorm2(n.p, m.p, sd.p) m.p1 <- mean(p1) sd.p1 <- sd(p1) p2 <- rnorm2(n.p, m.p+4, sd.p) m.p2 <- mean(p2) sd.p2 <- sd(p2) n.s <- 25 se.z1 <- c(sqrt(var(p1)/n.s)) se.z2 <- c(sqrt(var(p2)/n.s)) se.z1 se.z2 x.p1 <- seq(mean(p1)-5*se.z1, mean(p2)+5*se.z1, length.out = 500) x.p2 <- seq(mean(p2)-5*se.z1, mean(p2)+5*se.z1, length.out… text/html 2025-11-30T22:56:21+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/summary_of_hypothesis_testing/code03?rev=1764543381&do=diff 위의 code1 script와 동일 내용 text/html 2025-11-30T15:51:08+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/summary_of_hypothesis_testing/code2?rev=1764517868&do=diff rm(list=ls()) rnorm2 <- function(n,mean,sd){ mean+sd*scale(rnorm(n)) } n.p <- 100000 m.p <- 100 sd.p <- 10 p1 <- rnorm2(n.p, m.p, sd.p) m.p1 <- mean(p1) sd.p1 <- sd(p1) p2 <- rnorm2(n.p, m.p+4, sd.p) m.p2 <- mean(p2) sd.p2 <- sd(p2) n.s <- 900 se.z1 <- c(sqrt(var(p1)/n.s)) se.z2 <- c(sqrt(var(p2)/n.s)) se.z1 se.z2 x.p1 <- seq(mean(p1)-5*se.z1, mean(p2)+5*se.z1, length.out = 500) x.p2 <- seq(mean(p2)-5*se.z1, mean(p2)+5*se.z1, length.ou… text/html 2025-11-30T22:50:30+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/summary_of_hypothesis_testing/output01?rev=1764543030&do=diff > rm(list=ls()) > > rnorm2 <- function(n,mean,sd){ + mean+sd*scale(rnorm(n)) + } > > n.p <- 100000 > m.p <- 100 > sd.p <- 10 > p1 <- rnorm2(n.p, m.p, sd.p) > m.p1 <- mean(p1) > sd.p1 <- sd(p1) > > p2 <- rnorm2(n.p, m.p+4, sd.p) > m.p2 <- mean(p2) > sd.p2 <- sd(p2) > > n.s <- 25 > se.z1 <- c(sqrt(var(p1)/n.s)) > se.z2 <- c(sqrt(var(p2)/n.s)) > se.z1 [1] 2 > se.z2 [1] 2 > > x.p1 <- seq(mean(p1)-5*se.z1, + mean(p2)+5*se.z1, + length.out = 500) > x.p2 <- seq(mean(p… text/html 2025-11-30T22:55:29+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/summary_of_hypothesis_testing/output03?rev=1764543329&do=diff > rm(list=ls()) > > rnorm2 <- function(n,mean,sd){ + mean+sd*scale(rnorm(n)) + } > > n.p <- 100000 > m.p <- 100 > sd.p <- 10 > p1 <- rnorm2(n.p, m.p, sd.p) > m.p1 <- mean(p1) > sd.p1 <- sd(p1) > > p2 <- rnorm2(n.p, m.p+4, sd.p) > m.p2 <- mean(p2) > sd.p2 <- sd(p2) > > n.s <- 25 > se.z1 <- c(sqrt(var(p1)/n.s)) > se.z2 <- c(sqrt(var(p2)/n.s)) > se.z1 [1] 2 > se.z2 [1] 2 > > x.p1 <- seq(mean(p1)-5*se.z1, + mean(p2)+5*se.z1, + length.out = 500) > x.p2 <- seq(mean(p… text/html 2025-11-30T15:50:14+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/summary_of_hypothesis_testing/output2?rev=1764517814&do=diff > > rm(list=ls()) > > rnorm2 <- function(n,mean,sd){ + mean+sd*scale(rnorm(n)) + } > > n.p <- 100000 > m.p <- 100 > sd.p <- 10 > p1 <- rnorm2(n.p, m.p, sd.p) > m.p1 <- mean(p1) > sd.p1 <- sd(p1) > > p2 <- rnorm2(n.p, m.p+4, sd.p) > m.p2 <- mean(p2) > sd.p2 <- sd(p2) > > n.s <- 900 > se.z1 <- c(sqrt(var(p1)/n.s)) > se.z2 <- c(sqrt(var(p2)/n.s)) > se.z1 [1] 0.3333333 > se.z2 [1] 0.3333333 > > x.p1 <- seq(mean(p1)-5*se.z1, + mean(p2)+5*se.z1, + length.out = 500) …