Ho usato lm()
per adattarsi a più modelli di regressione, per più variabili di risposta (~ 1 milione) in R. Es.Ottieni gli errori standard residui di un oggetto "mlm" restituito da `lm()`
allModels <- lm(t(responseVariablesMatrix ~ modelMatrix)
Questo restituisce un oggetto di classe "mlm", che è simile a un oggetto enorme contenente tutti i modelli. Voglio ottenere la somma residua dei quadrati per ogni modello, che posso fare utilizzando:
summaries <- summary(allModels)
rss1s <- sapply(summaries, function(a) return(a$sigma))
Il mio problema è che penso che la funzione "Riepilogo" calcola un sacco di altre cose, troppo, ed è quindi piuttosto lento. Mi chiedo se c'è un modo più veloce di estrarre solo la somma residua dei quadrati per il modello?
Grazie!
Questo doesn rispondi davvero alla domanda Stai utilizzando la funzione di riepilogo che il richiedente ha detto è piuttosto lenta. Chiedono se c'è un modo più rapido per estrarre solo la somma residua dei quadrati per il modello? ' – Ren
La mia colpa. Grazie – eFinance