Вопросы по теме 'mle'

Максимальное правдоподобие для многомерного нелинейного нормального распределения
Я пытаюсь написать код для определения g, sigma и lambda с максимальной вероятностью для следующего Уравнения функции . Мне пришлось ограничить интервалы, чтобы гарантировать, что уравнение не принимает логарифм отрицательного числа, однако...
76 просмотров
r mle
schedule 04.10.2021
Использовать преобразование уравнения в функции Matlab?
Я пытаюсь написать очень простую функцию оценки среднего правдоподобия, которая полагается на функцию fminsearch() из Optimization Toolbox. Я пытаюсь передать дескриптор функции в качестве входных данных и оптимизировать сумму neg. журнала...
38 просмотров
function matlab optimization anonymous-function mle
schedule 01.12.2021
Оценка максимального правдоподобия при моделировании с R
Я генерирую 10000 значений из своего дистрибутива и сохраняю в объекте val_1 . Затем я сохранил val_1 в объекте pop. Затем оценил параметр того же распределения в исследовании повторной выборки с номером. итераций 5000 и размера выборки 100. Код...
187 просмотров
r statistics simulation mle
schedule 17.11.2021
Как включить оптимизацию с помощью optim(), возвращающей близкую оценку?
Во-первых, мне нужно уточнить, что я прочитал следующие сообщения, но моя проблема все еще не может быть решена: R optim() L-BFGS-B нуждается конечные значения 'fn' - Вейбулл Оптимизация optim() в R ( L-BFGS-B требует конечных значений...
787 просмотров
r nonlinear-optimization mle
schedule 27.02.2022
Степенной закон, подобранный функцией `fitdistr()` в пакете `fitdistrplus`
Я генерирую некоторые случайные переменные, используя функцию rplcon() в пакете poweRlaw data <- rplcon(1000,10,2) Теперь я хочу знать, какие известные распределения лучше всего соответствуют данным. Логнорм? опыт? гамма? сила закона?...
2494 просмотров
r mle fitdistrplus power-law
schedule 17.03.2022
Почему не выходит beta.fit?
import numpy as np import matplotlib.pyplot as plt from scipy import stats observed = [0.294, 0.2955, 0.235, 0.2536, 0.2423, 0.2844, 0.2099, 0.2355, 0.2946, 0.3388, 0.2202, 0.2523, 0.2209, 0.2707, 0.1885, 0.2414, 0.2846, 0.328, 0.2265, 0.2563,...
929 просмотров
python statistics distribution mle
schedule 28.03.2022
MLE с векторами в качестве аргументов
Я пытаюсь использовать функцию mle() в R, но кажется, что аргумент начального списка может быть только списком отдельных элементов, а не векторов. Есть ли способ передать список векторов? Ниже то, что у меня есть сейчас: inits <-...
120 просмотров
r mle
schedule 04.04.2022
Ошибка оценки параметра распределения Вейбулла
Я использовал следующую функцию для оценки трехпараметрического распределения Вейбулла. library(bbmle) library(FAdist) set.seed(16) xl=rweibull3(50, shape = 1,scale=1, thres = 0) dweib3l <- function(shape, scale, thres) { -sum(dweibull3(xl ,...
161 просмотров
r statistics distribution mle weibull
schedule 28.04.2022
Использование fminsearch для оценки параметров
Я пытаюсь найти оценку максимального правдоподобия журнала для распределения Гаусса, чтобы оценить параметры. Я знаю, что в Matlab есть встроенная функция, которая делает это, подбирая распределение Гаусса, но мне нужно сделать это с помощью logMLE,...
526 просмотров
matlab mle fminsearch
schedule 24.05.2022
`optimize()`: оценка максимального правдоподобия скорости экспоненциального распределения
Я действительно борюсь с пониманием расчетов MLE в R. Если у меня есть случайная выборка размера 6 из распределения exp(λ) , результаты наблюдений: x <- c(1.636, 0.374, 0.534, 3.015, 0.932, 0.179) Я рассчитал MLE следующим образом...
1318 просмотров
r mle
schedule 12.06.2022
Расчет уровней возврата на основе GPD в различных пакетах R
Я выполняю анализ экстремальных значений для метеорологических данных, а точнее для данных об осадках, доступных в мм/сутки. Я использую подход порогового превышения для оценки параметров обобщенного распределения Парето методом максимального...
2782 просмотров
r quantile mle
schedule 12.06.2023
Как определить начальные точки метода максимального правдоподобия
В данный момент занимаюсь монтажом раздатки. Я использовал функцию fitdistr , но у меня возникла проблема с определением начальных точек для MLE. Например, я хочу сопоставить свои данные (матрица осадков-13149 на 1) с гамма-распределением....
1039 просмотров
r mle
schedule 11.03.2023
Ограниченный MLE с Python
Я делаю реализацию MLE с Python. Моя функция логарифмического правдоподобия имеет 5 оцениваемых параметров, и два из них имеют ограничение, что они должны быть между 0 и 1. Я могу реализовать MLE, используя модуль GenericLikelihoodModel из пакета...
1390 просмотров
python statsmodels mle
schedule 12.05.2023
R optim(){fExtremes} получает матрицу Гессе 0
Я использую R {fExtremes}, чтобы найти лучшие параметры распределения GEV для моих данных (вектор). но получите следующее сообщение об ошибке Ошибка вsolve.default(fit$hessian): подпрограмма Лапака dgesv: система точно в единственном числе:...
1484 просмотров
r optimization math hessian mle
schedule 03.03.2023
Weibull MLE в R
Логарифмическая функция правдоподобия выглядит следующим образом . Я написал эту функцию в R, возвращающую отрицательное значение, умноженное на логарифмическую вероятность. Я сгенерировал вектор с помощью rWeibull с параметрами shape=1.5,...
3585 просмотров
r mle weibull
schedule 18.07.2023
подгонка обратного гауссова распределения к данным в R
Я пытаюсь использовать функцию fitdist в R, чтобы подогнать данные к трем различным распределениям с максимальной вероятностью для их сравнения. Lognormal и Weibull работают нормально, но я борюсь с Inverse Gaussian. Мне нужно указать начальные...
1686 просмотров
r gaussian mle
schedule 14.03.2023
Поставка mle с градиентами в R
Моя цель - предоставить градиенты моей целевой функции журнала (плотности) в mle(), поскольку это, по-видимому, ускоряет процесс сходимости и делает его более стабильным (согласно...
76 просмотров
r gradient estimation derivative mle
schedule 29.11.2022
построение логарифмической функции правдоподобия в R
Привет, я получил ошибку длины x и y в моем коде. datasim — это смоделированная выборка размером 1000. Пожалуйста, помогите мне с этим. x <- datasim loglik <- function(theta){ k<- theta[1] lambda<- theta[2] out <-...
1617 просмотров
r statistics log-likelihood mle
schedule 08.01.2023
Использование mle2 для оценки параметров с ошибкой и предсказанием
Я использую mle2 для оценки параметров нелинейной модели, и мне нужны оценки ошибок вокруг оценок параметров (стандартная ошибка). Кроме того, я хотел бы использовать модель для прогнозирования с помощью новых данных, и у меня возникают проблемы...
178 просмотров
r predict hessian non-linear mle
schedule 02.04.2024