The Pearson Score Statistic for Multinomial-Poisson Models

The score statistic S² is commonly used for general likelihood-based inference. Pearson’s Chi-squared statistic X² = ∑(O ? E)²/E is ubiquitous in contingency table inference. Because tests and confidence intervals based on S² have been shown to work well in practice and theory and because X² has such a simple and intuitively appealing form, it is of interest to know when S² is identical to X² and when X² has an approximate Chi-squared distribution. Toward these ends, this paper gives a simple proof that S² = X² for the broad class of multinomial-Poisson distributions when the alternative hypothesis is unrestricted in a certain sense. This paper also gives a sufficient condition under which the null distribution of the Pearson score statistic is approximately Chi-squared. Several examples illustrate the utility of the results and counter-examples highlight the importance of the sufficient conditions of the results.