I have a negative binomial model where the outcome is the number of complications and my main predictors are three different continuous factor scores that relate to where a disease is found in the body. I found an interaction between one of the factor scores (factor3) and how long, in year, there was a delay from first experiencing symptoms and subsequently receiving surgery. I tried to make an estimate statement but all the mean estimates come back as zero and there's no exp findings. I am looking at this using the mean factor score of each quartile of the variable (so the mean of quartile 1 of factor3 is -1.18, mean of q2 is -0.33, etc) and for diagnostic delay I want to know the difference between people with no diagnostic delay because their disease was found incidentially (diagnostic delay of 0) and those who have an average diagnostic delay in my population (9.8 years). I am using sas 9.4. Any help is greatly appreciated. My simplified code is below proc genmod data =data;
model num_complications = factor1 factor2 factor3 white BMI
diagnosticdelay alcohol_use smoking
factor3*diagnosticdelay
/type3 dist=negbin;
estimate 'average q1 by incidential findings' intercept 1 diagnosticdelay 0 factor3 -1.18 diagnosticdelay*factor3 0/ exp;
estimate 'average q2 by incidential findings' intercept 1 diagnosticdelay 0 factor3 -0.33 diagnosticdelay*factor3 0/ exp;
estimate 'average q3 by incidential findings' intercept 1 diagnosticdelay 0 factor3 0.17 diagnosticdelay*factor3 0/ exp;
estimate 'average q4 by incidential findings' intercept 1 diagnosticdelay 0 factor3 1.30 diagnosticdelay*factor3 0/ exp;
estimate 'average q1 by average delay' intercept 1 diagnosticdelay 9.8 factor3 -1.18 diagnosticdelay*factor3 -113.3272/ exp;
estimate 'average q2 by average delay' intercept 1 diagnosticdelay 9.8 factor3 -0.33 diagnosticdelay*factor3 -3.234/ exp;
estimate 'average q3 by average delay' intercept 1 diagnosticdelay 9.8 factor3 0.17 diagnosticdelay*factor3 1.66/ exp;
estimate 'average q4 by average delay' intercept 1 diagnosticdelay 9.8 factor3 1.30 diagnosticdelay*factor3 12.74/ exp;
run;
... View more