Почему мои остатки для линейной модели в R имеют на одно меньшее значение, чем набор данных, который я использовал для его создания?

Question

Я создал линейную модель в R, используя данные округа из ACS. В моем наборе данных содержится 3140 записей, и все они имеют свои соответствующие коды fips. Я пытаюсь составить карту остатков из моей линейной модели, но у меня есть только 3139 остатков. Кто-нибудь знает, что делает R при создании линейной модели, которая отвечает за это, и как я могу это исправить, чтобы я мог создать эту карту? Спасибо!

В ответ на предложение проверить АН, я запустил это:

which(completedata$fipscode == NA)
integer(0)

R код, если это поможет:

sectorcodes <- read.csv("sectorcodes1.csv") #ruralubrancode, median hh income
sectorcodesdf <- data.frame(sectorcodes)
religion <- read.csv("Religion2.csv")
religiondf <- data.frame(religion)

merge1 <- merge(sectorcodesdf,religiondf, by = c('fipscode'))
merge1df <- data.frame(merge1)

family <- read.csv("censusdataavgfamsize.csv") #avgfamilysize
familydf <- data.frame(family)

merge2 <- merge(merge1df, familydf, by = c('fipscode'))
merge2df <- data.frame(merge2)

gradrate <- read.csv("censusdatahsgrad.csv")
gradratedf <- data.frame(gradrate)
evenmoredata2 <- merge(gradrate,merge2df, by=c("fipscode"))

#write.csv(evenmoredata2, file = "completedataset.csv")
completedata <- read.csv("completedataset.csv")
completedatadf <- data.frame(completedata)

lm8 <- lm(completedatadf$hsgrad ~ completedatadf$averagefamilysize*completedatadf$Rural_urban_continuum_code_2013*completedatadf$TOTADH*completedatadf$Median_Household_Income_2016)
summary(lm8)

library(blscrapeR)
library(RgoogleMaps)
library(choroplethr)
library(acs)
attach(acs)
require(choroplethr)

dataframe1 <- data.frame(completedatadf$fipscode,completedatadf$averagefamilysize)
names(dataframe1) <- c("region","value")

dataframe2 <- data.frame(completedata$fipscode,completedata$hsgrad)
names(dataframe2) <- c("region","value")

residdf <- data.frame(lm8$residuals)
dataframe3 <- data.frame(completedata$fipscode,lm8$residuals)
names(dataframe3) <- c("region","value")
county_choropleth(dataframe1)
county_choropleth(dataframe2)
county_choropleth(dataframe3)

Когда я пытаюсь запустить dataframe3, появляется сообщение об ошибке:

dataframe3 <- data.frame(completedata$fipscode,lm8$residuals)
Error in data.frame(completedata$fipscode, lm8$residuals) : 
arguments imply differing number of rows: 3140, 3139

Answer 1

Это может быть вызвано NA в ответе. Например, используя встроенный фрейм данных BOD, обратите внимание, что в этом примере есть 5 невязок, а в b: 6 строк:

b <- BOD
b[3, 2] <- NA
nrow(b)
## [1] 6
fm <- lm(demand ~ Time, b)
resid(fm)
##          1          2          4          5          6 
## -0.3578947 -0.2657895  1.6184211 -0.6894737 -0.3052632 

Мы можем справиться с этим, указав na.action = na.exclude при запуске lm. Обратите внимание, что теперь есть 6 невязок, из которых дополнительным является NA.

fm <- lm(demand ~ Time, b, na.action = na.exclude)
resid(fm)
##          1          2          3          4          5          6 
## -0.3578947 -0.2657895         NA  1.6184211 -0.6894737 -0.3052632 

Answer 2

Попробуйте

data.frame(na.omit(completedata$fipscode), lm8$residuals) 

Возможно, ваши данные имеют NA значения.