I am wondering if guidance/suggestions can be offered about the comparability of results for negative binomial regression from GENMOD and COUNTREG in SAS 9.3. I am using a Statistics Canada public use microdata file with non-integer weights. I know that GENMOD uses GLM-style coding for classification variables but is this also true for COUNTREG? (It appears to be so but I have seen no documentation to that effect though an example strongly hints that this is so.) Resulting parameter estimates look different (signs are same but magnitudes different in the 1 st or 2 nd decimal place). I am using defaults in both procedures (Newton-Raphson in COUNTREG and I assume a maximum likelihood in GENMOD); both procedures converge. Though I am quite content to accommodate such differences, is there a general preference of one procedure over the other based on the default algorithms and types of data? More interestingly, COUNTEG has a ‘probcount’ keyword in its output statement of the "probability of the response variable taking particular values." But no such facility appears in GENMOD. GENMOD can provide the mean of "the predicted probability that the response variable is less than or equal to the value of _LEVEL_ if the multinomial model for ordinal data is used (in other words, Pr , where Y is the response variable)," which is manifestly not the same as the output from COUNTREG. If I wanted to get the equivalent of PROBCOUNT from GENMOD is there a macro that would enable me to do so?
... View more