Advertisement
Guest User

lefty

a guest
Mar 3rd, 2014
152
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 0.41 KB | None | 0 0
  1. data.frame(y = c(332928,358323,422244,430035,3168,2397,3099,5508,23802,26346,28546,33933,2286,2370,5667,9468,6024,5370,6126,7971,5979,4893,5007,6381,17985,18771,24081,26532),t=rep(c(1996,2001,2006,2013),7),type=rep(c("car","mb","bus","train","other","bike","walk"),each=4)) -> dat
  2.  
  3. dat$t96 <- dat$t-1996
  4. dat$type <- relevel(dat$type, "car")
  5.  
  6. model <- lm(log(y) ~ t96*type, data=dat)
  7. summary(model)
  8. anova(model)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement