Проблема в том, что NA не будет вычисляться в функции вычетов. Таким образом, у вас нет одинаковой длины. Свободен для вас о том, как вы хотите справиться с этим. Но если вы удалите все строки с помощью NA, это сработает
survey2=na.omit(survey)
survfit <- lm(Height~Wr.Hnd, data =survey2)
residual <- resid(survfit)
plot(survey2$Wr.Hnd, residual, ylab = "residuals", xlab = "Writing Handspan", main = "Residual plot")