Probability Plot with Confidence Interval

2011 9 28
2 个回答
回答已采纳
更新时间：2024 7 5
11 次查看（30 天）

 采纳的回答

UJJWAL
UJJWAL 2011-9-28

1 个投票

Hi Jessica, Functions like normplot and probplot allow you to mark the probability distribution of your data with respect to a specified distribution like Normal or Rayleigh etc.
Confidence Interval is something that they do not evaluate so they do not include an option explicitly for them. You will have to specify the confidence interval values of the parameters of your distribution and then you can plot them. I think the below program will help :-
clc;
clear all;
a= random('normal',10,4,1,100); % Generate a normally distributed random numbers
pd = fitdist(a','normal'); % Creates a ProbDistUnivParam object by fitting the data to a normal distribution. For details look at Documentation
ci = paramci(pd); % This function calculates the values of the parameters based on a certain confidence interval. Here the by default the confidence interval is 95 percent
probplot(a);
h=probplot(gca,@(a,x,y)normcdf(a,x,y),[ci(1,1),ci(1,2)]);
set(h,'color','r','linestyle','-')
t= probplot(gca,@(a,x,y)normcdf(a,x,y),[ci(2,1),ci(2,2)]);
set(t,'color','g','linestyle','-.')
For further help reply back or mail. I hope it will help you.
HAPPY TO HELP
UJJWAL

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Jacob Misura
Jacob Misura 2017-6-12
Just found this question and this answer is amazing thank you for helping people!!!!
Dennis Craggs
Dennis Craggs 2021-1-9
编辑：Dennis Craggs 2021-1-9
This code worked for a normal distribution, but not for a lognormal. So I made changes from normal to lognormal as follows:
pd = fitdist(a,'lognormal');
ci = paramci(pd);
probplot('lognormal',a);
h = probplot(gca,@(a,x,y)logncdf(a,x,y),[ci(1,1),ci(1,2)]);
t = probplot(gca,@(a,x,y)logncdf(a,x,y),[ci(2,1),ci(2,2)]);
The confidence limit lines were not plotted around the best fit line. Any idea what went wrong?
Diego
Diego 2024-7-5
Can you do it for a distribution gev?

请先登录，再进行评论。

更多回答（1 个）

arif rabbani
arif rabbani 2013-5-2

0 个投票

please describe the following term, h=probplot(gca,@(a,x,y)normcdf(a,x,y),[ci(1,1),ci(1,2)]);

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Dennis Craggs
Dennis Craggs 2021-1-9
gca is a reference to the last graphic created.
a is the original data set
x and y are used for the confidence intervals. I don't know the code, but the ci is a variance covariance matrix. h and t are used to reference the upper confidence bound and t for the lower confidence bound.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Half-Normal Distribution 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by