<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJS</journal-id><journal-title-group><journal-title>Open Journal of Statistics</journal-title></journal-title-group><issn pub-type="epub">2161-718X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojs.2015.56052</article-id><article-id pub-id-type="publisher-id">OJS-60252</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Bootstrap Confidence Intervals for Proportions of Unequal Sized Groups Adjusted for Overdispersion
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>livia</surname><given-names>Wanjeri Mwangi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ali</surname><given-names>Islam</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Orawo</surname><given-names>Luke</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Mathematics, Egerton University, Njoro, Kenya</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>owanjeri@yahoo.com(LWM)</email>;<email>asislam54@yahoo.com(AI)</email>;<email>orawo2000@yahoo.com(OL)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>13</day><month>10</month><year>2015</year></pub-date><volume>05</volume><issue>06</issue><fpage>502</fpage><lpage>510</lpage><history><date date-type="received"><day>17</day>	<month>August</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>8</month>	<year>October</year>	</date><date date-type="accepted"><day>13</day>	<month>October</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   Group testing is a method of pooling a number of units together and performing a single test on the resulting group. It is an appealing option when few individual units are thought to be infected leading to reduced costs of testing as compared to individually testing the units. Group testing aims to identify the positive groups in all the groups tested or to estimate the proportion of positives (p) in a population. Interval estimation methods of the proportions in group testing for unequal group sizes adjusted for overdispersion have been examined. Lately improvement in statistical methods allows the construction of highly accurate confidence intervals (CIs). The aim here is to apply group testing for estimation and generate highly accurate Bootstrap confidence intervals (CIs) for the proportion of defective or positive units in particular. This study provided a comparison of several proven methods of constructing CIs for a binomial proportion after adjusting for overdispersion in group testing with groups of unequal sizes. Bootstrap resampling was applied on data simulated from binomial distribution, and confidence intervals with high coverage probabilities were produced. This data was assumed to be overdispersed and independent between groups but correlated within these groups. Interval estimation methods based on the Wald, the Logit and Complementary log-log (CLL) functions were considered. The criterion used in the comparisons is mainly the coverage probabilities attained by nominal 95% CIs, though interval width is also regarded. Bootstrapping produced CIs with high coverage probabilities for each of the three interval methods. 
 
