MSCquartets: powerDivStat – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

powerDivStat

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)

MSCquartets

Analyzing Gene Tree Quartets under the Multi-Species Coalescent

v1.1.0

MIT + file LICENSE

Authors

Elizabeth Allman [aut], Hector Banos [aut], Jonathan Mitchell [aut], John Rhodes [aut, cre] (<https://orcid.org/0000-0001-9921-1091>)

Initial release