Делаем выбор модели в python/scipy

Во время выбора модели иногда требуется тест отношения правдоподобия или анализ с использованием BIC (байесовского информационного критерия). Хотя я определенно мог бы сделать это вручную, мне было интересно, есть ли какие-либо scipy-функции, предназначенные для этого?

Я задаю этот вопрос, потому что считаю, что должен существовать способ проведения такого анализа или, по крайней мере, функция для получения значения вероятности.

PS: я не думаю о подгонке одного распределения, вместо этого я думаю о том, чтобы посмотреть на некоторые одномерные данные, которые меняются со временем (т.е. предсказание модели меняется со временем).

Любая помощь будет оценена по достоинству!



Пример для этого вопроса:

У меня есть некоторые данные, которые выглядят так. Фиктивные данные

И вот, у меня есть две модели — одна с четырьмя параметрами, другая вложенная в нее модель с двумя параметрами (фиксация двух других).

Я хочу выполнить тест BIC / отношения правдоподобия, чтобы увидеть, будут ли два свободных параметра иметь существенное значение.


person Yuxiang Wang    schedule 06.04.2014    source источник


Ответы (1)


В статистических моделях вы можете выполнять тесты отношения правдоподобия и Вальда. Различные информационные критерии также доступны для всех моделей. Есть несколько других методов выбора модели, но мне нужно немного больше узнать о том, что вы делаете, чтобы дать конкретные ответы. Тем временем наша документация должна помочь http://statsmodels.sourceforge.net/devel/

person jseabold    schedule 06.04.2014
comment
Большое спасибо за ваш ответ! Я не обновлял свою тему, поэтому в ней есть пример. Я просмотрел статистические модели, но нашел только тесты с линейной регрессией (в рамках МНК), а не общие модели (например, пример, который я привел в вопросе). Не могли бы вы помочь? Действительно ценю это. Спасибо! - person Yuxiang Wang; 11.04.2014
comment
Можете ли вы написать настоящие модели или предоставить код, который вы используете для их соответствия? До сих пор не ясно (мне), что вы имеете в виду под своими параметрами. Являются ли эти коэффициенты моделью линейной регрессии, какой-то моделью временных рядов или чем-то еще? - person jseabold; 16.04.2014