因此,我需要从线性模型中获取R平方索引,我试图在循环中应用汇总函数,但只显示了第一四分位数,第三四分位数,均值,中位数等。
以下是代码,感谢您的事先帮助。
function_model <- function(n, n_sim, x_mean, x_sd, eps_mean, eps_sd, alpha, beta, ...) {
matrix_results <- matrix(ncol=2, nrow=n_sim)
for (i in 1:n_sim) {
x <- runif(n, min = 0, max = 6)
eps <- rnorm(n, mean=0, sd = x^2)
y <- alpha + beta * x + eps
matrix_results[i,] <- lm(y ~ x)$coefficients
}
z <- lm(y ~ x)
par(mfrow=c(1, 2))
hist(matrix_results[, 1], main="Histogram of Alpha", xlab="Alpha",...)
hist(matrix_results[, 2], main="Histogram of Beta", xlab="Beta",...)
}