</p></abstract><kwd-group><kwd>Group Testing</kwd><kwd> Overdispersion</kwd><kwd> Quasi-Likelihood</kwd><kwd> Confidence Interval</kwd><kwd> Bootstrapping</kwd><kwd> Coverage Probability</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Group testing started or originated with [<xref ref-type="bibr" rid="scirp.60252-ref1">1</xref>] during World War II as a frugal method of testing blood specimens of army inductees in order to detect the presence of infection. Group testing has been applied in many areas: plant disease assessment [<xref ref-type="bibr" rid="scirp.60252-ref2">2</xref>] , fisheries [<xref ref-type="bibr" rid="scirp.60252-ref3">3</xref>] , quality control, drug discovery and transmission of viruses by insect vectors [<xref ref-type="bibr" rid="scirp.60252-ref4">4</xref>] . It has also been used to screen the population for the presence of HIV/Aids antibody [<xref ref-type="bibr" rid="scirp.60252-ref5">5</xref>] . [<xref ref-type="bibr" rid="scirp.60252-ref1">1</xref>] proposed that, rather than testing each blood specimen individually, portions each of k specimens can be pooled and the pooled specimen can be tested. If the pooled specimen is free of infection, all k inductees are passed with no further tests, otherwise the remaining portions of each of the blood specimens are to be tested individually. If the prevalence of infection is low, the expected number of tests per inductee and thus the expected cost per inductee, would be reduced. [<xref ref-type="bibr" rid="scirp.60252-ref1">1</xref>] assumed that tests were perfect, that is, a negative reading indicated the group contained no defective item and a positive reading indicated the presence of at least one defective item.</p><p>Group tests save resources since many units are tested without testing them individually. One of their advantages as a method of estimation is they are time efficient. In most of the studies carried out, the experimental unit/group is a litter. It has been observed that there is a tendency of littermates to respond more alike than animals from different litters, the “litter effect”. This litter effect is also known as the extra-dispersion (over/under- dispersion) or the intra-litter correlation. These litters may be of equal or unequal sizes. The concern here is with methods of establishing the confidence interval for the proportion of defective units p in unequal groups with adjustment for overdispersion using bootstrap technique.</p><p>Overdispersion is the phenomenon of having greater variability than predicted by the random component of the model; this is common in the modeling of binomial distribution for group testing [<xref ref-type="bibr" rid="scirp.60252-ref6">6</xref>] . Overdispersion is a very common feature in applied data analysis because naturally populations are frequently heterogeneous contrary to the assumptions implicit within widely used simple parametric models. The over-dispersion parameter is common in proportions in biology, toxicology, medicine, genetics and other similar fields and is important in making inference regarding the regression parameters on the mean [<xref ref-type="bibr" rid="scirp.60252-ref7">7</xref>] .</p><p>Maximum likelihood estimation gives a unified approach to estimation, which is well-defined, in binomial distribution and many other problems. The maximum likelihood as an estimator has been studied and seconded as an approach for using the proposed extended Beta-Binomial (BB) model to analyze over/under-dispersed proportions [<xref ref-type="bibr" rid="scirp.60252-ref8">8</xref>] . They estimated this parameter using maximum likelihood estimation (MLE). This concurred with the observation that the MLEs may be biased when the sample size n or the total Fisher information is small [<xref ref-type="bibr" rid="scirp.60252-ref9">9</xref>] . This bias is usually ignored in practice, the justification being that it is small compared with the standard errors.</p><p>Studies on point estimation in terms of bias and efficiency and the test for presence of overdispersion for both counts and proportions data have been done [<xref ref-type="bibr" rid="scirp.60252-ref10">10</xref>] . Interval estimation methods in which groups were of unequal sizes were carried out using interval based on the MLE, the Logit and CLL [<xref ref-type="bibr" rid="scirp.60252-ref11">11</xref>] . The quasi-likelihood approach to correct for overdispersion is used by [<xref ref-type="bibr" rid="scirp.60252-ref6">6</xref>] , he examines interval estimation methods based on functions of the MLE, the Logit and CLL functions.</p></sec><sec id="s2"><title>2. Literature Review</title><sec id="s2_1"><title>2.1. Bootstrapping</title><p>The bootstrap method was introduced by [<xref ref-type="bibr" rid="scirp.60252-ref12">12</xref>] . It is a very general resampling procedure for estimating the distributions of statistics based on independent observations. It has spread fast in statistical sciences within the last decades. The primary task of a statistician is to summarize a sample based study and generalize the finding to the parent population in a scientific manner. Of course, a summary statistic like the sample mean will fluctuate from sample to sample and a statistician would like to know the magnitude of these fluctuations around the corresponding population parameter.</p><p>To understand bootstrap, suppose it were possible to draw repeated samples (of the same size) from the population of interest, a large number of times. Then, one would get a fairly good idea about the sampling distribution of a particular statistic from the collection of its values arising from these repeated samples. But, that does not make sense as it would be too expensive and defeat the purpose of a sample study. The purpose of a sample study is to gather information cheaply in a timely fashion. The idea behind bootstrap is to use the data of a sample study at hand as a “surrogate population”, for the purpose of approximating the sampling distribution of a statistic, that is, to resample (with replacement) from the sample data at hand and create a large number of “phantom samples” known as bootstrap samples. These samples can be used to obtain more improved estimation of unknown parameter (s) of a probability model.</p><p>Recently, study has been done on the construction of bootstrap confidence intervals for the overdispersion parameter in equal proportions in Beta Binomial using Maximum Likelihood Estimator, Method of Moments Estimator and Quasi-likelihood [<xref ref-type="bibr" rid="scirp.60252-ref13">13</xref>] . Bootstrapping has been seen to provide improved CIs and is useful even when very little is known about the underlying distributions [<xref ref-type="bibr" rid="scirp.60252-ref14">14</xref>] .</p></sec><sec id="s2_2"><title>2.2. Estimation of Proportions</title><sec id="s2_2_1"><title>2.2.1. Confidence Intervals Based on the Maximum Quasi-Likelihood Estimate</title><p>Overdispersion causes one to underestimate the variance of parameter estimates. A quasi-likelihood approach can be employed to correct for the overdispersion phenomenon which occurs with binary data [<xref ref-type="bibr" rid="scirp.60252-ref15">15</xref>] . For binary data, an alternative method for handling overdispersion is to use models for which a binomial parameter itself has a beta distribution [<xref ref-type="bibr" rid="scirp.60252-ref16">16</xref>] . However, the MLE of the beta-binomial method has been shown to be less flexible than that of the quasi-likelihood approach, and so cannot be recommended for general use. For comparisons of CI methods in previous literature, the overdispersion phenomenon needs to be taken into account for constructing the CI of the proportion. Therefore, the objective is to use quasi-likelihood methodology to evaluate several asymptotic CI methods and construct bootstrap CI based on these methods for groups of different sizes, adjusted for overdispersion in group testing.</p><p>Suppose for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x5.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x6.png" xlink:type="simple"/></inline-formula> groups of size<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x7.png" xlink:type="simple"/></inline-formula>, are tested, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x8.png" xlink:type="simple"/></inline-formula> of these groups test positive. The binomial parameter for the distribution of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x9.png" xlink:type="simple"/></inline-formula>, which is the probability the groups tests positive is, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x10.png" xlink:type="simple"/></inline-formula>, where p is the probability an individual unit tests positive. Therefore considering overdispersion, the expected value and the variance of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x11.png" xlink:type="simple"/></inline-formula> are</p><disp-formula id="scirp.60252-formula332"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x12.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.60252-formula333"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x13.png"  xlink:type="simple"/></disp-formula><p>When groups of unequal size occur and overdispersion occurs together, the variance of the quasi-score function of p is</p><disp-formula id="scirp.60252-formula334"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x14.png"  xlink:type="simple"/></disp-formula><p>The maximum quasi-likelihood estimate of p denoted as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x15.png" xlink:type="simple"/></inline-formula>, is the solution to</p><disp-formula id="scirp.60252-formula335"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x16.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.60252-formula336"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x17.png"  xlink:type="simple"/></disp-formula><p>And</p><disp-formula id="scirp.60252-formula337"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x18.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x19.png" xlink:type="simple"/></inline-formula>is a variance function, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x20.png" xlink:type="simple"/></inline-formula> is the generalized Pearson statistic.</p><p>Under the usual limiting conditions on the<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x21.png" xlink:type="simple"/></inline-formula>, the asymptotic variance of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x22.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.60252-formula338"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x23.png"  xlink:type="simple"/></disp-formula><p>Then the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x24.png" xlink:type="simple"/></inline-formula> Wald type CI for p is <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x25.png" xlink:type="simple"/></inline-formula></p><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x26.png" xlink:type="simple"/></inline-formula> is the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x27.png" xlink:type="simple"/></inline-formula> quantile of the standard normal distribution.</p></sec><sec id="s2_2_2"><title>2.2.2. Confidence Intervals Based on the Logit</title><p>Asymptotic-likelihood methods to construct the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x28.png" xlink:type="simple"/></inline-formula> confidence intervals of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x29.png" xlink:type="simple"/></inline-formula> using the Logit trans- formations were presented [<xref ref-type="bibr" rid="scirp.60252-ref11">11</xref>] .</p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x30.png" xlink:type="simple"/></inline-formula> be the Fisher information that Y contains about the parameter<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x31.png" xlink:type="simple"/></inline-formula>. If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x32.png" xlink:type="simple"/></inline-formula> and h is differen-</p><p>tiable, then the information that Y contains about <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x33.png" xlink:type="simple"/></inline-formula> is <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x34.png" xlink:type="simple"/></inline-formula></p><p>For the logit transformation</p><disp-formula id="scirp.60252-formula339"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x35.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.60252-formula340"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x36.png"  xlink:type="simple"/></disp-formula><p>Thus the information of logit (p) is</p><disp-formula id="scirp.60252-formula341"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x37.png"  xlink:type="simple"/></disp-formula><p>After adjusting for overdispersion, the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x38.png" xlink:type="simple"/></inline-formula> approximate CI for p is</p><disp-formula id="scirp.60252-formula342"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x39.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.60252-formula343"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x40.png"  xlink:type="simple"/></disp-formula></sec><sec id="s2_2_3"><title>2.2.3. Confidence Intervals Based on the Complementary Log Log (CLL)</title><p>Hepworth also presented the CLL parameter transformations.</p><p>Suppose</p><disp-formula id="scirp.60252-formula344"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x41.png"  xlink:type="simple"/></disp-formula><p>And</p><disp-formula id="scirp.60252-formula345"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x42.png"  xlink:type="simple"/></disp-formula><p>The information of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x43.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.60252-formula346"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x44.png"  xlink:type="simple"/></disp-formula><p>Then the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x45.png" xlink:type="simple"/></inline-formula> approximate CI for p after allowing for overdispersion is given by</p><disp-formula id="scirp.60252-formula347"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-1240558x46.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x47.png" xlink:type="simple"/></inline-formula> (17)</p></sec></sec><sec id="s2_3"><title>2.3. Choice of Group Sizes</title><p>A lot of attention has been devoted in studies to problems involving the group sizes, k. The question being what values should be chosen for k, the number of units in each group. If k is too large, π is close to 1 and all groups are likely to test positive, also if k is too small, π is closer to 0 [<xref ref-type="bibr" rid="scirp.60252-ref17">17</xref>] . When k is large group testing becomes uninformative and expensive. Therefore it is suitable to choose group sizes that produce some positives and some negatives. We therefore ensured that the sample sizes used yielded the desirable level, that is, margin of error is minimized.</p></sec><sec id="s2_4"><title>2.4. Interval Width and Coverage Probability</title><p>Each Interval was examined for its Interval width and coverage probability. Confidence Intervals based on MLE ordering and Sterne’s technique for a dilution assay with 64 outcomes was compared [<xref ref-type="bibr" rid="scirp.60252-ref18">18</xref>] , those with the shorter widths were considered superior. The width of the intervals is considered an important criterion in comparing CI methods. Interval widths of the intervals were examined since the three methods; Wald, Logit and CLL do not differ a great deal with respect to coverage.</p><p>Coverage probability gives the proportion of the time that the interval contains p. Bootstrapping enabled us to estimate the coverage of the CIs with respect to the number of hypothetical repetitions of the entire procedure. The nominal coverage will be set at 0.95. The actual coverage of the three methods will be compared to the nominal coverage, greater coverage than the nominal is preferable and hence conservative.</p></sec></sec><sec id="s3"><title>3. Simulation</title><p>First, positive groups, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x48.png" xlink:type="simple"/></inline-formula>are simulated from the number of groups, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x49.png" xlink:type="simple"/></inline-formula>with their respective sizes <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x50.png" xlink:type="simple"/></inline-formula> for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x51.png" xlink:type="simple"/></inline-formula>. This is followed by bootstrapping Y<sub>i</sub>’s and then computing the estimate of p. Coverage probabilities of the three interval methods for the estimate of p for 500 bootstrap simulations are obtained, and the program is run repetitively in order to see how the three method’s coverage probabilities vary.</p><p>The nominal coverage probability is set at 0.95.</p><p>This is then compared to estimating p from the same data that is the Y<sub>i</sub>’s without bootstrapping and the program run repetitively so as to calculate the coverage probability.</p><p>The number of bootstraps B’s were varied for the same combination of p, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x52.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x53.png" xlink:type="simple"/></inline-formula>. Coverage probabilities are generated. This was to assess the performance of the coverage probabilities as the number of bootstrap simulations increase for each method.</p><p>Next p was varied for the same N<sub>i</sub>, k<sub>i</sub> and the same number of bootstrap simulations (500), and coverage probabilities are obtained.</p><p>Lastly, interval widths for the data are calculated for bootstrapping technique for the three interval methods and compared to those intervals got from each of the three methods without bootstrapping. This is in order to investigate whether bootstrapping gives more precise intervals.</p></sec><sec id="s4"><title>4. Discussion</title><p>Group testing was primarily used to inspect individual members from a large population [<xref ref-type="bibr" rid="scirp.60252-ref1">1</xref>] . Therefore probability used in this study ranged from 0.01 to 0.2. From <xref ref-type="table" rid="table1"><xref ref-type="table" rid="table">Table </xref>1</xref> it was seen that the CLL and the Logit gave more conservative coverage probabilities consistently for different bootstrap simulations for p = 0.02. All the three methods differ only slightly with CLL and Logit sometimes taking the same values of coverage probability. The coverage probabilities oscillate about the nominal level with Wald being less conservative than CLL and Logit. This was then compared to coverage probabilities of the three methods for the same data combination without applying bootstrapping technique. Similarly CLL outperformed the Wald interval method and at times the Logit Interval method as seen in <xref ref-type="table" rid="table2"><xref ref-type="table" rid="table">Table </xref>2</xref>. It is clearly seen that without bootstrapping technique, CIs generated have lower coverage probabilities than CIs generated after applying bootstrapping.</p><p>In <xref ref-type="table" rid="table3"><xref ref-type="table" rid="table">Table </xref>3</xref> and <xref ref-type="fig" rid="fig1">Figure 1</xref>, the number of bootstrap simulations was varied. It was seen that as the number of</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1"><xref ref-type="table" rid="table">Table </xref>1</xref></label><caption><title> Coverage probabilities for the three Interval methods for 500 simulations for p = 0.02</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Method</th><th align="center" valign="middle" >Wald</th><th align="center" valign="middle" >CLL</th><th align="center" valign="middle" >Logit</th></tr></thead><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.938</td><td align="center" valign="middle" >0.9461</td><td align="center" valign="middle" >0.946</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.954</td><td align="center" valign="middle" >0.948</td><td align="center" valign="middle" >0.948</td></tr><tr><td align="center" valign="middle" >p=0.02</td><td align="center" valign="middle" >0.942</td><td align="center" valign="middle" >0.9441</td><td align="center" valign="middle" >0.943</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.938</td><td align="center" valign="middle" >0.952</td><td align="center" valign="middle" >0.952</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.95</td><td align="center" valign="middle" >0.958</td><td align="center" valign="middle" >0.958</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.932</td><td align="center" valign="middle" >0.9481</td><td align="center" valign="middle" >0.948</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2"><xref ref-type="table" rid="table">Table </xref>2</xref></label><caption><title> Coverage probabilities for data for 500 simulations without bootstrapping for p = 0.02</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Method</th><th align="center" valign="middle" >Wald</th><th align="center" valign="middle" >CLL</th><th align="center" valign="middle" >Logit</th></tr></thead><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.812</td><td align="center" valign="middle" >0.862</td><td align="center" valign="middle" >0.862</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.822</td><td align="center" valign="middle" >0.873</td><td align="center" valign="middle" >0.87</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.832</td><td align="center" valign="middle" >0.912</td><td align="center" valign="middle" >0.91</td></tr><tr><td align="center" valign="middle" >p = 0.02</td><td align="center" valign="middle" >0.81</td><td align="center" valign="middle" >0.852</td><td align="center" valign="middle" >0.85</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.802</td><td align="center" valign="middle" >0.844</td><td align="center" valign="middle" >0.844</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.804</td><td align="center" valign="middle" >0.812</td><td align="center" valign="middle" >0.81</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3"><xref ref-type="table" rid="table">Table </xref>3</xref></label><caption><title> Coverage probabilities for the three methods for varying bootstrap simulations (B) for p = 0.02</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Number of</th><th align="center" valign="middle" ></th><th align="center" valign="middle" ></th><th align="center" valign="middle" ></th></tr></thead><tr><td align="center" valign="middle" >Simulations (B)</td><td align="center" valign="middle" >Wald</td><td align="center" valign="middle" >CLL</td><td align="center" valign="middle" >Logit</td></tr><tr><td align="center" valign="middle" >200</td><td align="center" valign="middle" >0.937</td><td align="center" valign="middle" >0.945</td><td align="center" valign="middle" >0.94</td></tr><tr><td align="center" valign="middle" >400</td><td align="center" valign="middle" >0.936</td><td align="center" valign="middle" >0.9486</td><td align="center" valign="middle" >0.948</td></tr><tr><td align="center" valign="middle" >600</td><td align="center" valign="middle" >0.9283</td><td align="center" valign="middle" >0.9417</td><td align="center" valign="middle" >0.9416</td></tr><tr><td align="center" valign="middle" >800</td><td align="center" valign="middle" >0.935</td><td align="center" valign="middle" >0.9512</td><td align="center" valign="middle" >0.951</td></tr><tr><td align="center" valign="middle" >1000</td><td align="center" valign="middle" >0.932</td><td align="center" valign="middle" >0.945</td><td align="center" valign="middle" >0.945</td></tr><tr><td align="center" valign="middle" >1500</td><td align="center" valign="middle" >0.94</td><td align="center" valign="middle" >0.951</td><td align="center" valign="middle" >0.952</td></tr><tr><td align="center" valign="middle" >2000</td><td align="center" valign="middle" >0.944</td><td align="center" valign="middle" >0.9556</td><td align="center" valign="middle" >0.9556</td></tr><tr><td align="center" valign="middle" >2500</td><td align="center" valign="middle" >0.9525</td><td align="center" valign="middle" >0.9565</td><td align="center" valign="middle" >0.956</td></tr><tr><td align="center" valign="middle" >3000</td><td align="center" valign="middle" >0.951</td><td align="center" valign="middle" >0.9562</td><td align="center" valign="middle" >0.9561</td></tr><tr><td align="center" valign="middle" >3500</td><td align="center" valign="middle" >0.93514</td><td align="center" valign="middle" >0.95457</td><td align="center" valign="middle" >0.9544</td></tr><tr><td align="center" valign="middle" >4000</td><td align="center" valign="middle" >0.93425</td><td align="center" valign="middle" >0.94525</td><td align="center" valign="middle" >0.9451</td></tr><tr><td align="center" valign="middle" >5000</td><td align="center" valign="middle" >0.9416</td><td align="center" valign="middle" >0.9556</td><td align="center" valign="middle" >0.9554</td></tr><tr><td align="center" valign="middle" >7000</td><td align="center" valign="middle" >0.9405</td><td align="center" valign="middle" >0.9521</td><td align="center" valign="middle" >0.952</td></tr><tr><td align="center" valign="middle" >10,000</td><td align="center" valign="middle" >0.9512</td><td align="center" valign="middle" >0.967</td><td align="center" valign="middle" >0.967</td></tr></tbody></table></table-wrap><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> The figure shows a plot of coverage probabilities for the three methods for varying bootstrap simulations for p = 0.02 and for a nominal coverage probability of 0.95</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-1240558x54.png"/></fig><p>bootstrap simulations increased, all the three methods became more precise, that is, their coverage probabilities hovered slightly above and below the nominal coverage probability.</p><p>Group testing works best when the prevalence rate or the proportion is sufficiently small. According to <xref ref-type="table" rid="table4"><xref ref-type="table" rid="table">Table </xref>4</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref>, we see that as the proportion increases to 0.1 the coverage probability becomes very low for all the three interval methods. Between 0.03 and 0.08, the coverage probabilities for the three methods are conservative.</p><p>From <xref ref-type="table" rid="table5"><xref ref-type="table" rid="table">Table </xref>5</xref> and <xref ref-type="table" rid="table6"><xref ref-type="table" rid="table">Table </xref>6</xref>, it is seen that the interval widths for both CLL and Logit are not differing a great deal for the same bootstrap simulations. CLL has shorter widths than logit. However, interval widths for the Wald for the three bootstrap simulations are considerably wider than those of the other two interval methods. When interval widths are compared without bootstrap simulations to those interval widths with bootstrap simulations we note that they are wider.</p><p>According to <xref ref-type="table" rid="table4"><xref ref-type="table" rid="table">Table </xref>4</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref>, all the three methods perform well as the number of bootstrap simulations increases, with the coverage probabilities of the CLL and Logit being oftentimes greater than the nominal coverage probability.</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> The figure above is a plot of coverage probabilities for the three methods for 500 bootstrap simulations for varying probability, that is, 0.01 &lt; p &lt; 0.2</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/3-1240558x55.png"/></fig><table-wrap id="table4" ><label><xref ref-type="table" rid="table4"><xref ref-type="table" rid="table">Table </xref>4</xref></label><caption><title> Coverage probabilities for varying p for the same 500 simulations for a nominal coverage of 0.95</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" ></th><th align="center" valign="middle" >Coverage Probability</th><th align="center" valign="middle" ></th></tr></thead><tr><td align="center" valign="middle" >Probability</td><td align="center" valign="middle" >Wald</td><td align="center" valign="middle" >CLL</td><td align="center" valign="middle" >Logit</td></tr><tr><td align="center" valign="middle" >0.01</td><td align="center" valign="middle" >0.904</td><td align="center" valign="middle" >0.922</td><td align="center" valign="middle" >0.922</td></tr><tr><td align="center" valign="middle" >0.03</td><td align="center" valign="middle" >0.92</td><td align="center" valign="middle" >0.934</td><td align="center" valign="middle" >0.93</td></tr><tr><td align="center" valign="middle" >0.04</td><td align="center" valign="middle" >0.982</td><td align="center" valign="middle" >0.986</td><td align="center" valign="middle" >0.984</td></tr><tr><td align="center" valign="middle" >0.06</td><td align="center" valign="middle" >0.988</td><td align="center" valign="middle" >0.99</td><td align="center" valign="middle" >0.99</td></tr><tr><td align="center" valign="middle" >0.08</td><td align="center" valign="middle" >0.99</td><td align="center" valign="middle" >0.992</td><td align="center" valign="middle" >0.99</td></tr><tr><td align="center" valign="middle" >0.1</td><td align="center" valign="middle" >0.988</td><td align="center" valign="middle" >0.992</td><td align="center" valign="middle" >0.99</td></tr><tr><td align="center" valign="middle" >0.2</td><td align="center" valign="middle" >0.842</td><td align="center" valign="middle" >0.857</td><td align="center" valign="middle" >0.857</td></tr></tbody></table></table-wrap><table-wrap id="table5" ><label><xref ref-type="table" rid="table5"><xref ref-type="table" rid="table">Table </xref>5</xref></label><caption><title> <xref ref-type="table" rid="table">Table </xref>of interval widths for the three interval methods for data without bootstrapping</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Interval Method</th><th align="center" valign="middle" >Y = (2,2,4,3,3)</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x56.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ></th></tr></thead><tr><td align="center" valign="middle" >Upper Confidence Limit</td><td align="center" valign="middle" >Lower Confidence Limit</td><td align="center" valign="middle" >Interval Width</td></tr><tr><td align="center" valign="middle" >Wald</td><td align="center" valign="middle" >0.02654614</td><td align="center" valign="middle" >0.01587489</td><td align="center" valign="middle" >0.01067125</td></tr><tr><td align="center" valign="middle" >CLL</td><td align="center" valign="middle" >0.02726701</td><td align="center" valign="middle" >0.01648786</td><td align="center" valign="middle" >0.01077915</td></tr><tr><td align="center" valign="middle" >Logit</td><td align="center" valign="middle" >0.02725691</td><td align="center" valign="middle" >0.01648266</td><td align="center" valign="middle" >0.01077426</td></tr></tbody></table></table-wrap><table-wrap id="table6" ><label><xref ref-type="table" rid="table6"><xref ref-type="table" rid="table">Table </xref>6</xref></label><caption><title> <xref ref-type="table" rid="table">Table </xref>of interval widths for the three interval methods after bootstrapping</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Interval Method</th><th align="center" valign="middle" >Y = (2,2,4,3,3)</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1240558x57.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" ></th></tr></thead><tr><td align="center" valign="middle" >Upper Confidence Limit</td><td align="center" valign="middle" >Lower Confidence Limit</td><td align="center" valign="middle" >Interval Width</td></tr><tr><td align="center" valign="middle" >Wald</td><td align="center" valign="middle" >0.02673085</td><td align="center" valign="middle" >0.016212005</td><td align="center" valign="middle" >0.0106108</td></tr><tr><td align="center" valign="middle" >CLL</td><td align="center" valign="middle" >0.02743541</td><td align="center" valign="middle" >0.01672070</td><td align="center" valign="middle" >0.01071471</td></tr><tr><td align="center" valign="middle" >Logit</td><td align="center" valign="middle" >0.02742548</td><td align="center" valign="middle" >0.01671553</td><td align="center" valign="middle" >0.01070995</td></tr></tbody></table></table-wrap></sec><sec id="s5"><title>5. Conclusion</title><p>Bootstrap methods produce good confidence intervals by an order of magnitude upon the accuracy of standard intervals [<xref ref-type="bibr" rid="scirp.60252-ref19">19</xref>] . In this paper we implement bootstrapping to the three interval methods: Wald, CLL and Logit. The comparison results in this study show that bootstrapping produces more precise confidence intervals than standard intervals for all the three procedures. Likewise the coverage probabilities from the three methods are consistent with the CLL being most conservative. Wald is more liberal. As the number of bootstrap simulations increases the more conservative the coverage becomes.</p></sec><sec id="s6"><title>Cite this paper</title><p>Olivia WanjeriMwangi,AliIslam,OrawoLuke, (2015) Bootstrap Confidence Intervals for Bootstrap Confidence Intervals for Adjusted for Overdispersion. Open Journal of Statistics,05,502-510. doi: 10.4236/ojs.2015.56052</p></sec></body><back><ref-list><title>References</title><ref id="scirp.60252-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Dorfman, R. (1943) The Detection of Defective Members of Large Populations. The Annals of Mathematical Statistics, 14, 436-440. &lt;/br&gt;http://dx.doi.org/10.1214/aoms/1177731363</mixed-citation></ref><ref id="scirp.60252-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Fletcher, J., Russell, A. and Butler, R. (1999) Seed-Borne Cucumber Mosaic Virus in New Zealand Lentil Crops: Yield Effects and Disease Incidence. New Zealand Journal of Crop and Horticulture Science, 27, 197-204.&lt;/br&gt;http://dx.doi.org/10.1080/01140671.1999.9514097</mixed-citation></ref><ref id="scirp.60252-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Worlund, D. and Taylor, G. (1983) Estimation of Disease Incidence in Fish Populations. Canadian Journal of Fisheries and Aquatic Science, 40, 2194-2197. &lt;/br&gt;http://dx.doi.org/10.1139/f83-254</mixed-citation></ref><ref id="scirp.60252-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Swallow, W.H. (1985) Group Testing for Estimating Infection Rates and Probabilities of Disease Transmission. Phytopathology, 75, 882-889. &lt;/br&gt;http://dx.doi.org/10.1094/phyto-75-882</mixed-citation></ref><ref id="scirp.60252-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Kline, R.L., Bothus, T., Brookmeyer, R., Zeyer, S. and Quinn, T. (1989) Evaluation of Human Immunodeficiency Virus Seroprevalence in Population Surveys Using Pooled Sera. Journal of Clinical Microbiology, 27, 1449-1452.</mixed-citation></ref><ref id="scirp.60252-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Liu, S.-C., Chiang, K.-S., Lin, C.-H. and Deng, T.-C. (2010) Confidence Interval Procedures for Proportions Estimated by Group Testing with Groups of Unequal Size Adjusted for Overdispersion. Journal of Applied Statistics, 38, 1467-1482. &lt;/br&gt;http://dx.doi.org/10.1080/02664763.2010.505953</mixed-citation></ref><ref id="scirp.60252-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Paul, S. and Islam, A.S. (1998) Joint Estimation of the Mean and Dispersion Parameters in the Analysis of Proportions: A Comparison of Efficiency and Bias. Canadian Journal of Statistics, 26, 83-94. &lt;/br&gt;http://dx.doi.org/10.2307/3315675</mixed-citation></ref><ref id="scirp.60252-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Saha, K.K. and Paul, S.R. (2005) Bias-Corrected Maximum Likelihood Estimator of the Negative Binomial Dispersion Parameter. Biometrics, 61, 179-185. &lt;/br&gt;http://dx.doi.org/10.1111/j.0006-341X.2005.030833.x</mixed-citation></ref><ref id="scirp.60252-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Williams, D.A. (1975) The Analysis of Binary Responses for Toxicological Experiments Involving Reproduction and Teratogenicity. Biometrics, 31, 946-952. &lt;/br&gt;http://dx.doi.org/10.2307/2529820</mixed-citation></ref><ref id="scirp.60252-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Saha, K.K. (2008) Semi-Parametric Estimation for the Dispersion Parameter in the Analysis of over- or Underdispersed Count Data. Journal of Applied Statistics, 35, 1383-1397. &lt;/br&gt;http://dx.doi.org/10.1080/02664760802382459</mixed-citation></ref><ref id="scirp.60252-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Hepworth, G. (2005) Confidence Intervals for Proportions Estimated by Group Testing with Groups of Unequal Size. Journal of Agricultural, Biological, and Environmental Statistics, 10, 478-497. &lt;/br&gt;http://dx.doi.org/10.1198/108571105X81698</mixed-citation></ref><ref id="scirp.60252-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Efron, B. (1979) Bootstrap Methods: Another Look at the Jackknife. The Annals of Statistics, 7, 1-26.&lt;/br&gt;http://dx.doi.org/10.1214/aos/1176344552</mixed-citation></ref><ref id="scirp.60252-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Siele, D. (2012) Interval Estimation for the Beta Binomial Dispersion Parameter. MSc Thesis, Egerton University, Njoro.</mixed-citation></ref><ref id="scirp.60252-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Saha, K. and Sen, D. (2009) Improved Confidence Interval for the Dispersion Parameter in Count Data Using Profile Likelihood. Technical Report No. 4/09.</mixed-citation></ref><ref id="scirp.60252-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Williams, D.A. (1982) Extra-Binomial Variation in Logistic Linear Models. Journal of the Royal Statistical Society. Series C (Applied Statistics), 31, 144-148. &lt;/br&gt;http://dx.doi.org/10.2307/2347977</mixed-citation></ref><ref id="scirp.60252-ref16"><label>16</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Crowder</surname><given-names> M.J. </given-names></name>,<etal>et al</etal>. (<year>1978</year>)<article-title>Beta-Binomial Anova for Proportions. Applied Statistics, 27, 34-37</article-title><source> Mathematical Statistics</source><volume> 14</volume>,<fpage> 436</fpage>-<lpage>440</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.60252-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Hepworth, G. (1996) Exact Confidence Intervals for Proportions Estimated by Group Testing. Biometrics, 52, 1134-1146.&lt;/br&gt; http://dx.doi.org/10.2307/2533075</mixed-citation></ref><ref id="scirp.60252-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Loyer, M.W. and Hamilton, M.A. (1984) Interval Estimation of the Density of Organisms Using a Serial-Dilution Experiment. Biometrics, 40, 907-916. &lt;/br&gt;http://dx.doi.org/10.2307/2531142</mixed-citation></ref><ref id="scirp.60252-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Efron, B. and DiCiccio, T. (1996) Bootstrap Confidence Intervals. Statistical Science, 11, 189-228.&lt;/br&gt;http://dx.doi.org/10.1214/ss/1032280214</mixed-citation></ref></ref-list></back></article>