umontreal.iro.lecuyer.probdist
Class BivariateTDist

java.lang.Object
  umontreal.iro.lecuyer.probdist.ContinuousDistribution2Dim
      umontreal.iro.lecuyer.probdist.BivariateTDist

public class BivariateTDist
extends ContinuousDistribution2Dim

Extends the class ContinuousDistribution2Dim for the standard bivariate t distribution. The correlation between X and Y is r and the number of degrees of freedom is ν. Its probability density is

f (x, y) = (1 + (x² -2ρxy + y²)/(ν(1 - ρ²)))^-(ν+2)/2/(2π(1-r^2)^1/2),

and the corresponding distribution function (the cdf) is

T_ν(x, y, r) = ∫_-∞^xdx∫_-∞^ydy f (x, y)/(2π(1 - r^2)^1/2).

We also define the upper distribution function called barF as

bar(T)_ν(x, y, r) = ∫^∞_xdx∫^∞_ydy f (x, y)/(2π(1 - r^2)^1/2).

Field Summary

Fields inherited from class umontreal.iro.lecuyer.probdist.ContinuousDistribution2Dim

decPrec

Constructor Summary

BivariateTDist(int nu, double rho)
Constructs a BivariateTDist object with correlation ρ = rho and ν = nu degrees of freedom.

Method Summary

double	barF(double x, double y) .
static double	barF(int nu, double x, double y, double rho) Computes the standard upper bivariate t distribution.
double	cdf(double x, double y) .
static double	cdf(int nu, double x, double y, double rho) Computes the standard bivariate t distribution using the method described in.
double	density(double x, double y) Returns f (x, y), the density of (X, Y) evaluated at (x, y).
static double	density(int nu, double x, double y, double rho) Computes the standard bivariate t density function with correlation ρ = rho and ν = nu degrees of freedom.

Methods inherited from class java.lang.Object

equals, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

BivariateTDist

public BivariateTDist(int nu,
                      double rho)

Constructs a BivariateTDist object with correlation ρ = rho and ν = nu degrees of freedom.

Method Detail

density

public double density(double x,
                      double y)

Description copied from class: ContinuousDistribution2Dim

Returns f (x, y), the density of (X, Y) evaluated at (x, y).

Specified by:

density in class ContinuousDistribution2Dim

Parameters:

x - value x at which the density is evaluated

y - value y at which the density is evaluated

Returns:

density function evaluated at (x, y)

cdf

public double cdf(double x,
                  double y)

Description copied from class: ContinuousDistribution2Dim

. Computes and returns the distribution function F(x, y):

F(x, y) = P[X <= x, Y <= y] = ∫_-∞^xds∫_-∞^ydt f (s, t).

Specified by:

cdf in class ContinuousDistribution2Dim

Parameters:

x - value x at which the distribution function is evaluated

y - value y at which the distribution function is evaluated

Returns:

distribution function evaluated at (x, y)

barF

public double barF(double x,
                   double y)

Description copied from class: ContinuousDistribution2Dim

. Returns the upper cumulative distribution function bar(F)(x, y):

bar(F)(x, y) = P[X >= x, Y >= y] = ∫^∞_xds∫^∞_ydt f (s, t).

The default implementation computes bar(F)(x, y) = 1 + F(x, y) - F(∞, y) - F(x,∞).

Overrides:

barF in class ContinuousDistribution2Dim

Parameters:

x - value x at which the upper distribution is evaluated

y - value y at which the upper distribution is evaluated

Returns:

upper distribution function evaluated at (x, y)

density

public static double density(int nu,
                             double x,
                             double y,
                             double rho)

Computes the standard bivariate t density function with correlation ρ = rho and ν = nu degrees of freedom.

cdf

public static double cdf(int nu,
                         double x,
                         double y,
                         double rho)

Computes the standard bivariate t distribution using the method described in. The code for the cdf was translated directly from the Matlab code written by Alan Genz and available from his web page at http://www.math.wsu.edu/faculty/genz/homepage (the code is copyrighted by Alan Genz and is included in this package with the kind permission of the author). The correlation is ρ = rho and the number of degrees of freedom is ν = nu.

barF

public static double barF(int nu,
                          double x,
                          double y,
                          double rho)

Computes the standard upper bivariate t distribution.