Чи коли-небудь є причина не використовувати ортогональні многочлени під час регресії?


13

Взагалі мені цікаво, чи краще коли-небудь не використовувати ортогональні поліноми під час регресії із змінними вищого порядку. Зокрема, мені цікаво використання R:

Якщо poly()з raw = FALSEвиробляє те ж значення, що і підігнаного poly()з raw = TRUE, і polyз raw = FALSEвирішуєте деякі з проблем , пов'язаних з поліноміальними регрессиями, то слід poly()з raw = FALSE завжди можна використовувати для установки поліноміальних регрессий? За яких обставин було б краще не використовувати poly()?

Відповіді:


16

Коли-небудь причина? Впевнений; ймовірно, кілька.

Розглянемо, наприклад, де мене цікавлять значення необроблених коефіцієнтів (скажімо, порівнюйте їх з гіпотезованими значеннями), а колінеарність не є особливою проблемою. Це майже та сама причина, чому я часто не маю на увазі центр у звичайній лінійній регресії (що є лінійним ортогональним многочленом)

Це не речі, з якими ви не можете мати справу через ортогональні поліноми; це більше питання зручності, але зручність - це велика причина, чому я роблю багато речей.

Однак, я прихиляюсь до ортогональних многочленів у багатьох випадках, підходячи до поліномів, оскільки вони мають певні переваги.


чи можна порівняти коефіцієнти, отримані в результаті ортогональної поліноміальної регресії, до гіпотезованих значень?
user2374133

2
Так. Ви можете перетворити їх назад на коефіцієнти, що маються на увазі, і стандартні помилки, наприклад, із "необроблених" поліномів.
Glen_b -Встановіть Моніку

2
Частіше перетворення з ортогональної поліноміальної основи на мономічну основу є погано обумовленим процесом (для високих ступенів; низький ступінь перетворення не надто поганий), тож якщо апріорі цікавляться коефіцієнтами мономічної бази, будь-який чисельну стабільність, яку ви отримали від використання ортогональних многочленів, викидається у вікно при перетворенні, тож ви можете також із самого початку використовувати одночлени. Caveat emptor , звичайно.
JM не є статистиком

1
@JM Дякую, це прекрасний момент. На щастя, в статистичні додатки в наші дні було б дуже рідко розміщувати більше, ніж досить багаточлен поліноміального порядку (моя звичайна порада полягає в тому, що, якщо немає вагомих теоретичних причин, щоб підняти ступінь вищої ступеня три-чотири, слід переглянути різні підходи - яка альтернатива можливо, найкраще залежить від обставин, але такі речі, як сплайни, наприклад, можуть бути придатними для деяких ситуацій.).
Glen_b -Встановіть Моніку

13

Тому що якщо ваша модель залишає R, коли вона виростає, ви повинні пам'ятати, щоб упакувати її константи центрування та нормалізації, і тоді вона повинна весь час їх перетягувати. Уявіть, що один день натрапите на нього з жорстким кодом у SQL, і жах від усвідомлення того, що це їх помилково!

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.