Weekly Lab 3

library(readxl)

## Warning: package 'readxl' was built under R version 3.4.4

WeeklyLab3Data <- read_excel("~/Downloads/Weekly Lab datafiles-20180616/WeeklyLab3Data.xlsx")
View(WeeklyLab3Data)

library(readxl)
data <- read_excel("~/Downloads/Weekly Lab datafiles-20180616/WeeklyLab3Data.xlsx")

d <- density(data$Rating)
plot(d)

shapiro.test(data$Rating)

## 
##  Shapiro-Wilk normality test
## 
## data:  data$Rating
## W = 0.84697, p-value = 0.007589

library(moments)
agostino.test(data$Rating)

## 
##  D'Agostino skewness test
## 
## data:  data$Rating
## skew = -1.0248, z = -2.0462, p-value = 0.04074
## alternative hypothesis: data have a skewness

data$Rating2 <- data$Rating^2 
agostino.test(data$Rating2) 

## 
##  D'Agostino skewness test
## 
## data:  data$Rating2
## skew = -0.50677, z = -1.07270, p-value = 0.2834
## alternative hypothesis: data have a skewness

data$Ad <- factor(data$Ad)
data$Day <- factor(data$Day)
data$Audience <- factor(data$Audience)
data$Time <- factor(data$Time)

bartlett.test(data$Rating2~data$Ad)

## 
##  Bartlett test of homogeneity of variances
## 
## data:  data$Rating2 by data$Ad
## Bartlett's K-squared = 6.0618, df = 2, p-value = 0.04827

tapply(data$Rating2, data$Ad, var)

##         1         2         3 
##  678.5667 1322.9667  102.6667

model <- aov(Rating2~Day+Day/Audience+Time+Ad, data = data) 

library(xtable)
table <- xtable(model)

qqnorm(model$residuals)

TukeyHSD(model, "Ad")

##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = Rating2 ~ Day + Day/Audience + Time + Ad, data = data)
## 
## $Ad
##          diff        lwr       upr     p adj
## 2-1 -45.00000 -77.274460 -12.72554 0.0100098
## 3-1 -15.83333 -48.107793  16.44113 0.3846714
## 3-2  29.16667  -3.107793  61.44113 0.0749544

tapply(data$Rating2, data$Ad, mean)

##        1        2        3 
## 80.16667 35.16667 64.33333

tapply(data$ Rating2, data$Ad, sd)

##        1        2        3 
## 26.04931 36.37261 10.13246