# NOT RUN {
# data
Y <- subset(ifit, id == 1)[,-1]
# fit model with alias (var_estimate also works)
fit <- var_estimate(Y, progress = FALSE)
plts <- plot(summary(fit))
plts$pcor_plt
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab