set.seed(100)
simulated_data <- reg_simulation1(N = 50,
ni_lambda = 8,
beta = c(0.5,0.5,0.5),
beta_b = 1.5,
dsq = 0.1,
sigmasq = 0.5,
delta = 0.6,
nu = 5.89)
y <- simulated_data$y
X <- simulated_data$X
print(head(y))
print(head(X))
Run the code above in your browser using DataLab