Описание тега fitdistrplus

Вопросы с тегом
R package: Help to Fit of a Parametric Distribution to Non-Censored or Censored Data
0 ответов

R fitdistrplus: проблема в построении графика - аргумент 'by' слишком мал

Я хотел бы понять проблему с построением объекта fitdistrplus: ftnb = tryCatch(expr = fitdist(unlist(x[,ftr]),distr="nbinom",discrete = TRUE), error = function(err) {}) if (class(ftnb) %in% c("fitdist")) { print(summary(ftnb)) plot(ftnb)) } Резюме м…
03 фев '16 в 00:20
0 ответов

Проблема при использовании fitdist внутри команды tryCatch в R

Я хотел бы понять поведение tryCatch в моем примере кода: require(fitdistrplus) a = as.data.frame(table(rpois(n = 1000,lambda = 6))) fitdist(data = a$Freq,distr = "nbinom",discrete = TRUE) result = tryCatch(expr = (fitdist(data = a$Freq, distr = "nb…
18 дек '15 в 17:54
2 ответа

Масштабирование для оценки масштаба и параметра формы с помощью функции fitdist (пакет fitdistrplus)

Как указано в заголовке, у меня проблема с масштабированием fitdist функция в R (fitdistrplus пакет). Пожалуйста, посмотрите на следующий код: # Initialize arrays for storing result fit_store_scale <- rep(NA, 3) fit_store_shape <- rep(NA, 3) #…
27 окт '15 в 14:44
0 ответов

R fitdist/optim Код ошибки 1

Кто-нибудь знает, есть ли дополнительная информация, когда fitdist функция не удалась с error code 1, Я думаю, это означает: error Code 1 указывает на то, что максимальный предел итерации был достигнут внутри optimAlghorithmus. Что я могу изменить в…
29 июн '17 в 15:38
0 ответов

Стандартные ошибки оценки и размера "NA" (не NaN) при подборе отрицательного биномиального распределения с использованием fitdist()

Это мой первый пост, поэтому, если есть какие-либо улучшения или дополнительная информация, которую я могу дать, пожалуйста, дайте мне знать. Я использую R 3.4.3 на Mac под управлением 10.11.6. Я работаю с данными подсчета. Меня интересует параметр …
30 мар '18 в 14:47
0 ответов

R - пример распределения смеси с помощью fitdistrplus

Я определил смесь LogNormal и Обобщенного Парето (актуар) следующим образом: dlnormgenpar<-function(x, w_lnorm, meanlog=0, sdlog=1, par_shape1=1, par_shape2=1, par_scale=1) { w_par=1-w_lnorm dln=dlnorm(x,meanlog = meanlog,sdlog = sdlog) dpar=dgen…
26 ноя '15 в 15:29
1 ответ

В легенде отсутствует функция ppcomp в пакете fitdistrplus в R

У меня есть большие данные с более чем 20 миллионами значений. Благодаря конфиденциальности и воспроизводимости кодов, я использую mydata для их замены. set.seed(1234) mydata <- rlnorm(28000000,3.14,1.3) Я хочу найти, какие известные дистрибутивы…
05 июл '16 в 16:04
0 ответов

Нахождение остатков в пакете fitdistrplus в R

У меня есть код, где я использую fitdistrplus, чтобы приспособить распределение Вейбулла к различным наборам данных и вывести исходные данные на CDF. ##Performs Weibull fit using the Maximum GOF Estimations zFIT<-fitdist(x,"weibull",keepdata = TR…
27 окт '16 в 12:57
2 ответа

Подгонка нормального распределения в R

Я использую следующий код, чтобы соответствовать нормальному распределению. Ссылка для набора данных для "b" (слишком велика для прямой публикации): ссылка для б setwd("xxxxxx") library(fitdistrplus) require(MASS) tazur <-read.csv("b", header= TR…
10 окт '16 в 15:45
2 ответа

Как исключить НС? (функция фитдиста)

У меня есть кадр данных 100x2 DFN. Бег fitdist в столбце DFN$Lret выдает сообщение об ошибке "Функция mle не смогла оценить параметры с кодом ошибки 100". Я понял, что причина в том, что последний ряд содержит NA. Отсюда бегаю fitdist исключая NA, и…
07 май '16 в 17:27
1 ответ

Создание графиков фитдиста с помощью ggplot2

Я установил нормальное распределение fitdist функция от fitdistrplus пакет. С помощью denscomp, qqcomp, cdfcomp а также ppcomp мы можем построить histogram against fitted density functions, theoretical quantiles against empirical ones, the empirical…
27 июл '15 в 12:27
0 ответов

R optim() L-BFGS-B требует конечных значений 'fn'

Я использую MLE и метод l-моментов для гамма-распределения, экспоненциального распределения и распределения Вейбулла. Я использовал fitdistr() в гамма-распределении. У меня есть сообщение об ошибке, в котором говорится, что L-BFGS-B нужны конечные з…
29 июн '18 в 14:42
3 ответа

Установка распределения Gumbel с помощью fitdistrplus

Я пытаюсь воспроизвести код из этого ответа, однако у меня проблемы с этим. Я использую дистрибутив gumbel из пакета VGAM а также fitdistrplus, Проблема возникает при выполнении: fit = fitdist(data1, 'gumbel', start = list(location = 0, scale = 1)) …
18 мар '16 в 14:41
1 ответ

Получите степень свободы для стандартизированного распределения T с MLE

Прежде всего, я заранее благодарю всех за то, что прочитали это. Я пытаюсь подогнать стандартизированное распределение T-Student (т. Е. T-Student со стандартным отклонением = 1) на серии данных; то есть: я хочу оценить степени свободы с помощью оцен…
18 май '16 в 20:17
1 ответ

Как я могу получить фотографии по одной, а не показывать их вместе?

Когда я ввел код plot(x.logis), на выходе четыре изображения, показывающие вместе, но я ищу только одну из них. Как мне построить их отдельно? Вот мой код: x<-c(2800,3260,66.5,195,420,840,1380,469,260,50,209.8,370,27,420,157) y<-log10(x) # loa…
11 июн '16 в 05:15
0 ответов

Как построить усеченные распределения (truncdist) с помощью fitdistrplus?

Я пытаюсь построить кривые соответствия для усеченных распределений из пакета fitdistrplus, используя его функцию plot. library(fitdistrplus) library(truncdist) library(truncnorm) dataNum <- c(433.6668, 413.0450, 435.9952, 449.7559, 457.3629, 498…
17 авг '17 в 04:31
1 ответ

Составить несколько объектов FitDist в одном и том же сюжете с разными цветами?

У меня есть список fitdist объекты, которые я сохранил, используя этот кусок кода. norm_dist_res <- list() for(i in 1:10) { x <- 1+(8000*(i-1)) y <- 8000*i print (x) print(y) norm_dist_res[[i]] = norm_dist_res[[i]] <- fitdist(data=as.vec…
19 апр '16 в 11:39
0 ответов

Создать функцию для интегрирования 2 распределений по x

Я хотел бы создать функцию, чтобы интегрировать произведение двух распределений (PDF и CDF) по x. Это делается для того, чтобы выполнить вероятностную оценку риска и определить "ожидаемый общий риск". Мой код у меня до сих пор, который не работает, …
03 авг '18 в 08:38
0 ответов

Не удалось оценить параметры Вейбулла-Модифицированный Обобщенный

Я пытаюсь оценить параметры обобщенного распределения, модифицированного Вейбуллом, но у меня возникает следующая ошибка: "Функция mle не смогла оценить параметры, с кодом ошибки 100" Что мне делать? Модифицированное Вейбуллом обобщенное распределен…
11 окт '18 в 01:43
0 ответов

Пользовательская функция плотности с помощью fitdistrplus

Можно ли использовать fitdistrplus для настройки пользовательских функций плотности? Например: Как оценить параметры из следующей функции плотности?
12 май '16 в 00:12