alpha0 <- 0.03
alpha1 <- 0.02
iuse0 <- 4
iuse1 <- 4
phi0 <- -4
phi1 <- 1
tc0 <- c(3,6,9,12)
tc1 <- c(6,12,18,24)
t0 <- c(0.3,0.6,0.9,1)
t1 <- (1:4)/4
rho <- 0
cvecList0 <- gbounds(t=t0,iuse=iuse0,
alpha=alpha0,phi=phi0)
cvec0 <- cvecList0$bd
cvecList1 <- gbounds(t=t1,iuse=iuse1,
alpha=alpha1,phi=phi1)
cvec1 <- cvecList1$bd
result <- sErrRphInt2(cvec0, cvec1,
t0, t1, tc0, tc1, rho)
print(result)
Run the code above in your browser using DataLab