Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Proc nlmixed GEE for ZINB

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

☑ This topic is **solved**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 09-01-2023 10:04 AM
(460 views)

Hello,

I am looking to fit a GEE ZINB model using proc nlmixed. I have turned to proc nlmixed as neither proc genmod or proc GEE can handle zero-inflated distributions for GEE analysis.

Following SAS sources, such as note 2 here , I believe I can set up the nlmixed statement correctly as needed. However, I am not sure how to specify clustering. In general, my code looks like:

proc nlmixed data = data; /* linear predictor for zero inflation */ eta_zero = a0; prob_zero = 1/(1 + exp(-eta_zero)); /* linear predictor for negative binomial */; eta_mean = b0 + b1*dTLN + b2*dNLLN + b3*dT; mean = exp(eta_mean + myoff); IF y= 0 THEN lglk = LOG(prob_zero + (1-prob_zero)*(1+(phi*mean))**(1/phi)); ELSE lglk = LOG(1-prob_zero) + y*LOG(phi*mean) - (y+(1/phi))*LOG(1+(phi*mean)) + lgamma(y+(1/phi)) - lgamma(1/phi) - lgamma(y+1); model y ~ general(lglk);

random ??? ; run;

I understand proc nlmixed can utilize a 'random' statement which can include a subject identifier, but I do not want to include any random effects (I am trying to avoid a linear mixed effects model). I would just like to account for clustering/correlation as a subject can provide more than one response. Note, the data is not a longitudinal study setup.

Is there a way to manipulate the random statement to allow for some correlation structure but not to estimate any random effects, setting up for interpretation as a GEE? Or is there another statement/line in nlmixed I can do this?

Thanks!

- Tags:
- GEE
- Proc nlmixed
- zinb

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

A marginal, GEE type model uses R-side random effects and is not supported in NLMIXED which only allows for G-side effects. That is, the equivalent to the _RESIDUAL_ keyword in the RANDOM statement is not available in NLMIXED like it is in GLIMMIX. The best you can do would be to add random effects for your subjects as you suggest.

3 REPLIES 3

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

No, I'm just noting, syntactically and statistically, what NLMIXED lacks and would need to fit that model

**Available on demand!**

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.