Differences

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

--- c:ma:2019:lecturer [2019/10/01 23:52] – [t.test: mtcars] hkimscil
+++ c:ma:2019:lecturer [2019/10/04 11:28] (current) – [anova: mtcars] hkimscil
@@ Line 230: / Line 230: @@
    return(out)
 }
-stats <- stats4each(mtcars$hp, mtcars$cyl)
+library(MASS)
+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)
 stats
 sswithin <- sum(stats[5,])
-sstotal <- var(mtcars$hp)*(length(mtcars$hp)-1)
+sstotal <- var(y)*(length(y)-1)
 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
-dftotal <- length(mtcars$hp)-1
 dfwithin
@@ Line 253: / 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(mtcars$hp~mtcars$cyl, data=mtcars)
+mod <- aov(y~x, data=tempd)
 summary(mod)