The chi-square value is presented in the goodness-of-fit table, so getting this into a dataset with ODS OUTPUT should not be difficult. I do worry that the goodness of fit applies to a model with some terms in it, so you may have to change to: proc genmod; class field; freq observ_freq; model class.insect =field / dist=binomial; run; And of course that makes it impossible to find out if different distributions are better applied in different fields. So, try your method first, and if it gives no results in the goodness of fit table, you could try regressing on a dummy variable, such as regressing on a constant. Steve Denham
... View more