在線性模型R中包含NA

Question

我有一個關於如何在迴歸之後將NA值包含在向量中的問題。

這裏是我的代碼：

m <- c(100, 108, 105, 120) 
n <- c(5, NA, 10, 8) 

b <- lm(n~m)$fitted.values 

# This returns: 

     1  3  4 
6.961538 7.384615 8.653846 

我需要列2被包括在內，無論是與零或NA，所以它由而不是3 4列對如何做到這一點的任何建議？

Answer 1

爲了您的提取擬合值也將包括NA（即4列），可以使用na.exclude內lm()

na.exclude確實casewise缺失相對於預測和準則，在這種情況下，功能fitted()將填補輸出與NA

所以，你的代碼看起來像這樣

m <- c(100, 108, 105, 120) 
n <- c(5, NA, 10, 8) 

# ordinary linear model 
b <- lm(n~m, na.action = na.exclude) 

# extract elements from a linear model fit 
fitted(b) 

# And output of NA column 
fitted(b) 
     1  2  3  4 
6.961538  NA 7.384615 8.653846