R: Wyodrębnij pełne przypadki/uwzględnione obserwacje z modelu liniowego lub zmiennych formuły

Question

Po uruchomieniu m1 <- lm(f1, data=DT) chcę zapisać obserwacje, które są zawarte (w rodzaju "obs < - complete.cases (m1)", ale coś, co działa), więc że mogę uruchomić drugą regresję na tych samych obserwacjach: m2 <- lm(f2, data=DT[obs]).

Alternatywnie, chciałbym uzyskać obserwacje, które są kompletne dla danego zbioru zmiennych zdefiniowanego przez obiekt formuły. Rozważmy ten pseudokod podobny do R:

f1 <- as.formula("y ~ x1 + x2 + x3") 
f2 <- as.formula("y ~ x1 + x2") 
obs <- complete.cases(DT[,list(all.vars(f1)]) 
m2 <- lm(f2, data=DT[obs]) 

Jak to zrobić? W pierwszym przypadku, lm już wykonuje pracę pośrednio; jak mogę to wydobyć? W drugim, all.vars zwraca wektor znaków; w jaki sposób poprawnie utworzyć niecytowaną listę, którą zrozumie DT (data.table)?

Answer 1

Od data.table v1.9.5, na.omit ma argument cols.

na.omit(DT, cols = all.vars(f)) 

Answer 2

Zakładając, że na.action twojego lm() jest domyślnym na.omit, dlaczego nie zadzwonić na.omit na oryginalne dane?

Sprawdź plik pomocy dla na.omit w poszukiwaniu alternatyw.