## Not run:
# set.seed(1345)
# n <- 120
# w <- rnorm(n)
# x <- arima.sim(n, model = list(ar=0.8));
# y <- arima.sim(n, model = list(ar=0.8));
# z <- lag(x,-1) + rnorm(n,sd=2) # dependence at lag 1
# # UNIVARIATE VERSION
# res1 <- Srho.test.ts(w, lag.max = 5, B = 10, bw='reference', method='integral'
# ,ci.type="perm") # independence
# res2 <- Srho.test.ts(x, lag.max = 5, B = 10, bw='reference', method='integral'
# ,ci.type="perm") # dependence
#
# # BIVARIATE VERSION
# res3 <- Srho.test.ts(x, y, lag.max = 5, B = 10, bw='reference', method='integral'
# ,ci.type="mbb") # independence
# res4 <- Srho.test.ts(x, z, lag.max = 5, B = 10, bw='reference', method='integral'
# ,ci.type="mbb") # dependence
# ## End(Not run)Run the code above in your browser using DataLab