# Create sample data for the example
chain_data <- as.data.frame(matrix(rnorm(500 * 4), ncol = 4))
colnames(chain_data) <- c("param1", "param2", "param3", "param4")
# Run the convergence check
conv_results <- check_gaussian_convergence(chain_data)
print(conv_results)
# The plot method for this object can be used to create convergence plots.
# plot(conv_results)
Run the code above in your browser using DataLab