- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Posted 10-18-2018 03:50 PM
(749 views)
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
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
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