close Warning: Can't synchronize with repository "(default)" (/var/svn/mms does not appear to be a Subversion repository.). Look in the Trac log for more information.

Changes between Version 1 and Version 2 of Combinations


Ignore:
Timestamp:
Aug 25, 2010, 8:13:14 AM (14 years ago)
Author:
Pedro Gea
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Combinations

    v1 v2  
    1313
    1414Comúnmente la solución de la combinación se plantea como aquélla que satisfaciendo
    15 el sistema de ecuaciones es más cercana (según la distacia de Mahalanobis) a
     15el sistema de ecuaciones es más cercana (según la [http://es.wikipedia.org/wiki/Distancia_de_Mahalanobis distacia de Mahalanobis]) a
    1616las medias a priori de las variables aleatorias.
    1717
     
    4444wA == wB
    4545}}}
     46
     47=== Solución determinista ===
     48
     49El método utilizado para encontrar la solución determinista
     50(no como variables aleatorias) del sistema es reducir la
     51[http://es.wikipedia.org/wiki/Distancia_de_Mahalanobis distacia de Mahalanobis]
     52de la solución (zA, zB) al punto de partida (muA, muB):
     53
     54{{{
     55Dist(Mu, Sigma2, Z) = Sqrt((Mu-Z)'·Sigma2**-1·(Mu-Z))
     56}}}
     57donde Mu y Z son dos vectores columna y Sigma2 una matriz de covarianzas:
     58{{{
     59Mu' = (muA, muB)
     60Z' = (zA, zB)
     61Sigma2 = ((sigmaA**2, 0), (0, sigmaB**2))
     62}}}
     63
     64La solución de minimizar la distancia:
     65{{{
     66Dist(Mu, Sigma2, Z)
     67}}}
     68sujeta al sistema de ecuaciones:
     69{{{
     70B·Z == C
     71}}}
     72viene dad por:
     73{{{
     74Z = Sigma2·B' · (B·Sigma2·B')**-1 · (C-B·Mu) + Mu
     75}}}
     76
     77En nuestro ejemplo, las matrices B y C son:
     78{{{
     79B = (1, -1)
     80C = (0)
     81}}}
     82de modo que la solución encontrada es:
     83{{{
     84zA = muA + sigmaA**2 * (muB-muA) / (sigmaA**2+sigmaB**2)
     85zB = muB - sigmaB**2 * (muB-muA) / (sigmaA**2+sigmaB**2)
     86}}}
     87y simplificando:
     88{{{
     89z = zA = zB = (muA * sigmaB**2 + muB * sigmaA**2) / (sigmaA**2 + sigmaB**2)
     90}}}
     91Nótese que este resultado puede verse como una media ponderada con las inversas
     92de las varianzas de los valores muA y muB:
     93{{{
     94z = (muA * sigmaA**-2 + muB * sigmaB**-2) / (sigmaA**-2 + sigmaB**-2)
     95}}}
     96
     97