data(veWaningData)
set.seed(1234)
ind <- sample(1:nrow(x = veWaningData), 2500, FALSE)
# NOTE: This sample size is chosen for example only -- larger data sets
# should be used.
res <- veWaning(data = veWaningData[ind,],
L = 52,
lag = 6,
modelGam1 = ~ X1+X2+A+A:X1+A:X2,
modelGam2 = ~ X1+X2,
modelEntry = ~ X1+X2,
modelPsiGam1 = ~ X1+X2,
modelPsiGam2 = ~ X1+X2,
gFunc = 'piece',
v = c(15,30))
plot(x = res)
Run the code above in your browser using DataLab