Hello,

I'm currently working on a paper, and for one section of it need I BCa confidence intervals.

My goal is to do exactly what is done in this article (from page 793: "Comparison of the Hedging Performance Across Subperiods.

Basically, I have two time-series each divided into two subperiods. The first time-series is the ln daily change of an unhedged position for time t [z] and t+1 [x]. The second is the ln daily change of a hedged position for time t [y] and t+1 [w].

I want to investigate the change in hedging efficiency between the two periods through this formula:

stat = [1 - ( var(w) / var(x) ) ] - [1 - ( var(y) / var(z) ) ]

What I need is to get 10 000 of stat through a bootstrap, and create BCas.

I have managed to get my 10 000 numbers with the following code:

```
# Unhedged t
BootstrapVAR_UHt <- function(X=UHt){
x.boot <- sample(X, size = length(X), replace = T)
var(x.boot)
}
# Unhedged t+1
BootstrapVAR_UHt1 <- function(X=UHt1){
x.boot <- sample(X, size = length(X), replace = T)
var(x.boot)
}
# Hedged t
BootstrapVAR_Ht <- function(X=Ht){
x.boot <- sample(X, size = length(X), replace = T)
var(x.boot)
}
# Hedged t+1
BootstrapVAR_Ht1 <- function(X=Ht1){
x.boot <- sample(X, size = length(X), replace = T)
var(x.boot)
}
N <- 5000
stat <- numeric(N)
for (i in 1:N){
Boot.UHt = BootstrapVAR_UHt (X=UHt)
Boot.UHt1 = BootstrapVAR_UHt1 (X=UHt1)
Boot.Ht = BootstrapVAR_Ht (X=Ht)
Boot.Ht1 = BootstrapVAR_Ht1 (X=Ht1)
stat[i] = (1 - (Boot.Ht1 / Boot.UHt1) ) - (1 - (Boot.Ht / Boot.UHt) )
}
```

However, in order to use boot.ci() to get my confidence interval, I guess I first have to achieve the above through the boot() function, as this saves the results in the correct way for boot.ci().

I have almost broken both my computer and my mind trying to achieve the above through boot(), but I fall short.

This is my first time using both this site and R, so please comment if I'm missing any vital information in the post (and be kind).

Thanks for reading!

EDIT:

The variables are defined like this, where "STATA" is my spreadsheet and "S4B_58" is unhedged and "HS4B_58" is hedged. Rows 2412:2833 is *t* while rows 2836:3066 is *t+1*.

```
UHt <- as.numeric(na.omit(STATA[["S4B_58"]] [c(2412:2833)] ))
UHt1 <- as.numeric(na.omit(STATA[["S4B_58"]] [c(2836:3066)] ))
Ht <- as.numeric(na.omit(STATA[["HS4B_58"]] [c(2412:2833)] ))
Ht1 <- as.numeric(na.omit(STATA[["HS4B_58"]] [c(2836:3066)] ))
```