Page 233 - Computational Statistics Handbook with MATLAB
P. 233

220                        Computational Statistics Handbook with MATLAB


                                2.  Sample  with replacement from the original  sample to get
                                          * b
                                               ,
                                           ,
                                   x * b  =  ( x 1 … x n * b  . )
                                3. Calculate the same statistic using the sample in step 2 to get the
                                   bootstrap replicates,  θ ˆ *b  .
                                4. Repeat steps 2 through 3, B times.
                                5. Using the bootstrap replicates, calculate  θ ˆ  *  .
                                                     θ
                                                     ˆ
                                6. Estimate the bias of   using Equation 6.17.
                             Example 6.10
                             We return to the forearm data of Example 6.9, where now we want to esti-
                             mate the bias in the sample skewness. We use the same bootstrap replicates
                             as before, so all we have to do is to calculate the bias using Equation 6.17.
                                % Use the same replicates from before.
                                % Evaluate the mean using Equation 6.15.
                                meanb = mean(thetab);
                                % Now estimate the bias using Equation 6.17.
                                biasb = meanb - theta;
                             We have an estimated bias of -0.011. Note that this is small relative to the stan-
                             dard error.

                              In the next chapter, we discuss another method for estimating the bias and
                             the standard error of a statistic called the jackknife. The jackknife method is
                             related to the bootstrap. However, since it is based on the reuse or partition-
                             ing of the original sample rather than resampling, we do not include it here.



                                      onfide
                                       onfiden
                             Bootstrap  CC onfideonfide  nnceInteInte  rvval  als  ss s
                                      C
                             BootstrapC
                                              ce
                                             nceInteceInte
                             BootstrapBootstrap
                                                    rr vvalal
                                               r
                             There are several ways of constructing confidence intervals using the boot-
                             strap. We discuss three of them here: the standard interval, the bootstrap-t
                             interval and the percentile method. Because it uses the jackknife procedure,
                                                                                 will be presented
                             an improved bootstrap confidence interval called the  BC a
                             in the next chapter.
                                                   ddeenceIntervaInterva
                                                   ii
                              a
                                                     nce
                             BootstrapStand
                             BootstrapStandStand
                                            d
                                            rdConfdConfi
                             Bootstrap
                             Bootstrap Stand  ar aarr dConfConf  idde  enceIntervanceInterva l  ll l
                             The bootstrap standard confidence interval is based on the parametric form
                             of the confidence interval that was discussed in Section 6.2. We showed that
                                      ⋅
                             the  1 –(  α) 100%   confidence interval for the mean can be found using
                                                (  1 – α 2) σ     ( α 2) σ 
                                                                    ⁄
                                                     ⁄
                                           P X –  z   ------- <  µ <  X –  z  -------  =  1 – α  .  (6.20)
                                                       n               n 
                             © 2002 by Chapman & Hall/CRC
   228   229   230   231   232   233   234   235   236   237   238