Statistics Toolbox™

copulapdf - Copula probability density function

Syntax

Y = copulapdf('Gaussian',U,rho)
Y = copulapdf('t',U,rho,NU)
Y = copulapdf(family,U,alpha)

Description

Y = copulapdf('Gaussian',U,rho) returns the probability density of the Gaussian copula with linear correlation parameters rho, evaluated at the points in U. U is an n-by-p matrix of values in [0,1], representing n points in the p-dimensional unit hypercube. rho is a p-by-p correlation matrix. If U is an n-by-2 matrix, rho may also be a scalar correlation coefficient.

Y = copulapdf('t',U,rho,NU) returns the probability density of the t copula with linear correlation parameters rho and degrees of freedom parameter NU, evaluated at the points in U. U is an n-by-p matrix of values in [0,1]. rho is a p-by-p correlation matrix. If U is an n-by-2 matrix, rho may also be a scalar correlation coefficient.

Y = copulapdf(family,U,alpha) returns the probability density of the bivariate Archimedean copula determined by family, with scalar parameter alpha, evaluated at the points in U. family is 'Clayton', 'Frank', or 'Gumbel'. U is an n-by-2 matrix of values in [0,1].

Example

u = linspace(0,1,10);
[U1,U2] = meshgrid(u,u);
F = copulapdf('Clayton',[U1(:) U2(:)],1);
surf(U1,U2,reshape(F,10,10))
xlabel('u1')
ylabel('u2')

See Also

copulacdf, copulaparam, copularnd, copulastat

copulaparam

copulastat