На этот вопрос уже есть ответ здесь:

У меня есть фиктивная переменнаяblack гдеblack==0 белый иblack==1 черный. Я пытаюсь соответствовать линейной моделиlm дляblack==1 только категория, однако выполнение кода ниже дает мне неправильные коэффициенты. Есть ли способ в R запустить модель сif утверждение, похожее на стату?

library(foreign)
df<-read.dta("hw4.dta")
attach(df)
black[black==0]<-NA
model3<-lm(rent~I(income^2)+income+black)