# An analysis of the data from Table 2.7 similar to the SAS analysis

a <- c(2.2, 1.6, 0.8, 1.8, 1.4, 0.4, 0.6, 1.5, 0.5)

v <- c(0.3, 0.0, 0.6, 0.0, -0.3, 0.2)

s <- c(0.1, 0.1, 0.2, -0.4, 0.3, 0.1, 0.1, -0.5)

treatment <- c(rep('amiodarone',length(a)),rep('vehicle',length(v)),rep('saline',length(s)))

rabbitears <- data.frame(treatment,c(a,v,s))

names(rabbitears) <- c("treatment","temperaturediff")

rabbitears$treatment <- as.factor(rabbitears$treatment)    

boxplot(temperaturediff ~ treatment, data=rabbitears ,main="Rabbit data")

anova(lm(temperaturediff ~ treatment, data=rabbitears))

windows()
par(mfrow=c(2,2))  
plot(lm(temperaturediff ~ treatment, data=rabbitears))
par(mfrow=c(1,1))