## ----setup,echo=FALSE,eval=TRUE---------------------------------------------------------------------
knitr::opts_knit$set(progress = TRUE, verbose = TRUE, concordance = TRUE)


## ----load_data--------------------------------------------------------------------------------------
options(contrasts=c("contr.sum","contr.poly") )
options(digits=4) # formats the output
load(url("http://pnb.mcmaster.ca/bennett/psy710/datasets/rose.Rdata") )


## ----r0,echo=TRUE,eval=TRUE-------------------------------------------------------------------------
rose$ageGroup <- factor(rose$age, labels = "g", ordered = FALSE)
y.means <- tapply(rose$score, rose$ageGroup, mean)
x.group <- sort(unique(rose$age))


## ----cor, echo=TRUE,eval=TRUE-----------------------------------------------------------------------
cor(x.group, y.means)


## ---------------------------------------------------------------------------------------------------
rose.aov.01 <- aov(score ~ ageGroup, data = rose)
summary(rose.aov.01)


## ---------------------------------------------------------------------------------------------------
rose.lm.01 <- lm(score ~ ageGroup, data = rose)
anova(rose.lm.01)


## ---------------------------------------------------------------------------------------------------
contr.poly(n = 5,scores=c(8,9,10,11,12))


## ---------------------------------------------------------------------------------------------------
rose$orderedAge <- factor(rose$age, labels = "a", ordered = TRUE)
class(rose$orderedAge)
class(rose$ageGroup)
unique(rose$orderedAge)
unique(rose$ageGroup)


## ---------------------------------------------------------------------------------------------------
rose.aov.02 <- aov(score ~ orderedAge, data = rose)
summary(rose.aov.02)


## ---------------------------------------------------------------------------------------------------
summary(rose.aov.02, split = list(orderedAge = list(1, 2, 3, 4)))


## ---------------------------------------------------------------------------------------------------
summary(rose.aov.02,
        split = list(orderedAge = list(linear = 1, quadratic = 2,cubic = 3, quartic = 4)))


## ---------------------------------------------------------------------------------------------------
summary(rose.aov.02, split = list(orderedAge = list(linear = 1, nonlinear = 2:4)))


## ---------------------------------------------------------------------------------------------------
summary(rose.aov.02, split = list(orderedAge = list(lin = 1, quad=2, other = 3:4)))


## ---------------------------------------------------------------------------------------------------
coef(rose.aov.02)
confint(rose.aov.02) # confidence intervals


## ---------------------------------------------------------------------------------------------------
with(rose, tapply(score,orderedAge,mean) ) # calc group means