Наложение ограничения на кусочную модель для обеспечения непрерывности в R

Question

У меня есть следующая подогнанная модель без ограничения:

reg <- lm(y ~ indi_x + x + inter) 

, где indi_x = переменная индикатора для x> 14 и inter = переменная взаимодействия для indi_x и x.

Я хочу наложить ограничение, что indi_x + (inter * 14) = 0, чтобы соответствовать двум сегментам в x = 14. Я использовал функцию I() в lm, но не получаю выход я хочу.

Спасибо!

Answer 1

Если я правильно понимаю, у вас есть два уклона, которые объединены в x = 14, и вы хотите вывести отдельные уклоны (и, возможно, общий перехват?)

Это будет сделано:

reg <- lm(y ~ 1 + x + x : I(x > 14))

Обратите внимание, что x * I(x > 14) теперь является изменением на склоне. Таким образом, абсолютный наклон второго отрезка равен slope_2 - slope_1.