# Simulating fake data to illustrate the fit of a Poisson regression
n = 100
x = rnorm(n=n,mean=3)
a = .6 ; b = -1
l = a*x+b
mu = exp(l)
y = rpois(n=n,lambda = mu)
# End Simulating fake data to illustrate fiiting of a Poisson regression


plot(x,y)

# Fitting of Poisson regression
dat = data.frame(nb.cycles=y,dist=x)

glm.res = glm(formula = (nb.cycles ~ dist),
              family = poisson(link="log"),
              data = dat)

summary(glm.res)

plot(glm.res)