BookmarkSubscribeRSS Feed
Melk
Lapis Lazuli | Level 10

I am running a count model on number of treatments given and comparing them between 2 diseases. I have excess zeros, but there is no extra "Excess" zeros so I am using a negative binomial model and not a zero inflated model. One of the disease groups is small at about 20 patients. when I run an estimate statement it says "non-est" but I do get MLE parameter estimates. Can I just exponentiate the estimates to get the rate per group although the estimate statement that is doing this is not running?

 

proc genmod data = dat;
  class group;
  model numtx = group/ type3 link=log dist=negbin offset=loghospdays;
  estimate 'Asthma' group 1 0 / exp;
run;

1 REPLY 1
PGStats
Opal | Level 21

To get the group mean counts use

lsmeans group / ilink;

To get the ratio of the group mean counts use

 

lsmeans group / diff exp;
PG

sas-innovate-white.png

Register Today!

Join us for SAS Innovate 2025, our biggest and most exciting global event of the year, in Orlando, FL, from May 6-9. Sign up by March 14 for just $795.

Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

SAS Training: Just a Click Away

 Ready to level-up your skills? Choose your own adventure.

Browse our catalog!

Discussion stats
  • 1 reply
  • 709 views
  • 2 likes
  • 2 in conversation