anova_note
Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| anova_note [2025/09/18 21:03] – [with 2 levels] hkimscil | anova_note [2025/09/18 21:05] (current) – [with more than 3 levels] hkimscil | ||
|---|---|---|---|
| Line 108: | Line 108: | ||
| </ | </ | ||
| ====== with more than 3 levels ====== | ====== with more than 3 levels ====== | ||
| - | |||
| < | < | ||
| # | # | ||
| Line 114: | Line 113: | ||
| # | # | ||
| rm(list=ls()) | rm(list=ls()) | ||
| - | # set.seed(101) | ||
| rnorm2 <- function(n, | rnorm2 <- function(n, | ||
| ss <- function(x) { | ss <- function(x) { | ||
| Line 120: | Line 118: | ||
| } | } | ||
| - | n <- 21 | + | set.seed(11) |
| + | n <- 31 | ||
| na <- nb <- nc <- nd <- n | na <- nb <- nc <- nd <- n | ||
| mean.a <- 98 | mean.a <- 98 | ||
| - | mean.b <- 100 | + | mean.b <- 99 |
| - | mean.c <- 101 | + | mean.c <- 102 |
| - | mean.d <- 105 | + | mean.d <- 103 |
| A <- rnorm2(na, mean.a, sqrt(900/ | A <- rnorm2(na, mean.a, sqrt(900/ | ||
| Line 157: | Line 156: | ||
| br <- seq(floor(min.x), | br <- seq(floor(min.x), | ||
| # Example bin width of 1 | # Example bin width of 1 | ||
| + | |||
| hist(A, breaks=br, | hist(A, breaks=br, | ||
| Line 263: | Line 263: | ||
| # graph 로 이해 | # graph 로 이해 | ||
| - | x <- rf(500000, df1 = df.bet, df2 = df.wit) | + | x <- rf(50000, df1 = df.bet, df2 = df.wit) |
| - | y.max <- max(df(x, | + | y.max <- max(df(x, df1=df.bet, df2=df.wit)) |
| hist(x, | hist(x, | ||
| Line 339: | Line 339: | ||
| summary(lm.res) | summary(lm.res) | ||
| summary(a.res) | summary(a.res) | ||
| - | |||
| </ | </ | ||
| + | |||
| + | {{: | ||
| + | {{: | ||
| + | {{: | ||
anova_note.1758196980.txt.gz · Last modified: 2025/09/18 21:03 by hkimscil