# Create an example dataset
n <- 10
x <- cbind(runif(n), rbinom(n, 1, 0.5))
model <- NormalGLM$new()
params <- list(beta = c(2, 3), sd = 1)
y <- model$sample_yx(x, params = params)
data <- dplyr::tibble(x = x, y = y)
# Fit the model to the data
model$fit(data, params_init = params, inplace = TRUE)
# Resample from the model given data
resample_param(data, model)
Run the code above in your browser using DataLab