Why does lognpdf not integrate to 1?

1 回表示 (過去 30 日間)
Jessica
Jessica 2014 年 10 月 25 日
コメント済み: Jessica 2014 年 10 月 25 日
Hello all,
We know by definition that a pdf must integrate to 1, but I'm finding in my code that this is not the case for a pdf I create using lognpdf and I am wondering why. Can you please help me?
increment = 0.01 ;
x = 0:increment:1 ;
mean_val = 0.25 ;
sigma_val = 0.69 ;
y = lognpdf(x,mean_val,sigma_val) ;
>> trapz(x,y)
ans =
0.358554747437567
You can see that the integral (approximated by the trapezoidal method) is significantly below 1. Why?

採用された回答

the cyclist
the cyclist 2014 年 10 月 25 日
Doesn't it integrate to 1 over the interval [0,infinity], not [0,1]?
  1 件のコメント
Jessica
Jessica 2014 年 10 月 25 日
Ugh. You're right. I had been working with Beta distributions the whole time before this, so my frame of reference was skewed. That's silly.
Now, how to get the lognormal pdf to behave the way I want it to between 0 and 1...anyway, thank you.

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeNumerical Integration and Differentiation についてさらに検索

製品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by