R lm Захватывать условия взаимодействия, но не категориальную переменную - PullRequest
1 голос
/ 24 апреля 2020

Я хотел бы оценить следующую регрессионную модель: y = b0 + b1 * X + b2 * x * dummy

, где y и x непрерывны, а dummy является категориальной (фиктивная переменная).

Другими словами, я хотел бы, чтобы моя оценочная модель оценила три коэффициента: bo, b1 и b2.

Я пробовал следующее ...

lm(y ~ x + x * dummy, data)

, но добавляет в модель переменную dummy и оценивает коэффициент dummy .

к тому, что я хочу сделать, но он преобразует термин взаимодействия в двоичную переменную (true / false).

lm(y ~ x + I(!x * dummy), data)

Для репликации рассмотрим следующий пример:

data <- tibble(y=rnorm(10), x=runif(10), dummy=ifelse(x>.5,1,0))
lm(y ~ x + x * dummy, data)
lm(y ~ x + I(!x * dummy), data)

Спасибо

1 Ответ

1 голос
/ 24 апреля 2020

Здесь:

> summary(lm(y ~ x+ x : dummy, data))

Call:
lm(formula = y ~ x + x:dummy, data = data)

Residuals:
     Min       1Q   Median       3Q      Max 
-0.61312 -0.15558 -0.00354  0.23965  0.47351 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)
(Intercept)  0.06755    0.36162   0.187    0.857
x            0.94953    1.18299   0.803    0.449
x:dummy     -1.10220    0.88112  -1.251    0.251

Residual standard error: 0.4148 on 7 degrees of freedom
Multiple R-squared:  0.2645,    Adjusted R-squared:  0.05438 
F-statistic: 1.259 on 2 and 7 DF,  p-value: 0.3412
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...