load("bengal.dat")
K <- 20; t0 <- 101; Delta <- t0/K
breaks <- c(0:K)*Delta
(y <- hist(bengal-1877,breaks=breaks,plot=FALSE)$counts)
t <- Delta*(c(0:(K-1))+0.5)
log.Delta <- rep(1,K)*log(Delta)
summary(glm(y~1+t+offset(log.Delta),family=poisson))


c <- cos(2*pi*t)
s <- sin(2*pi*t)
summary(glm(y~t+offset(log.Delta)+s+c,family=poisson))