Модель гауссовой смеси - Сингулярная матрица

Question

Модель гауссовой смеси - Сингулярная матрица

1

Во время вычисления модели гауссовой смеси я должен вычислить pdf() многомерного гауссовского распределения. Я делаю это, используя метод scipy multivariate_normal(). Теперь во время вычислений я сталкиваюсь с ситуацией, когда я получаю сообщение об ошибке

LinAlgError: сингулярная матрица

Установка параметра allow_singular scipy multivariate_normal()= True обходит эту ошибку, и я получаю довольно значимый результат:

Итак, как scipy "разрешает" сингулярность и каковы последствия для вычислений? → Если никаких последствий, строго говоря, этот параметр может быть установлен по умолчанию TRUE (но это не так).

2Obe 27 июль 2018, в 13:15

Источник

0

«Во время расчета гауссовой модели смеси ...» Показать код. Без этого все, что может сделать каждый, - это спекулировать. Кому-то будет легче помочь вам, если вы создадите минимальный, полный и проверяемый пример .
Warren Weckesser 27 июль 2018, в 12:53
0

@ Уоррен Векессер, спасибо за ваш ответ. На самом деле вопрос не в том, «как его кодировать», а в том, «что делает Сципи на заднем плане и каковы могут быть последствия для моего результата».
2Obe 27 июль 2018, в 12:56
1

Это правда. Мне любопытно, что происходит в вашем коде, который приводит к числовой единственной матрице, но, как вы сказали, это всего лишь преамбула, приводящая к актуальному вопросу о allow_singular .
Warren Weckesser 27 июль 2018, в 14:21

Показать ещё 1 комментарий

Теги:

python

scipy

linear-algebra

gmm

expectation-maximization

1 ответ

Ещё вопросы

«Во время расчета гауссовой модели смеси ...» Показать код. Без этого все, что может сделать каждый, - это спекулировать. Кому-то будет легче помочь вам, если вы создадите минимальный, полный и проверяемый пример .
@ Уоррен Векессер, спасибо за ваш ответ. На самом деле вопрос не в том, «как его кодировать», а в том, «что делает Сципи на заднем плане и каковы могут быть последствия для моего результата».
Это правда. Мне любопытно, что происходит в вашем коде, который приводит к числовой единственной матрице, но, как вы сказали, это всего лишь преамбула, приводящая к актуальному вопросу о allow_singular .

2Obe · Answer 1 · 2018-07-27T12-39-00.000Z

Этот ответ не отвечает на вопрос о том, как scipy допускает сингулярность, а также не то, что влияет на это поведение на результаты, но дает немного больше информации о том, что происходит во время вычислений:

EDIT: я опубликовал подробный ответ о том, что происходит, когда мы сталкиваемся с проблемами сингулярности при расчетах GMM здесь

Может быть, вы хотите поместить ответ в блокнот Jupyter, чтобы увидеть обозначение LATEX, или посмотрите следующий ответ ( stats.stackexchange.com/questions/359347/gaussian-mixture-model )