c:ma:2019:lecturer
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
c:ma:2019:lecturer [2019/10/03 06:58] – [anova: mtcars] hkimscil | c:ma:2019:lecturer [2019/10/04 11:28] (current) – [anova: mtcars] hkimscil | ||
---|---|---|---|
Line 231: | Line 231: | ||
} | } | ||
- | data <- mtcars | + | library(MASS) |
- | x <- mtcars$gear | + | |
- | y <- mtcars$mpg | + | tempd <- iris |
+ | x <- tempd$Species | ||
+ | y <- tempd$Sepal.Width | ||
+ | |||
+ | tempd <- mtcars | ||
+ | x <- tempd$gear | ||
+ | y <- tempd$mpg | ||
+ | |||
+ | tempd <- mtcars | ||
+ | x <- tempd$am | ||
+ | y <- tempd$mpg | ||
+ | |||
+ | |||
+ | x <- factor(x) | ||
+ | dfbetween <- nlevels(x)-1 | ||
stats <- stats4each(y, | stats <- stats4each(y, | ||
Line 242: | Line 256: | ||
ssbetween <- sstotal-sswithin | ssbetween <- sstotal-sswithin | ||
- | round(sswithin, | + | round(sswithin, |
- | round(ssbetween, | + | round(ssbetween, |
- | round(sstotal, | + | round(sstotal, |
dfwithin <- sum(stats[4, | dfwithin <- sum(stats[4, | ||
- | dfbetween <- 3-1 | ||
dftotal <- length(y)-1 | dftotal <- length(y)-1 | ||
Line 258: | Line 271: | ||
mstotal <- sstotal / dftotal | mstotal <- sstotal / dftotal | ||
- | round(mswithin, | + | round(mswithin, |
- | round(msbetween, | + | round(msbetween, |
- | round(mstotal, | + | round(mstotal, |
fval <- round(msbetween/ | fval <- round(msbetween/ | ||
fval | fval | ||
- | siglevel <- pf(q=36.18, df1=2, df2=29, lower.tail=FALSE) | + | siglevel <- pf(q=fval, df1=dfbetween, df2=dfwithin, lower.tail=FALSE) |
siglevel | siglevel | ||
- | mod <- aov(mtcars$hp~mtcars$cyl, data=mtcars) | + | mod <- aov(y~x, data=tempd) |
summary(mod) | summary(mod) | ||
c/ma/2019/lecturer.1570053495.txt.gz · Last modified: 2019/10/03 06:58 by hkimscil