set.seed(12)
library(VGAM)
library(rlang)
library(dplyr)
library(ggplot2)
X <- matrix(data = sample(1:3, 1200 * 5, replace = TRUE), nrow = 1200, ncol = 5)
b <- sample(1:3, 5, replace = TRUE)
e <- rlaplace(1200, 0, 1)
Y <- X %*% b + e
alpha <- 0.05
K <- 10
result <- LLIC(X, Y, alpha, K)
plot_LLIC(X, Y, result)
plots <- plot_LLIC(X, Y, result)
print(plots$Actual_vs_Pred)
print(plots$Coef)
Run the code above in your browser using DataLab