MultinomCI {DescTools}  R Documentation 
Confidence intervals for multinomial proportions are often approximated by single binomial confidence intervals, which might in practice often yield satisfying results, but is properly speaking not correct. This function calculates simultaneous confidence intervals for multinomial proportions either according to the methods of Sison and Glaz, Goodman, Wald, Wald with continuity correction or Wilson.
MultinomCI(x, conf.level = 0.95, sides = c("two.sided", "left", "right"),
method = c("sisonglaz", "cplus1", "goodman", "wald", "waldcc",
"wilson", "qh", "fs"))
x 
A vector of positive integers representing the number of occurrences of each class. The total number of samples equals the sum of such elements. 
conf.level 
confidence level, defaults to 0.95. 
sides 
a character string specifying the side of the confidence interval, must be one of 
method 
character string specifing which method to use; can be one out of

Given a vector of observations with the number of samples falling in each class of a multinomial distribution,
builds the simultaneous confidence intervals for the multinomial probabilities according to the method proposed by the mentioned authors.
The R code for Sison and Glaz (1995) has been translated from thes SAS code written by May and Johnson (2000). See the references for the other methods (qh = QuesensberryHurst, fs = FitzpatrickScott).
Some approaches for the confidence intervals can potentially yield negative results or values beyond 1. These would be reset such as not to exceed the range of [0, 1].
A matrix with 3 columns:
est 
estimate 
lwr.ci 
lower bound of the confidence interval 
upr.ci 
upper bound of the confidence interval 
The number of rows correspond to the dimension of x.
Pablo J. Villacorta Iglesias <pjvi@decsai.ugr.es>
Department of Computer Science and Artificial Intelligence, University of Granada (Spain)
(SisonGlaz)
Andri Signorell <andri@signorell.net> (Goodman, Wald, Wilson, FitzpatrickScott, QuesensberryHurst)
Fitzpatrick, S. and Scott, A. (1987). Quick simultaneous confidence interval for multinomial proportions. Journal of American Statistical Association 82(399): 875878.
Glaz, J., Sison, C.P. (1999) Simultaneous confidence intervals for multinomial proportions. Journal of Statistical Planning and Inference 82:251262.
Goodman, L. A. (1965) On Simultaneous Confidence Intervals for Multinomial Proportions Technometrics, 7, 247254.
May, W.L., Johnson, W.D.(2000) Constructing twosided simultaneous confidence intervals for multinomial proportions for small counts in a large number of cells. Journal of Statistical Software 5(6) . Paper and code available at https://www.jstatsoft.org/v05/i06.
Quesensberry, C.P. and Hurst, D.C. (1964). Large Sample Simultaneous Confidence Intervals for Multinational Proportions. Technometrics, 6: 191195.
Sangeetha, U., Subbiah, M., Srinivasan, M. R. (2013) Mathematical Analysis of propensity of aberration on the methods for interval estimation of the multinomial proportions. IOSR Journal of Mathematics, eISSN: 22785728,pISSN: 2319765X, Volume 7, Issue 4 (Jul.  Aug. 2013), PP 2328
Sison, C.P and Glaz, J. (1995) Simultaneous confidence intervals and sample size determination for multinomial proportions. Journal of the American Statistical Association, 90:366369.
Wald, A. Tests of statistical hypotheses concerning several parameters when the number of observations is large, Trans. Am. Math. Soc. 54 (1943) 426482.
Wilson, E. B. Probable inference, the law of succession and statistical inference, J.Am. Stat. Assoc. 22 (1927) 209212.
# Multinomial distribution with 3 classes, from which a sample of 79 elements
# were drawn: 23 of them belong to the first class, 12 to the
# second class and 44 to the third class. Punctual estimations
# of the probabilities from this sample would be 23/79, 12/79
# and 44/79 but we want to build 95% simultaneous confidence intervals
# for the true probabilities
MultinomCI(c(23, 12, 44), conf.level=0.95)
x < c(35, 74, 22, 69)
MultinomCI(x, method="goodman")
MultinomCI(x, method="sisonglaz")
MultinomCI(x, method="cplus1")
MultinomCI(x, method="wald")
MultinomCI(x, method="waldcc")
MultinomCI(x, method="wilson")
# compare to
BinomCI(x, n=sum(x))
# example in Goodman (1965)
MultinomCI(x = c(91,49,37,43),conf.level = 0.95,method="goodman")
# example from Sison, Glaz (1999) in Sangeetha (2013)  Table 2
x < c(56, 72, 73, 59, 62, 87, 58)
do.call(cbind, lapply(c("wald", "waldcc", "wilson",
"qh", "goodman", "fs", "sisonglaz"),
function(m) round(MultinomCI(x, method=m)[,1], 3)))