我需要从函数迭代R平方

时间:2018-11-15 18:08:03

标签: r

因此,我需要从线性模型中获取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",...)
}

0 个答案:

没有答案