Erstellen einer benutzerdefinierten Ergebnistabelle aus vielen linearen Modellen

Question

Ich arbeite in einem analytischen Labor, in dem viele Daten generiert werden. Wir modellieren Anpassungen und interessieren uns für Modellkoeffizienten (Achsenabschnitt, a, b), r_squared und Reststandardfehler. Die Anzahl von bis zu 70 abhängigen Variablen (Antworten) ist ziemlich groß. Ich möchte aus mylms eine Tabelle erzeugen, in der die Zeilen die Antworten a bis f und die Spalten die erwähnten Modellparameter sind.

Ich war nicht in der Lage, die Lösungen anzupassen ich im Netz zu meinem Fall gefunden:

• By selfmade function
• By only calling coefficients

Ein Teil der Lösung Zusammenfassung ist sicherlich (mylms) wo Die Parameter sind gespeichert. Hier

ist das vereinfachte Beispiel möchte ich eine Lösung finden zu:

#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

Sie dies tun können:

mylms$residuals 

oder den Bereich der Residuen zum Beispiel zu bekommen:

apply(mylms$residuals,2,range) 

oder um die Quantile zum Beispiel zu erhalten:

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