2015-01-22 2 views
0

я вписываюсь гамма-распределение на эмпирической функции распределения с помощью $ fitdist $ функции:R: Как получить значения подгонки от распределения?

fit = fitdist(data=empdistr,distr="gamma") 

Затем я использую $ denscomp $ функции для сравнения данных с подобранными значениями:

dc = denscomp(fit) 

Но я хотел бы извлечь из $ fit $ или из $ dc $ фактические установленные значения, т. е. точки гамма-плотности (с установленными параметрами), которые отображаются в $ denscomp $ function.

Есть ли у кого-нибудь представление о том, как я могу это сделать.

Заранее благодарен!

ответ

3

Использование dgamma для прогнозирования плотности для данной квантили:

dgamma(x, coef(fit)[1], coef(fit)[2]) 
+0

Спасибо за ваш ответ. Однако я хотел бы найти способ прямо получить хороший «x» в функции dgamma. – user3771535

+1

Вы, кажется, смущены. Что такое «хороший» x »? 'denscomp' просто определяет' dgamma' для квантилей между 0 и 35. – Roland

+0

Мой вопрос о правильном «x», который нужно поставить, следующий: если я хочу подражать сюжету, создаваемому dencomp, я должен написать hist (empdistr) и plot (dgamma (x, coef (fit) [1], coef (fit) [2])), но какое значение я должен дать переменной x? Я попробовал x = 0:35, но он не работает (это своего рода смещение). – user3771535

Смежные вопросы