1 / 24

Applications of bootstrap method to finance

Applications of bootstrap method to finance. Chin-Ping King. Population distribution function F empirical distribution function(EDF) F n. F (x 1, x 2,…, x n ) where x = (x 1, x 2,…, x n )

berne
Download Presentation

Applications of bootstrap method to finance

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Applications of bootstrap method to finance Chin-Ping King

  2. Population distribution function Fempirical distribution function(EDF) Fn F (x1, x2,…,xn) where x= (x1, x2,…,xn) Fn(x1*, x2*,…,xn*) where x*= (x1*, x2*,…,xn*) Probability of elements of population which occur in x : P Probability of elements of EDF which occur in x*: Pn Pn~ Bi(P, (p*(1-p)/n)) By Weak Law of Large Number and Central Limit Theorem (n)1/2(Fn - F) d N(0, p*(1-p))

  3. Estimation of standard deviation and bias Estimator of θ : θ’ θ’ =s(x) Standard deviation: se={∑nj=1 [θ’j- θ’(.)]2/(n-1)} se(.)= ∑nj=1 θ’j/n Bias: bias=E[θ’]- θ Root mean square error of an estimator θ’ for θ: E[(θ’- θ)2]= se2*{1+(1/2)*(bias/se)} 2

  4. Nonparametric bootstrap The bootstrap algorithm for estimating standard errors(or bias) 1. Select B independent bootstrap samples x*1, x*2 , …, x*B, each consisting of n data drawn with replacement from x . Total possible number of distinct bootstrap samples is C(2n-1,n) . 2. Evaluate the bootstrap replication corresponding to each bootstrap sample θ’*(b)=s(x*b ) b=1,2,…,B 3. Estimate the standard error (or bias) by the sample standard deviation (or bias) of the B replications: se’B = se={∑Bb=1 [θ ’*(b)- θ’*(.)]2/(B-1)} θ’* (.)= ∑Bb=1 θ ’*(b)/B Bias’B=θ’*(.)- θ’

  5. A Schematic diagram of the nonparametric bootstrap Unknown Observed Random Empirical Bootstrap Population Sample Distribution Sample Distribution F x= (x1, x2,…,xn) Fn x*= (x1*, x2*,…,xn*) θ’ =s(x) θ’*(b)=s(x*b ) Statistic of interestBootstrap Replication

  6. Parametric bootstrap Function form of population probability distribution F has been known, but parameters in population probability distribution F are not known Parametric estimate of populationprobability distribution : Fpar We draw B samples of size n from the parametric estimate of estimate of the population probability distribution Fpar: Fparx*= (x1*, x2*,…,xn*)

  7. Error in bootstrap estimates mi= the ith moment of the bootstrap distribution of θ’ Var(se’B) = Var(m21/2 ) + E[(m2(△+2))/4B] △= m4/m22-3, the kurtosis of the bootstrap distribution of θ’ Var(m21/2 ):sample variation , it approaches zero as the sample size n approaches infinity E[(m2(△+2))/4B]:resampling variation, it approaches zero as B approaches infinity

  8. Confidence intervals based on bootstrap percentiles (1-α) Percentile interval: [θ’%low , θ’%up ]= [θ’*(α/2)B , θ’*(1-(α/2))B ] θ’*(α/2)B : 100*(α/2)th empirical percentile, or B *(α/2)th value in the ordered list of the B replications of θ’* θ’*(1-(α/2))B : 100*(1-(α/2))th empirical percentile, or B *(1-(α/2))th value in the ordered list of the B replications of θ’*

  9. Percentile interval lemma Suppose the transformation ψ’=t(θ’) perfectly normalize the distribution of θ’: ψ’~N(ψ, c2) For some standard deviation c. Then the percentile interval based on θ’ equals [t-1(ψ’-z(1-(α/2))*c), t-1(ψ’-z(α/2)*c)] Example: θ’ =exp(x) x ~ N(0,1) ψ’=t(θ’)=logθ’

  10. Coverage performance Results of 300 confidence interval realizations for θ’ =exp(x) Method % miss left % miss right Standard normal Interval 1.2 8.8 Bootstrap percentile Interval 4.8 5.2 miss left: left endpoint >1 Miss right: right endpoint <1

  11. Transformation-respecting property The percentile interval for any (monotone) parameter transformation ψ’=t(θ’) is simply the percentile interval for θ’ mapped by t(θ’) : [ψ’%low , ψ’%up ]= [t(θ ’%low) , t(θ ’%up) ]

  12. Better bootstrap confidence intervals (1-α) BCa interval: [θ ’low , θ ’up ]= [θ’*(α1) , θ’*(α2) ] α1and α2 are obtained by standard normal cumulative distribution function of some correction formulas for bootstrap replications. BCa interval is transformation respecting.

  13. Accuracy of bootstrap confidence interval For (1- α )coverage, approximate confidence interval points θ ’low and θ ’up are called first order accurate if: Pr(θ ≦θ ’low )= (α/2 )+ O(n-1/2) Pr(θ ≧θ ’up)= (α /2)+ O(n-1/2 ) And second order accurate if Pr(θ ≦θ ’low )= (α/2 )+ O(n-1) Pr(θ ≧θ ’up)= (α/2 )+ O(n-1) Percentile interval : first order accurate. BCa interval : second order accurate.

  14. Calibration of confidence interval points • Generate B bootstrap samples x*1, x*2 , …, x*B. For each sample b=1,2,…,B: 1a) Compute a λ-level confidence interval point θ’*λ (b) for a grid of values of λ. Where θ’*λ (b) can be θ’*(b)-z1-λ *se ’*(b) . 2. For each λ compute p’ (λ)=#{θ’ ≦ θ’*λ (b) }/B. 3. Find the value of λ satisfying p’ (λ)= α/2

  15. Calibration of percentile interval and BCainterval Once calibration of percentile interval: second order accurate Pr(θ ≦θ ’low )= (α/2 )+ O(n-1) Pr(θ ≧θ ’up)= (α/2 )+ O(n-1) Once calibration of BCa interval: third order accurate Pr(θ ≦θ ’low )= (α/2 )+ O(n-3/2) Pr(θ ≧θ ’up)= (α /2)+ O(n-3/2 )

  16. Computation of the bootstrap test statistics • Draw B samples of size n with replacement from x. • Evaluate ϕ(.) on each sample, ϕ(x*b) where ϕ(.) is test statistics b=1,2,…,B 3. Approximate P-value by P-value=#{ϕ(x*b) ≧ ϕobs}/B or P-value=#{ϕ(x*b) ≦ ϕobs}/B Where ϕobs= ϕ(x) the observed value of test statistics

  17. Asymptotic refinement Asymptotically normal test statistics ϕ ϕd N(0,σ2) ϕ ~ Gn(u,F) Gn(u,F): exact cumulative distribution Gn(u,F)=Pr(|ϕ| ≦u|F) Gn(u,F) φ(u) as n approaches infinity (assume σ=1) φ(u): standard normal cumulative distribution

  18. An asymptotic test is based on φ(u) φ(u)- Gn(u,F)=O(n-1) G*n(u):bootstrap cumulative distribution A bootstrap test is based on G*n(u) G*n(u)-Gn(u,F)= O(n-3/2)

  19. Reality test for data snooping Forecasting model: lk Benchmark model: l0 dk=lk- l0 H0 :maxk=1,2,…,nE(dk) ≦0 Data: 1000 daily closing stock prices of UMC Benchmark model: random walk with drift lnPt = a + lnPt-1 + εt

  20. Forecasting model : lnPt = a + ΔlnPt-1 + εt where ΔlnPt = lnPt – lnPt-1 V=(1/B) ∑Bb=1 d1(b) Quantile of bootstrap distribution Statistics V for V Critical value 0.000808 -1.9874504* The difference is significant, so reject H0 Forecasting model beat random walk model

  21. Inference when a nuisance parameter is not identified the null hypothesis Threshold Autoregressive (TAR)model: α10 + α11yt -1+ ε1tyt -1 ≦η yt = α20 + α21yt -1+ ε2tyt -1 > η η : threshold value H0:time series is linear H1: time series is TARprocess

  22. Data: monthly data of U.S. dollar/Sweden krona exchange rate from January 1974 to December 1998 U.S. dollar/Sweden krona Bootstrap P-value 0.0200 Reject H0 U.S. dollar/Sweden krona exchange rates follow TAR process

  23. Bootstrap percentile confidence interval

  24. Thanks for listening

More Related