Ошибка в оценке параметра распределения дирихле в r

Когда я использую Dirichlet.mle для следующей матрицы, это приводит к:

Error in while ((conv > convcrit) & (iter < maxit)) { :

missing value where TRUE/FALSE needed

conv кажется NaN.

    A = matrix( 
            c(7.319874130189419E-4, 7.319874130189419E-4, 0.9355592727661133,
 7.319874130189419E-4, 0.04137665405869484, 7.319874130189419E-4, 
7.319874130189419E-4, 0.014280209317803383, 7.319874130189419E-4, 
7.319874130189419E-4, 7.319874130189419E-4, 7.319874130189419E-4, 
7.319874130189419E-4, 7.319874130189419E-4, 7.319874130189419E-4,7.319874130189419E-4, 7.319874130189419E-4, 
0.9355592727661133, 7.319874130189419E-4, 7.319874130189419E-4, 
7.319874130189419E-4, 7.319874130189419E-4, 0.04137665405869484, 
7.319874130189419E-4, 7.319874130189419E-4, 7.319874130189419E-4, 7.319874130189419E-4, 7.319874130189419E-4, 7.319874130189419E-4, 0.014280209317803383),  
            ncol=15,             
            byrow = TRUE)

    > sirt::dirichlet.mle(A)

Как решить проблему?

Двенадцать из пятнадцати столбцов в A постоянны. В вашем наборе данных нет информации об изменениях в этих компонентах, поэтому неудивительно, что расчет MLE не выполняется. Вам нужно больше данных.

user2554330 10.09.2018 11:27

Но разве это не ошибка? Разве мы не можем сказать, что когда все (очень немногие) выборки согласны в значении столбца, тогда это MLE для столбца (в соответствии с определением максимального правдоподобия)? Если да, то почему детали реализации (метод, используемый для оценки) должны мешать нам получить MLE?

Shayan 11.09.2018 07:32

Нет, параметры оценки максимального правдоподобия, а не данные. Вы можете с большой уверенностью сказать, что ваши данные не были получены из распределения Дирихле (это не привело бы к подобным связям), поэтому знать, какие значения параметров дают наибольшую вероятность, бессмысленно. Возможно, этот пакет сможет их найти, но я бы не стал винить авторов, если бы они написали его только для работы с наборами данных, где на самом деле дистрибутив Дирихле.

user2554330 11.09.2018 18:34
1
3
60
0

Другие вопросы по теме