# Probability Plot with Confidence Interval

27 visualizzazioni (ultimi 30 giorni)
Jessica Lam il 28 Set 2011
Commentato: Dennis Craggs il 9 Gen 2021
Hi
I would like to plot one probability plot with the confidence interval, like the plot in the below link
But I just found the following Matlab build-in functions , normplot() and probplot() which do not include the confidence interval.
Any suggestion would be greatly appreciated! Many thanks in advance!
Jessica
##### 0 CommentiMostra -1 commenti meno recentiNascondi -1 commenti meno recenti

Accedi per commentare.

### Risposta accettata

UJJWAL il 28 Set 2011
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','-.')
HAPPY TO HELP
UJJWAL
##### 2 CommentiMostra 1 commento meno recenteNascondi 1 commento meno recente
Dennis Craggs il 9 Gen 2021
Modificato: Dennis Craggs il 9 Gen 2021
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?

Accedi per commentare.

### Più risposte (1)

arif rabbani il 2 Mag 2013
please describe the following term, h=probplot(gca,@(a,x,y)normcdf(a,x,y),[ci(1,1),ci(1,2)]);
##### 1 CommentoMostra NessunoNascondi Nessuno
Dennis Craggs il 9 Gen 2021
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.

Accedi per commentare.

### Categorie

Scopri di più su Half-Normal Distribution in Help Center e File Exchange

### Community Treasure Hunt

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

Start Hunting!