清空筛选条件
清空筛选条件

Fitdist gives wrong answer

1 次查看（过去 30 天）
Emmeirrt
Emmeirrt 2024-1-5
评论： Star Strider 2024-1-29
Ran in:
fitdist is givin me mu of 6.43686 with this code. Am i doing something wrong. I'm trying to get mu and sigma values to compare this data with log normal distribution
Add=16
Add = 16
%%info given in question 4
River_flood_level=[516 678 802 683 491 582 625 531 451 726 652 556 470 570 774 670 432 604 881 707 735 614 698 403 548 504 758 1000 419];
River_flood_level=River_flood_level+Add;
%%solution of question 4
fitdist(River_flood_level.','Lognormal')
ans =
LognormalDistribution Lognormal distribution mu = 6.43686 [6.35256, 6.52116] sigma = 0.22162 [0.175873, 0.29973]
  1 个评论
Paul
Paul 2024-1-5
Ran in:
Why do you think the result is incorrect? What are you expecting mu to be?
The fitted distribution looks pretty reasonable given the limited amount of data provided.
Add=16;
%%info given in question 4
River_flood_level=[516 678 802 683 491 582 625 531 451 726 652 556 470 570 774 670 432 604 881 707 735 614 698 403 548 504 758 1000 419];
River_flood_level=River_flood_level+Add;
%%solution of question 4
p = fitdist(River_flood_level.','Lognormal');
figure
histogram(River_flood_level,'Normalization','pdf')
hold on
plot(400:1200,pdf(p,400:1200))

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Star Strider
Star Strider 2024-1-5
Ran in:
The μ and σ values are in log units. Calculatee their exponentials —
Add=16
Add = 16
%%info given in question 4
River_flood_level=[516 678 802 683 491 582 625 531 451 726 652 556 470 570 774 670 432 604 881 707 735 614 698 403 548 504 758 1000 419];
River_flood_level=River_flood_level+Add;
%%solution of question 4
df = fitdist(River_flood_level.','Lognormal')
df =
LognormalDistribution Lognormal distribution mu = 6.43686 [6.35256, 6.52116] sigma = 0.22162 [0.175873, 0.29973]
Mean = exp(df.mu)
Mean = 624.4451
Sigma = exp(df.sigma)
Sigma = 1.2481
figure
hf = histfit(River_flood_level.', 10, 'Lognormal');
hold on
Line = findobj(hf,'Type','line');
ymu = interp1(Line.XData, Line.YData, exp(df.mu));
plot(exp(df.mu)*[1 1], [0 ymu], ':r', 'LineWidth',3)
hold off
Ax = gca;
Ax.XTickLabelRotation = 60;
.
  6 个评论
显示 4更早的评论
Emmeirrt
Emmeirrt 2024-1-29
Thank you! ı just needed the info that they are in lognormal. I fixed the rest myself. Sorry for not answering early
Star Strider
Star Strider 2024-1-29
As always, my pleasure!
No worries!
Thank you for the follow-up (and for Accepting my Answer).

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Lognormal Distribution 的更多信息

标签

产品


版本

R2023b

Community Treasure Hunt

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

Start Hunting!

Translated by