# Poissonin regressio
# Ks. Altham s. 38
#


##########################

Onnettomuudet Cambridgessa 1978-81
Altham, s. 39
##########################

# y = onnettomuuksien lkm
y <- scan()
11 9 4 4 20 4
 
# v = arvioitu liikennemäärä
v <- scan()
2206 3276 1999 1399 2276 1417
 
rd <- c(1,1,1,2,2,2) 
# rd <- c("TR","TR","TR","MR","MR","MR")
# rd = 1 on Trumpington rd, rd = 2 on Mill rd
ToD <- c(1,2,3,1,2,3); ToD <- rep(c(1,2,3),2)
# ToD = 1,2,3 tarkoittaa aikoja 7-9.30, 9.30-15, 15-18.30.
# Huom. että muuttujan nimen "t" käyttö saattaa johtaa viheilmoitukseen,
# koska t() on R:n funktio (transpoosi)

lv <- log(v); RD <- factor(rd); ToD <- factor(ToD)