7 may. 2011

Distancia de Mahalanobis_rev 4

Os he hablado de un video de Youtube  que calcula la distancia de Mahalanobis a partir de esta tabla
mercurio peces
muestraedadlong.pesomg/kg
12831130.068.12
22428143.0127.89
32820136.089.03
43234130.578.28
52215125.0134.08
62637147.5135.31
72419135.0130.48
82822125.086.48
92426127.0129.47
103021139.082.43
112220121.5127.41
123038150.571.21
132417120.0132.06
142620125.090.85
centroid.26.2924.86132.50105.90

Los pasos para el desarrollo de la Distancia de Mahalanobis son sencillos y el cálculo de las operaciones que realiza los tenéis en la sección de Matrices en Youtube, que dicho sea de paso estoy sorprendido de las muchas visitas que tienen, en especial el de la matriz inversa.
Como sabemos podemos calcular la distancia de Mahalanobis a partir de la formula:
D2 = X.S-1.XT

Siendo X, la matriz centrada de la matriz original, XT su transpuesta y S-1 la inversa de la matriz de covarianzas.
Estas operaciones dan lugar a una matriz:

5.57-0.71-1.02-0-1.122-2.286-2.7720.832-0.11-2.391.52.775-2.25652.0319
-0.712.867-0.25-3.10.5992.52831.2702-2.0460.063-0.7880.211.71-1.3626-0.977
-1.02-0.252.683-1.20.222-1.9241.09460.45-2.173.414-1.640.459-0.67480.5231
-0.04-3.09-1.27.76-3.1471.3834-2.3681.7630.9270.5-2.081-0.6581.3821-1.175
-1.120.5990.222-3.12.382-1.2431.4491-0.3420.337-0.5221.718-2.1741.21790.6431
-2.292.528-1.921.38-1.2436.37070.5894-2.3141.544-1.375-1.1251.408-0.3483-3.179
-2.771.271.095-2.41.4490.58942.1366-0.909-0.431.337-0.308-1.0260.569-0.621
0.83-2.050.451.76-0.342-2.314-0.9091.535-0.330.741-0.111-0.9840.70410.9876
-0.110.063-2.170.930.3371.5438-0.433-0.332.022-2.6631.509-1.4951.3833-0.559
-2.39-0.793.4140.5-0.522-1.3751.33750.741-2.665.134-3.1990.313-0.3699-0.172
1.50.21-1.64-2.11.718-1.125-0.308-0.1111.509-3.1993.084-1.5350.88881.1161
2.771.710.459-0.7-2.1741.4076-1.026-0.984-1.490.313-1.5355.469-4.0553-0.211
-2.26-1.36-0.671.381.218-0.3480.5690.7041.383-0.370.889-4.0553.1587-0.233
2.03-0.980.523-1.20.643-3.179-0.6210.988-0.56-0.1721.116-0.211-0.23251.8163


En la diagonal esta el D2 para cada muestra, por lo que haciendo la raiz cuadrada calculamos la distancia de Mahalanobis para cada muestra.
    D2     D
5.56522.3591
2.86721.6933
2.68271.6379
7.76142.7859
2.38241.5435
6.37072.5240
2.13661.4617
1.53511.2390
2.02231.4221
5.13422.2659
3.08371.7561
5.46862.3385
3.15871.7773
1.81631.3477




No hay comentarios:

Publicar un comentario