# \donttest{
# Simulate data
set.seed(123)
sim <- simTVP(theta = c(0.2, 0, 0), beta_mean = c(1.5, -0.3, 0))
data <- sim$data
# Estimate model
res <- shrinkTVP(y ~ x1 + x2, data = data[1:199, ])
# Calculate LPDS
LPDS(res, data[200,])
# }
Run the code above in your browser using DataLab