model1 = AR1(phi = .99, sigma = 1)
model2 = WN(sigma2=1)
col1 = gen.gts(model1, N = 1000)
col2 = gen.gts(model2, N = 1000)
testMat = cbind(col1, col2, col1+col2)
testLts = lts(testMat, unit = 'sec', process = c('AR1', 'WN', 'AR1+WN'))
plot(testLts)
Run the code above in your browser using DataLab