MathGroup Archive 2013

[Date Index] [Thread Index] [Author Index]

Search the Archive

MultinormalDistributions an a RobustMultinomialDistribution

  • To: mathgroup at smc.vnet.net
  • Subject: [mg131303] MultinormalDistributions an a RobustMultinomialDistribution
  • From: J Jesus Rico-Melgoza <ricomelgozajjesus at gmail.com>
  • Date: Tue, 25 Jun 2013 21:15:09 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-outx@smc.vnet.net
  • Delivered-to: mathgroup-newsendx@smc.vnet.net

Dear members of the mathgroup

I have been working with multinormal distributions with dimensions about
24. In many cases, the covariance matrices are of deficient rank and
defining distributions via MultinormalDistribution becomes defficult. To
solve this issue I have define my own function like this

RobustMultinormalDistribution[x_, =CE=BC_, =CE=A3_] := Module[{aux, fac, d},
  d = Length[x];
  aux = 0.5 (x - =CE=BC).PseudoInverse[=CE=A3]. (x - =CE=BC) ;
  fac = (2.0 Pi)^(d/2) PseudoDet[=CE=A3]^(1/2);
  Exp[-aux]/fac
  ]
PseudoDet[a_] := Product[Eigenvalues[a][[i]], {i, 1, MatrixRank[a]}]

It works very well, I think!

However, I would like to use it in all the ways other distributions are
used with in mathematica, such as in RandomVariate and PDF functions.

Could anyone suggest how to  do function definition?
Best regards
Jesus Rico-Melgoza




  • Prev by Date: Re: What is f[1]? Advanced question
  • Next by Date: Re: Rotate Dateticks in DateListPlot
  • Previous by thread: Re: Making (a*\[Theta]11)/(1+a^2) into (a/(1+a^2))*\[Theta]11
  • Next by thread: Rearranging the Terms of a Polynomial