Установите распределение пуассонов на данные

Я построил гистограмму и хотел бы разместить распределение пуассонов на гистограмме. Для этого я передал вектор координаты гистограммы x и y функции poissfit() для оценки лямбда.Установите распределение пуассонов на данные

К примеру, вот что я сделал:

expecteddata = cat(2,x,y) 
[l,lci] = poissfit(expecteddata)

Мой вывод выглядит так:

l = 

    44.3766 0.0130 


lci = 

    42.8887 0.0003 
    45.8645 0.0724

я предполагаю лямбда я заинтересован в для построения бы be 0.013 (Я думаю, что моя лямбда настолько мала, потому что моя гистограмма является частотной гистограммой). Я сюжет пуассоновскую PDF с помощью:

x = 0:50 
y = poisspdf(x,0.013); 
plot(x,y,'r')

И я получаю в результате сюжета накладывается:

Однако, я думаю, что это встроено распределение выглядит немного странно. Кажется, это не очень «пуассон». Кто-нибудь знает, что я делаю что-то неправильно?

источник

2016-07-07 interstellar

Я думаю, что ваш диапазон x слишком велик (т. Е. Размер шага данных слишком велик, поэтому вы не видите изящество распределения), 0:50, в то время как ваша лямбда находится на уровне 0,013, возможно, попробуйте что-то вроде 0: 0.01 : 10 и посмотреть, лучше ли результат. – GameOfThrows

@GameOfThrows вместо использования 'x = 0: 50', я фактически использовал значения x из моей гистограммы. – interstellar

да, я думаю, что это слишком низкое разрешение (если вы знаете, что я имею в виду), чтобы ваше распределение выглядело правильно, вам действительно нужно более высокое разрешение, чтобы что-то вроде 0: 0.01: 10. – GameOfThrows

«Я построил гистограмму и хотел бы разместить распределение пуассонов на гистограмме».

Насколько я понимаю, вам необходимо подобрать распределение пуассонов к существующей гистограмме измеренных данных. Я считаю, что вы можете использовать функцию fitdist().

Например, если ваши данные x.

[n,bin]=hist(x,100); 
m=n/trapz(bin,n); 
bar(bin,m,'w'); 
hold on 
pd=fitdist(x,'poisson'); 
y=pdf(pd,bin); 
plot(bin,y,'k'); 
hold off;

предоставит вам гистограмму с распределенной кривой пуассона, установленной на нее.

источник

2016-07-07 15:50:14

Когда я пробую этот метод, я получаю красную строку в 0 ... – interstellar

Ohhhhh. Я вижу свою проблему. Я использовал частоту для x в 'fitdist', когда я должен был использовать данные, прежде чем вычислять частоту. – interstellar

Установите распределение пуассонов на данные

ответ

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