У меня есть набор данных, который в основном числовой, но у одного есть три варианта. Рассматриваемая переменная называется 'engine' из набора данных GASTURBINES, и варианты являются традиционными, aeroderiv и advanced. Мне было поручено написать полную модель второго порядка с переменной отклика y = heatrate и независимыми переменными RPM, CPRATIO и ENGINE. Модель выглядит довольно сумасшедшей:
E (у) = B0 + B1X1 + B2X2 + B3X3 + B4X4 + B5 (X1 ^ 2) + B6 (X2 ^ 2) + B7X1X2 + B8X1X3 + B9X1X4 + B10X2X3 + B11X2X4 + B11X2X4+ B12 (X1 ^ 2) X3 + B13 (X1 ^ 2) X4 + B14 (X2 ^ 2) X3 + B15 (X2 ^ 2) X4
Где X1 = об / мин, X2 = CPRATIO, X3 = 1если двигатель традиционный (иначе X3 = 0) и X4 = 1, если двигатель является аэродинамическим (иначе X4 = 0).
Я использовал это для создания фиктивных переменных:
X3=ifelse(GASTURBINE$ENGINE == "Traditional", 1, 0)
X4=ifelse(GASTURBINE$ENGINE == "Aeroderiv", 1, 0)
Но когда я запускаю регрессию, мне не хватает бета-оценок для следующего: X3, X4, об / мин: X4, CPRATIO: X4 и CPR2: X4.
Вот регрессия, которую я использую:
p5.34=lm(HEATRATE~RPM+CPRATIO+ENGINE +RPM2 + CPR2 + RPM*CPRATIO + RPM*X3+RPM*X4+CPRATIO*X3 +CPRATIO*X4 +RPM2*X3 + CPR2*X4 +CPR2*X3+CPR2*X4, data=GASTURBINE)
summary(p5.34)
Неправильно ли я установил фиктивные переменные? Почему я не получаю бета-оценки для этих терминов?