Some NA's. Maybe due to glm.my()
Some error on adding offset into both fit and predictions