powerDivStat {MSCquartets}

Power divergence statistic of Cressie & Read

Description

Computes any of the family of power-divergence statistics for multinomial data of Cressie and Read (1984), to compare observed and expected counts. Includes Likelihood Ratio and Chi-squared statistics as special cases.

Usage

powerDivStat(obs, expd, lambda)

Arguments

obs

observation vector

expd

expected vector

lambda

statistic parameter (e.g., 0=Likelihood Ratio, 1=Chi-squared)

Value

value of statistic

References

Cressie N, Read TRC (1984). “Multinomial Goodness-Of-Fit Tests.” J. Royal Stat. Soc. B, 46(3), 440-464.

Examples

obs=c(10,20,30)
expd=c(20,20,20)
powerDivStat(obs,expd,0)