p-Values of Multinomial Goodness-of-Fit Tests by Fourier Series Inversion R package MultinomFourier version 1.0 Computes p-values of multinomial goodness-of-fit tests based on an additive statistic: the log-likelihood ratio, Pearson's chi-squared, the Cressie-Read power divergence, and the exact multinomial test. The main function, pval_fourier(), writes the p-value as a Fourier series in the moment-generating function of the statistic and sums it until an adaptive stopping rule is met, which makes small p-values of large instances computationally feasible. The function pval_exact() computes the exact p-value by branch pruning with bounds from a dynamic program, and pval_flexible() chooses between the two according to the size of the instance. Software Imports: parallel Suggests: knitr, rmarkdown, ggplot2 Pablo Rivas <pablo.rivas.g@ug.uchile.cl> Comprehensive R Archive Network (CRAN) Charles Thraves [aut], Pablo Rivas [aut, cre], Felipe Subiabre [aut] GPL (>= 2) 2026-10-29 application/tgz https://CRAN.R-project.org/package=MultinomFourier doi:10.32614/CRAN.package.MultinomFourier