where b is a vector of regression parameters to be estimated. Furthermore, it is
assumed that yi is Poisson distributed having expectation E(exp(x′
ib)). The Poisson
regression model can be fitted by maximizing its loglikelihood function, which is
often done by an iteratively reweighted least squares algorithm. Besides the model
parameters bkℓ, also standard errors σkl can be derived by this algorithm.