Thanks for this.
I am no expert on gllam, but I ran side by side what you propose and
xtpoisson accident op_75_79 co_65_69 co_70_74 co_75_79, ///
fe i(ship) offset(ln_service) irr
I get coefficient estimates that are close, but not the same.
If I type
xtpoisson accident op_75_79 co_65_69 co_70_74 co_75_79, ///
re i(ship) offset(ln_service) irr
Then I get close to the results given by gllam. So I have one
question: can you get gllam to estimate a fixed effects specification
rather than a random effects one?
Regards,
Pierre
