Posted by
Maguin, Eugene on
May 20, 2019; 8:57pm
URL: http://spssx-discussion.165.s1.nabble.com/negative-binomial-regression-tp5737913p5737915.html
I haven't seen the UCLA example that Bruce referenced. But isn't there something unexpected about this statement
(582 wards are columns; variables are rows)
Gene Maguin
-----Original Message-----
From: SPSSX(r) Discussion <
[hidden email]> On Behalf Of keltre19
Sent: Monday, May 20, 2019 11:50 AM
To:
[hidden email]
Subject: negative binomial regression
urgent help needed: I am using SPSS 25 to perform a negative binomial regression on overdispersed data. I am trying to work out the negative binomial equation/interpretation (582 wards are columns; variables are rows) = Dependent variable is Deaths 2011 (numeric figure for each 582 wards) + Independent variables are: categorical treated as factors and the rest are % proportions treated as covariates.
My Q: controlling for demographics and health variables, it was hypothesised than an increase in an areas deprivation and farming activities will statistically increase the risk of ward mortality.
1. How do I make predictions of deaths as I am at my wits end trying to
figure it out!! Going around and round in circles and getting nowhere :C
2. The code I am running is not doing what I want it to. I want to look at
deprivation (NIMDM) at 1 (most deprived, and then at other intervals for
comparisons) and see what the predicted deaths are in wards and its influence this has on the other IVs. But when I run this code below it makes
no difference to running it without the /EMMEANS
CONTROL=NIMDMrank2010(1) line
3. Then trying to place them on a scatterplot for comparisons. But when I
run this code which of the 4 variables created do I plot (if any). i know XPPRED gives me the predicted log counts but how do I make them meaningful in a figure.
GENLIN deaths11_fig BY farms11code avSOfarm_abovBEL grass_NEW crops_NEW catttle11code sheep11code
pigs11code poultry11code (ORDER=DESCENDING) WITH pcnt65to100A
SEXmalprcnt11 tot_livalone11 Qbeldeg
FT31plushrs NIMDMrank2010 yesUNPcar LLTIyesTOT OCCprcnt51 LFAfarms11 FARMEDpcnt_NEW
/MODEL farms11code avSOfarm_abovBEL grass_NEW crops_NEW catttle11code sheep11code pigs11code
poultry11code pcnt65to100A SEXmalprcnt11 tot_livalone11 Qbeldeg FT31plushrs NIMDMrank2010 yesUNPcar
LLTIyesTOT OCCprcnt51 LFAfarms11 FARMEDpcnt_NEW INTERCEPT=YES
OFFSET=NL_popul11
DISTRIBUTION=NEGBIN(MLE) LINK=LOG
/CRITERIA METHOD=FISHER(1) SCALE=1 COVB=ROBUST MAXITERATIONS=1000
MAXSTEPHALVING=50
PCONVERGE=1E-006(ABSOLUTE) SINGULAR=1E-012 ANALYSISTYPE=3(WALD)
CILEVEL=95 CITYPE=WALD
LIKELIHOOD=FULL
/EMMEANS CONTROL=NIMDMrank2010(1)
/MISSING CLASSMISSING=EXCLUDE
/PRINT CPS DESCRIPTIVES MODELINFO FIT SUMMARY SOLUTION (EXPONENTIATED)
/SAVE MEANPRED XBPRED COOK DEVIANCERESID.
I am hoping someone on here knows the syntax or point-and-click method to solve this issue as I have exhausted youtube videos and reading SPSS syntax books and journals, perhaps someone knows the solution off-hand, thanks in advance for any advice, feel free to ask any questions, Kel
--
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
=====================
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