Hello everyone,

I am running an analysis to see if serum cholesterol together with sex, DBP, and Age are associated with number of heart attacks. I started off by checking for the assumption of (mean=variance) and noticed that its better to work with negative binomial regression due to overdispersion.

I went on to fit the model and zero inflated negative binomial regression turned out to be the best model I could work with.

Compared the -2LL and noticed the full model was better than the reduced model. However am stuck on how to proceed to the next step of dropping all insignificant interaction terms, which p-values to consider in the model and how i can fit my model.

Kindly help me out.

Thanks