Как избежать проблемы Log-Quartelive Log-Inf в MLE Function из пакета Stat4? -- r поле с участием dplyr поле с участием statistics поле с участием mle пол Связанный проблема

How to avoid the -Inf log-likelihood problem in mle function from stat4 package?


1
vote

проблема

русский

Я хочу максимизировать функцию вероятности логики, чтобы соответствовать некоторым данным, но функция MLE останавливается с этой ошибкой, когда логическая вероятность является -inf.

Ошибка в Optim (Start, f, метод = метод, гессиан = правда, ...): L-BFGS-B необходимо конечные значения «Fn»

Код:

 <код>   ll <- function(shape, scale) {   x <-     c(2,  12,  61,   9,  25, 110,   7,  27,  17,  15,  11,  17,  15,  25,   0)   - sum(dgamma(x, shape, scale, log = TRUE))   }      m <- stats4::mle(ll, start = list(shape = 2.15, scale = 2), method = "L-BFGS-B", lower = c(0.01, 0.01))   

Кто-нибудь знает, как избежать этой проблемы чистым образом?

Английский оригинал

I want to maximize the log-likelihood function to fit some data, but the mle function stop with this error when the log-likelihood is -Inf.

Error in optim(start, f, method = method, hessian = TRUE, ...) : L-BFGS-B needs finite values of 'fn'

The code is:

  ll <- function(shape, scale) {   x <-     c(2,  12,  61,   9,  25, 110,   7,  27,  17,  15,  11,  17,  15,  25,   0)   - sum(dgamma(x, shape, scale, log = TRUE))   }      m <- stats4::mle(ll, start = list(shape = 2.15, scale = 2), method = "L-BFGS-B", lower = c(0.01, 0.01)) 

Does anyone know how to avoid this problem in a clean manner?

</div
           
         
         

Список ответов

Связанный проблема

7  Ось X не соответствует барплоту  ( X axis does not match barplot ) 
Я использовал функцию barplot () для создания сложенного диаграммы из матрицы. Матрица выглядит так: <код> 1 0.989013 0.010987 2 0.999990 0.000010 3 0.999...

1  Есть ли не-Java внедрение регрессионной модели M5P библиотеки RWEKA?  ( Is there a non java implementation of the m5p regression model of the rweka libr ) 
Я ищу функцию, которая создает модель регрессии M5P, как <код> M5P функция <код> RWeka библиотеки (как <код> M5P функция на основе кода Java, который имеет...

0  R Язык: Как работать с динамически размером вектора?  ( R language how to work with dynamically sized vector ) 
Я изучаю программирование R, и пытаясь понять лучший подход к работе с вектором, когда вы не знаете, последний размер, он в конечном итоге будет. Например, в ...

2  Как изменить график разброса GGPlot2  ( How to modify ggplot2 scatter plot ) 
Вот мой набор данных: <код> df1 = data.frame(Count.amp = c(8,8,1,2,2,5,8), Count.amp.1 = c(4,4,2,3,2,5,4)) Я попробовал <код> library(ggplot2) qplot(C...

0  Скрепление финансовых таблиц с веб-страницы с R, Rvest, RCURL  ( Scraping financial tables from web page with r rvest rcurl ) 
Я пробую разбор финансовых таблиц с веб-страницы. Я продолжал. Но я не могу устраивать список, или данные. Карамент <код> library(rvest) link <- "http://www...

6  403 Ошибка при использовании Rvest для входа в веб-сайт для соскабливания  ( 403 error when using rvest to log into website for scraping ) 
Я пытаюсь высказать страницу на веб-сайте, который требует входа в систему и в целом получение ошибки 403. Я изменил код из этих 2 сообщений для моего сайта...

2  Вопросы с установкой пакета Caret of R в Archlinux  ( Issues with installing caret package of r in archlinux ) 
Я пытаюсь установить пакет <код> SCRIPT="""UPDATE IND_AFRO.DRIVER SET Emp_Id = 1000, update_user_id = 'RIBST-4059' WHERE Emp_Id IN (SELECT Emp_Id ...

2  Package Desolve Package Can Parames включают в себя матрицу?  ( Desolve package can parameters include a matrix ) 
Я пытаюсь кодировать модель Seir, которая является возрастной, стратифицированной; То есть в моих дифференциальных уравнениях у меня есть параметр для массово...

1  Замените отсутствующие значения в ячейке со значением из ячейки выше (N-1) с помощью петли  ( Replace missing values in a cell with a value from the cell above n 1 using a ) 
У меня есть файл данных с тысячами строк, у которых есть пробелы, которые я хочу заполнить значением. Мне нужно заменить пустые ячейки со значениями из них вы...

12  Частные члены в R Справочник класса R  ( Private members in r reference class ) 
Можно ли иметь частные поля элементов внутри класса R r. Играя с некоторыми из онлайн-примеров у меня есть: <код> > Account <- setRefClass( "ref_Account"...

0  Как сделать бесконечно рекурсивный список в R: путать [и [[  ( How to make an infinitely recursive list in r confuse and ) 
Редактировать: Этот вопрос глупо, я путающую [и [((спасибо @josilber), но я не могу удалить его. Как можно сделать бесконечно рекурсивный список, l == l [1]...

0  Ошибка при использовании списков  ( Error when using lists ) 
У меня есть файл .rdata со списком внутри. Однако я получаю сообщение об ошибке при попытке отобразить список. <код> ls(data$Test$t0) Error in list2env(lis...

5  Расчет дней в месяц между интервалом двух дат  ( Calculating days per month between interval of two dates ) 
У меня есть набор событий, которые каждый из которых имеет начало и окончание, но они проходят по объему ряда месяцев. Я хотел бы создать таблицу, которая пок...

2  dlyryr "не обещание" ошибка  ( Dplyr not a promise error ) 
У меня есть набор набора панели, для которого я создал отсталые переменные с помощью функции LAG (). Когда я пытаюсь рассчитать дельта для каждого TimePoint, ...

0  Извлечение данных из нижних слоев в растробрике  ( Extracting data from lower layers in a rasterbrick ) 
Итак, я извлекаю данные из растробрика, который я сделал, используя метод из этого вопроса: Как извлечь данные из растробрика? В дополнение к получению да...

Связанный проблема

7  Ось X не соответствует барплоту 
1  Есть ли не-Java внедрение регрессионной модели M5P библиотеки RWEKA? 
0  R Язык: Как работать с динамически размером вектора? 
2  Как изменить график разброса GGPlot2 
0  Скрепление финансовых таблиц с веб-страницы с R, Rvest, RCURL 
6  403 Ошибка при использовании Rvest для входа в веб-сайт для соскабливания 
2  Вопросы с установкой пакета Caret of R в Archlinux 
2  Package Desolve Package Can Parames включают в себя матрицу? 
1  Замените отсутствующие значения в ячейке со значением из ячейки выше (N-1) с помощью петли 
12  Частные члены в R Справочник класса R 
0  Как сделать бесконечно рекурсивный список в R: путать [и [[ 
0  Ошибка при использовании списков 
5  Расчет дней в месяц между интервалом двух дат 
2  dlyryr "не обещание" ошибка 
0  Извлечение данных из нижних слоев в растробрике