从许多线性模型创建自定义结果表 [重复]

Question

我在一个生成大量数据的分析实验室工作。我们进行模型拟合，并对模型系数（截距、a、b）、r_squared 和残差标准误差感兴趣。多达 70 个因变量（响应）的数量相当大。我想从 mylms 生成一个表，其中行是 a 到 f 的响应，列是提到的模型参数。

我无法根据我在网上找到的解决方案来适应我的情况：

• By selfmade function
• By only calling coefficients

解决方案的一部分肯定是存储参数的摘要（mylms）。

这是我想找到解决方案的简化示例：

#To generate a dataset
x <- c(1:30)
x2 <- x*x
a <- x2
b <- jitter(a, factor=10)
c <- jitter(b, factor=30)
d <- jitter(c, factor=40)
e <- jitter(d, factor=50)
f <- jitter(e, factor=60)
z <- cbind(a,b,c,d,e,f)

# The matrix y are the responses of x
y <- as.matrix(z)

# Define the weights for the regression
w <- 1/x

# create linear models
mylms <- lm(y~x+x2, weights=w)

Answer 1

你可以这样做：

mylms$residuals

或者获取残差的范围，例如：

apply(mylms$residuals,2,range)

或获取分位数例如：

apply(mylms$residuals,2,quantile,seq(0.1,0.9,by=0.1))