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

CFC_SAS_Communities_400x225.jpg

Call for content now open!

It's your turn to help shape SAS Innovate 2027. Share your expertise and inspire the SAS community.

Submit your proposal →

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
  • 1475 views
  • 2 likes
  • 2 in conversation