В статистике оценка максимального правдоподобия (MLE) - это метод оценки параметров статистической модели.
1 ответ

Как использовать функцию dmvnorm и mapply вместе

set.seed(1) ### i would like to do this dmvnorm(c(.5,.5), mean= c(2,15), matrix(c(3, 0, 0, 9), 2)) dmvnorm(c(.6,.6), mean= c(5,18), matrix(c(6, 0, 0, 15), 2)) ##### BUT using mapply instead... how can that be done? u1 = c(2,15) sigma1 = matrix(c(3, …
2 ответа

Ошибки, выполняющие оценку максимального правдоподобия по трехпараметрическому Weibull cdf

Я работаю с кумулятивным появлением мух с течением времени (взятых с нерегулярными интервалами) в течение многих лет (хотя сначала я просто пытаюсь сделать один год работы). Кумулятивное возникновение следует сигмоидальному шаблону, и я хочу создать…
22 июл '14 в 19:27
1 ответ

Выберите K в PCA python skitlearn

Я пытаюсь использовать пакет skitlearn для PCA . На веб-сайте документации, приведенном здесь http://scikit-learn.org/stable/modules/generated/sklearn.decomposition.PCA.html Говорят, что если n_components == 'mle', то mle используется для определени…
24 мар '16 в 15:31
0 ответов

Почему в mle2 возвращает subset() ошибку, а подмножество работает вручную

Обновление: для минимального примера прокрутите вниз (не воспроизводит ту же ошибку, но разные результаты) перед моим настоящим вопросом предостережение: я новичок в R, поэтому надеюсь, что в этом есть смысл. Я попытался выяснить проблему, используя…
26 ноя '14 в 14:18
3 ответа

GUI интерфейсы для вычислительно тяжелых программ

Я изучаю варианты наилучшего способа достижения следующего. У меня есть вычислительно сложная модель, построенная на R (она использует MLE в своей основе), и я хотел бы предоставить внешний интерфейс GUI для использования этой модели. Ввод в модель …
29 янв '14 в 16:17
2 ответа

Как зациклить все значения в кадре данных в качестве начального значения в maxLik

Я делаю оценку максимального правдоподобия, используя maxLik, что требует указания начальных значений. Вместо указания одного значения, есть ли способ, который позволяет мне использовать все значения из матрицы в качестве начального значения? Мой те…
24 сен '17 в 02:43
1 ответ

Численный метод для модели SARIMAX с использованием R

Мой друг в настоящее время работает над своим заданием по оценке параметра модели временного ряда, SARIMAX (сезонный ARIMA Exogenous), методом максимального правдоподобия (MLE). Используемые им данные относятся к месячному количеству осадков за 2000…
18 сен '16 в 18:11
0 ответов

Что такое скрытые (ненаблюдаемые) данные? а что такое скрытые данные в GMM?

Я изучаю EM алгоритм и GMM вместе На самом деле, я не могу понять это об алгоритме EM в Wiki The EM algorithm is used to find the maximum likelihood parameters of a statistical model in cases where the equations cannot be solved directly. Typically …
04 сен '14 в 13:00
0 ответов

Python ARMA MLE (Реализация алгоритмов из литературы)

обзор Я пытаюсь реализовать оптимизацию параметров авторегрессионного скользящего среднего (ARMA) с использованием оценки максимального правдоподобия (MLE) с помощью фильтра Калмана. Я знаю, что могу соответствовать моделям ARMA, используя пакет st…
01 июн '17 в 07:32
0 ответов

Оценка максимального правдоподобия модели обучения с подкреплением

Я пытаюсь смоделировать скорость обучения для задачи, которая похожа на задачу обучения по усилению бандитов. Единственная разница здесь заключается в том, что вознаграждение не является стохастическим (если сделан правильный выбор, агент награждает…
31 май '18 в 23:19
1 ответ

Оценка максимального правдоподобия в склеарне

В моем проекте я использую GridSearchCV в sklearn исчерпывающий поиск по указанным значениям параметров для модели, чтобы найти наилучшие возможные значения параметров. Я просто протестировал это в RandomForestClassifier и помог мне найти лучшее max…
01 май '17 в 16:13
1 ответ

Запись правильного нормального логарифмического правдоподобия в R

У меня есть проблема относительно следующей модели, где я хочу сделать вывод о µ и tau, u - известный вектор, а x - вектор данных. Логарифмическая вероятность У меня проблема с записью правдоподобия в R. x <- c(3.3569,1.9247,3.6156,1.8446,2.2196,…
0 ответов

Как найти максимум многомерной функции в R

У меня есть совместная вероятность того, что мне нужно максимизировать, что определяется примерно 25 различными переменными. Я надеялся, что найдется метод для нахождения комбинации значений переменных, который максимизирует эту функцию, не прибегая…
20 май '14 в 14:34
1 ответ

Оценка максимального правдоподобия в R нечисловой ошибке аргумента

Это моя функция правдоподобия для отрицательного биномиального распределения. Что я делаю не так, чтобы получить эту ошибку? Ошибка в k + z: нечисловой аргумент бинарного оператора Код R
20 мар '17 в 22:49
0 ответов

Неограниченная ошибка конечных разностей все еще существует даже при использовании mle2

Для предъявления по месту требования, Я выполнил анализ MLE для отдельного набора данных задачи двоичного выбора (каждый набор данных включает в себя несколько испытаний, в которых конкретный субъект должен выбрать вариант A или B) с помощью R. Снач…
11 янв '16 в 18:51
1 ответ

Как выбрать "стартовые" значения в функции stats4::mle в R?

Плохо, что я не могу предоставить воспроизводимый код моей проблемы, потому что я использую большой вектор. Но вопрос все еще довольно ясен. Я пытаюсь приспособить бета-дистрибутив с помощью функции mle: params <- mle(nll, start = list(alpha = 0.…
29 май '17 в 21:02
1 ответ

R optim(){fExtremes} получает 0 гессенскую матрицу

Я использую R {fExtremes}, чтобы найти лучшие параметры распределения GEV для моих данных (вектор). но получите следующее сообщение об ошибке Ошибка в solve.default(fit$hessian): Подпрограмма Лапака dgesv: система точно единственная: U[1,1] = 0 Я пр…
23 авг '16 в 21:34
1 ответ

Оценка максимального правдоподобия MATLAB

Привет я хотел бы сделать MLE оценки моих параметров, используя встроенные функции в Matlab. Вот что говорит Matlab: phat = mle(data,'distribution',dist) Я не знаю, как использовать вектор "данные". Предположим, у меня есть 340 наблюдений, дающих 0,…
03 фев '15 в 15:10
0 ответов

Функция вероятности нескольких данных и функция mle2 из пакета bbmle в R

Я написал пользовательскую функцию правдоподобия, которая подходит для модели с несколькими данными, которая объединяет данные по повторной поимке и телеметрии (sensu Royle et al. 2013 Methods in Ecology and Evolution). Функция правдоподобия разрабо…
27 апр '16 в 01:28
0 ответов

Оценка параметров настраиваемого распределения для MLE с использованием R

Я хочу оценить параметры моего пользовательского дистрибутива, используя mle (stats4). для этого кода R library(stats4) data=rweibull(1000, shape=5, scale = 1.5) rn=data n=length(rn) r=max(rn)+0.05 ll= function(b,d){-(sum(log(r-rn)+n*log(d)-n*log(ga…
19 июл '18 в 07:30