Nothing
# these look nicer if we parameterize differently in the model
airp.lm2 <- lm(pollution~-1 + location,airpollution)
contr2 <- rbind(
"hill - plains" = c(1,-1,0),
"suburb - urban" = c(1,1,-2))
###hop:3-6
summary(glht(airp.lm2,contr2))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.