Hi, I am building a model as follows: proc genmod data = aa; model totalmedcost= PDC age sex race education /dist=gamma link=log; run; Here, I need to check if the estimate of the major predictor PDC is appropriate or unbiased. One of methods is to plot the error of the model (deviance) against PDC to see whether the error is randomly scattered around PDC; another way is to examine the association between deviance and PDC. However, when I use PLOTS=RESDEV(XBETA), it only generates the plots of deviance residuals as a function of linear predictor, not necessarily the major predictor. Q1: How could I get the plot of deviance against the specific predictor I want after controlling all the covariates? Q2: How to get the deviance residual (DR) for each observation, so that I could run the regression of (DR) agaist the predictor PDC? Thanks for help! ... View more

Hi, I am building a multivariable model using proc genmod (gamma and a log link), and I have a major explanatory variable (X1, which is 5-level category variable) and a few covariates. When I got the output, I had estimates for intercept B0, B1(X1), and estimates for other covairables....but I could not find the overall mean value for the dependent variable in the whole model. Also, I want to calculate the group means by variable X1 in the model. Are there any ways to calculate the overall mean of dependent variable and the group means by X1 in the gamma and a log link model? Thanks a lot! ... View more