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.
- Timestamp:
-
Aug 25, 2010, 8:13:14 AM (14 years ago)
- Author:
-
Pedro Gea
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v1
|
v2
|
|
13 | 13 | |
14 | 14 | Comú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 |
| 15 | el sistema de ecuaciones es más cercana (según la [http://es.wikipedia.org/wiki/Distancia_de_Mahalanobis distacia de Mahalanobis]) a |
16 | 16 | las medias a priori de las variables aleatorias. |
17 | 17 | |
… |
… |
|
44 | 44 | wA == wB |
45 | 45 | }}} |
| 46 | |
| 47 | === Solución determinista === |
| 48 | |
| 49 | El 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] |
| 52 | de la solución (zA, zB) al punto de partida (muA, muB): |
| 53 | |
| 54 | {{{ |
| 55 | Dist(Mu, Sigma2, Z) = Sqrt((Mu-Z)'·Sigma2**-1·(Mu-Z)) |
| 56 | }}} |
| 57 | donde Mu y Z son dos vectores columna y Sigma2 una matriz de covarianzas: |
| 58 | {{{ |
| 59 | Mu' = (muA, muB) |
| 60 | Z' = (zA, zB) |
| 61 | Sigma2 = ((sigmaA**2, 0), (0, sigmaB**2)) |
| 62 | }}} |
| 63 | |
| 64 | La solución de minimizar la distancia: |
| 65 | {{{ |
| 66 | Dist(Mu, Sigma2, Z) |
| 67 | }}} |
| 68 | sujeta al sistema de ecuaciones: |
| 69 | {{{ |
| 70 | B·Z == C |
| 71 | }}} |
| 72 | viene dad por: |
| 73 | {{{ |
| 74 | Z = Sigma2·B' · (B·Sigma2·B')**-1 · (C-B·Mu) + Mu |
| 75 | }}} |
| 76 | |
| 77 | En nuestro ejemplo, las matrices B y C son: |
| 78 | {{{ |
| 79 | B = (1, -1) |
| 80 | C = (0) |
| 81 | }}} |
| 82 | de modo que la solución encontrada es: |
| 83 | {{{ |
| 84 | zA = muA + sigmaA**2 * (muB-muA) / (sigmaA**2+sigmaB**2) |
| 85 | zB = muB - sigmaB**2 * (muB-muA) / (sigmaA**2+sigmaB**2) |
| 86 | }}} |
| 87 | y simplificando: |
| 88 | {{{ |
| 89 | z = zA = zB = (muA * sigmaB**2 + muB * sigmaA**2) / (sigmaA**2 + sigmaB**2) |
| 90 | }}} |
| 91 | Nótese que este resultado puede verse como una media ponderada con las inversas |
| 92 | de las varianzas de los valores muA y muB: |
| 93 | {{{ |
| 94 | z = (muA * sigmaA**-2 + muB * sigmaB**-2) / (sigmaA**-2 + sigmaB**-2) |
| 95 | }}} |
| 96 | |
| 97 | |