Login  Register

Re: negative binomial regression

Posted by Bruce Weaver on May 21, 2019; 9:14pm
URL: http://spssx-discussion.165.s1.nabble.com/negative-binomial-regression-tp5737913p5737932.html

I still don't completely understand what you are trying to do.  But I suspect
that your EMMEANS subcommands need more to them than you show below to get
what you want.  Here is a GENLIN example from one of my syntax files.  I was
using GENLIN to estimate a binary logit model, not a negative binomial
model.  But perhaps the EMMEANS sub-commands will give you an idea about
what you need to do.

GENLIN collegegrad (REFERENCE=FIRST)
       BY WHITE male (order = descending)
       WITH age speduc
  /MODEL
  WHITE
  speduc speduc*speduc
  age age*age
  male
  WHITE*speduc
  speduc*age
  speduc*age*age
   INTERCEPT=YES
   DISTRIBUTION=BINOMIAL LINK=LOGIT
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(10) age(40) SCALE=TRANSFORMED
COMPARE=WHITE
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(10) age(50) SCALE=TRANSFORMED
COMPARE=WHITE
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(10) age(60) SCALE=TRANSFORMED
COMPARE=WHITE
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(15) age(40) SCALE=TRANSFORMED
COMPARE=WHITE
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(15) age(50) SCALE=TRANSFORMED
COMPARE=WHITE
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(15) age(60) SCALE=TRANSFORMED
COMPARE=WHITE
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(20) age(40) SCALE=TRANSFORMED
COMPARE=WHITE
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(20) age(50) SCALE=TRANSFORMED
COMPARE=WHITE
  /EMMEANS TABLES=male*WHITE CONTROL=speduc(20) age(60) SCALE=TRANSFORMED
COMPARE=WHITE
  /MISSING CLASSMISSING=EXCLUDE
  /PRINT CPS DESCRIPTIVES MODELINFO FIT SUMMARY SOLUTION (EXPONENTIATED).


I used the SCALE=TRANSFORMED option to get log-odds, by the way.  

I didn't show it here, but I also used OMS to write the EMMEANS results to
another dataset, thus making them available for graphing.  

HTH.


keltre19 wrote
> --- snip ---
>
> The rest of the code seems to work for me, however it is when I include
> the
> line  
> /EMMEANS CONTROL=NIMDMrank2010(1)  
>  it doesn’t seem to do what I need it to,
>
> --- snip ---





-----
--
Bruce Weaver
[hidden email]
http://sites.google.com/a/lakeheadu.ca/bweaver/

"When all else fails, RTFM."

NOTE: My Hotmail account is not monitored regularly.
To send me an e-mail, please use the address shown above.

--
Sent from: http://spssx-discussion.1045642.n5.nabble.com/

=====================
To manage your subscription to SPSSX-L, send a message to
[hidden email] (not to SPSSX-L), with no body text except the
command. To leave the list, send the command
SIGNOFF SPSSX-L
For a list of commands to manage subscriptions, send the command
INFO REFCARD
--
Bruce Weaver
bweaver@lakeheadu.ca
http://sites.google.com/a/lakeheadu.ca/bweaver/

"When all else fails, RTFM."

PLEASE NOTE THE FOLLOWING: 
1. My Hotmail account is not monitored regularly. To send me an e-mail, please use the address shown above.
2. The SPSSX Discussion forum on Nabble is no longer linked to the SPSSX-L listserv administered by UGA (https://listserv.uga.edu/).