Calculate the conditional distribution using copulas

22 views (last 30 days)
Hi all,
As I have determined the best fitted copula to a data set, I need to determine the conditional distributions of these copulas:
For this determination I want to calculate the PDF of the distributions of X and Y=y and determine het cumulative integral of the PDF's. The distribution I should use for the conditional distribution is 'Gumbel'.
gaus_cdf = copulacdf('Gaussian',[x(:),y(:)],0.6436);
clay_cdf = copulacdf('Clayton',[x(:),y(:)],a_clay);
gumb_cdf = copulacdf('Gumbel',[x(:),y(:)],d_gumb);
t_cdf = copulacdf('t',[x(:),y(:)],0.7055,nu);
frank_cdf = copulacdf('Frank',[x(:),y(:)],a_frank);
gaus_pdf = copulapdf('Gaussian',[x(:),y(:)],0.6436);
clay_pdf = copulapdf('Clayton',[x(:),y(:)],a_clay);
gumb_pdf = copulapdf('Gumbel',[x(:),y(:)],d_gumb);
t_pdf = copulapdf('t',[x(:),y(:)],0.7055,nu);
frank_pdf = copulapdf('Frank',[x(:),y(:)],a_frank);
surf(x,y,reshape(gumb_pdf,50,50))
%surf(x,y,ecdf)
xlabel('Temperature')
ylabel('Sunshine duration')
title('CDF Gumbel Copula')
Using the code above, the PDF's and the CDF's can be displayed in 3D, but I need the copulas depictered in 2D. Can anyone help me how to manage this correctly? The plot should like a bit like this:
In this plot, the conditional distribution of P(X=x | Y=y_50) (blue line) and P(X=x|Y=y_90) (pink line) should be plotted.
Thank you!

Answers (0)

Products


Release

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by