Home > brede > brede_pdf_t.m

brede_pdf_t

PURPOSE ^

brede_pdf_t - PDF of multivariate t-distribution

SYNOPSIS ^

function p = brede_pdf_t(t, dof, my, Sigma);

DESCRIPTION ^

 brede_pdf_t          - PDF of multivariate t-distribution

       function p = brede_pdf_t(t, dof, my, Sigma); 

       Input:  t      t-values
               dof    Degrees of freedom
               my     Mean value (Optional, default=0) 
               Sigma  Covariance (Optional, default=1)

       Output: p      Probability density value

       This function computes the probability density function for
       the multivariate Student's t-distribution.

       For a univariate input the standard t-value (t') can be
       computed as 

                t' = (t - my) / sqrt(Sigma)

       The Gaussian approximation is used for dof>300.

       Ref: Mardia, 1979, p57, exercise 2.6.5

       Example:
         x = (-5:0.1:5)';
         p1 = brede_pdf_t(x, 5, 0, 1);
         p2 = brede_pdf_t(x, 100);
         p3 = brede_pdf_gauss(x);
         p4 = brede_pdf_cauchy(x);
         figure, semilogy(x,p1,'-',x,p2,'--',x,p3,'x',x,p4,'o')
         legend('5 DOF', 'Many DOF', 'Gaussian', 'Cauchy')

       See also BREDE, BREDE_PDF_CAUCHY, BREDE_PDF_GAUSS,
                BREDE_CDF_T.

 $Id: brede_pdf_t.m,v 1.3 2004/10/12 23:14:55 fnielsen Exp $

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:
Generated on Fri 27-Nov-2009 18:11:22 by m2html © 2005