I was trying to fit neg binomial and poisson distribution parameters to a single data set (no predictors just a set of values) using Maximun Likelihood Estimation and then conduct a likelihood ratio test to see if the neg binomial fits the data better than the poisson. The problem is that I only have SAS Base 9.1 and SAS Stat, I don’t have SAS IML. The only option that I have was working with SAS Genmode procedure, but this is only for regression models. However, I found two things running the procedure with my data set as response variable and with no predictors. First, the dispersion parameter displayed in the “Analysis Of Maximum Likelihood Parameter Estimates” table for the negative binomial, gives the same number, for the distribution parameter r, I get using other methods (in R) that are meant to work with a single data set (no regression analysis). Second, I calculate the test statistic of the likelihood ratio test, using the loglikelihood values displayed in the genmode output and I get the same value for the test statistic too. My question is why this happens and can I make my distributional analysis using the genmode procedure in this way (with no predictors)?
