setwd("C:/Users/sharl/Desktop/USF/Fall 2021/LIS 4273 - Adv Stats & Analytics/Modules/Week 8 - ANOVA")
HighStress <- c(10, 9, 8, 9, 10, 8)
ModerateStress <- c(8, 10, 6, 7, 8, 8)
LowStress <- c(4, 6, 6, 4, 2, 2)
stress <- data.frame(cbind(HighStress, ModerateStress, LowStress))
stressStack <- stack(stress)
oneway.test(values ~ ind, data = stressStack)
##
## One-way analysis of means (not assuming equal variances)
##
## data: values and ind
## F = 17.525, num df = 2.0000, denom df = 9.2777, p-value = 0.0007065
aovResult <- aov(values ~ ind, data = stressStack)
summary(aovResult)
## Df Sum Sq Mean Sq F value Pr(>F)
## ind 2 82.11 41.06 21.36 4.08e-05 ***
## Residuals 15 28.83 1.92
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
library(ISwR)
data("zelazo")
attach(zelazo)
str(zelazo)
## List of 4
## $ active : num [1:6] 9 9.5 9.75 10 13 9.5
## $ passive: num [1:6] 11 10 10 11.8 10.5 ...
## $ none : num [1:6] 11.5 12 9 11.5 13.2 ...
## $ ctr.8w : num [1:5] 13.2 11.5 12 13.5 11.5
zelazo_df = data.frame(active, passive, none)
zStack = stack(zelazo_df)
oneway.test(values ~ ind, data = zStack)
##
## One-way analysis of means (not assuming equal variances)
##
## data: values and ind
## F = 1.7479, num df = 2.0000, denom df = 9.8852, p-value = 0.2239
aovResult2 = aov(values ~ ind, data = zStack)
summary(aovResult2)
## Df Sum Sq Mean Sq F value Pr(>F)
## ind 2 8.36 4.181 1.568 0.241
## Residuals 15 39.99 2.666