Differences

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

--- c:ma:2019:lecturer [2019/10/03 06:59] – [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, x)
@@ Line 242: / Line 256: @@
 ssbetween <- sstotal-sswithin
-round(sswithin,0)
+round(sswithin,2)
-round(ssbetween,0)
+round(ssbetween,2)
-round(sstotal,0)
+round(sstotal,2)
 dfwithin <- sum(stats[4,])
-dfbetween <- 3-1
 dftotal <- length(y)-1
@@ Line 258: / Line 271: @@
 mstotal <- sstotal / dftotal
-round(mswithin,0)
+round(mswithin,2)
-round(msbetween,0)
+round(msbetween,2)
-round(mstotal,0)
+round(mstotal,2)
 fval <- round(msbetween/mswithin,2)
 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
-mod <- aov(y~x, data=data)
+mod <- aov(y~x, data=tempd)
 summary(mod)