從多種線性模型創建自定義結果表

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))