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 (15 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 | |