作者gsuper (绿色苏打心)
看板Statistics
标题[问题] 一个运算方法看不懂
时间Sun Sep 6 19:45:11 2009
请问一下
The Mahalanobis distance is
simply quadratic multiplication
of mean difference & inverse of pooled covariance matrix.
诚心请教 simply quadratic mutiplication 是怎麽算?
-------------------------------------------------------
范例如下
-------------------------------------------------------
inverse of pooled covariance matrix | mean difference
|
X Y |
|
X 0.332 -0.047 | -2.3
|
Y -0.047 0.526 | -0.9
|
-------------------------------------------------------
经过 simply quadratic multiplication 运算後 得到的结果
1.41
------------------------------------------------------
我实在看不懂中间做了什麽运算
尤其是 2*2 的方阵搭配 2个数字的向量
是如何汇整成一个 1.41 的单一数字
------------------------------------------------------
simply quadratic multiplication
直翻是 : 简单 二次 乘法
辜狗後发现
可能翻成 "简单平方法" #应该不是最小平方法
可是这个 1.41 到底是怎麽算的还是不了解
-----------------------------------------
http://people.revoledu.com/kardi/tutorial/Similarity/MahalanobisDistance.html
连结的最下方可看到上述范例
求高人指点
--
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.239.247
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 19:46)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 19:46)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 19:47)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 19:48)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 19:49)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 19:53)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 19:56)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 20:09)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 20:10)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 20:13)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 20:13)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 20:23)
1F：→ casella123:不就是你的出处第一个公式吗？ 09/06 21:33
2F：→ casella123:(xbari-xbarj)^T=c(-2.3,-0.9) 09/06 21:35
3F：→ casella123:S^{-1} =inverse of pooled covariance matrix 09/06 21:36
可是第一个公式不是
横的向量 * 反矩阵 * 直的向量?
我用 R 去乘 得到得的结果是 "错误 : 非调和阵列"
还是不会算阿~~~sad
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 22:02)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 22:11)
4F：→ bmka:先不要用R, 你手算会算吗? 09/06 22:31
5F：→ gsuper:不会 我就是要求手算的算法阿 09/06 22:32
6F：→ gsuper:我会手算到最後一步 就差这步就能开始想马式距离的意义了 09/06 22:33
7F：→ bmka:那...你学过矩阵吗?没有的话, 查一下书吧 09/06 22:34
8F：→ bmka:这已经不是统计问题了 09/06 22:34
9F：→ gsuper:B大 我卡的点也不是矩阵运算阿 而是共变异矩阵搭配均差 09/06 22:37
10F：→ gsuper:是卡在看不出最後的运算的意义 不知道这是怎麽算的 09/06 22:37
11F：→ gsuper:更不晓得为什麽在算 希望能了解算法後 再回头想统计意义 09/06 22:40
12F：→ bmka:不就只是..把值代进你上面自己写的公式而已... = = 09/06 22:43
13F：→ gsuper:范例中的这6个数字 我凑来凑去也凑不出 1.41 的结果 09/06 22:43
14F：→ bmka:你倒底想问什麽啊 09/06 22:44
15F：→ bmka:那就是你的矩阵运算过程出错了,不然就是最後忘了开根号 09/06 22:45
是不是这样?
t(vector) * S^-1 * vector
m*n矩阵 1*2 * 2*2 * 2*1
会产生一个数字?
因为我用 R 跑不出来
还以为这样算是不对的
所以这样会产生一个数字出来?
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 22:52)
16F：→ bmka:先用手算 09/06 22:53
17F：→ bmka:你真的会矩阵计算?这很明显的出来一定是scaler啊 09/06 22:54
算出来了 洒花!!!!!!
> (t(z) %*% y %*% z)^0.5 z=mean difference
y=共变异反矩阵
我好开心
之前算不出来是因为把 %*% 用成 *
就是因为矩阵运算很烂
之前看这个矩阵运算指令的时候先跳过.............
谢谢 C大教我回去看第一行公式
还有 B大提醒我不要忘记开根号
谢谢
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 23:00)
18F：推 jangwei:这位同学,你标错重点了,你在算的是马氏距离,是距离的其中 09/06 22:58
19F：→ jangwei:一种而已,就如同你学过的汉明距离或欧氏距离 09/06 22:59
20F：→ jangwei:is後面那段文字是在解释如何算马氏距离 09/06 23:00
21F：→ jangwei:simply quadratic multiplication并不是一个专有名词。Orz 09/06 23:23
22F：→ jangwei:也许你该加强你的英文阅读能力和基本的矩阵运算能力。 09/06 23:23
跨领域该加强的地方太多了 (叹~~)
数学和电脑的书对我来讲 还要一段时间才能适应他们的语言
先前光是搞懂一个 "向量"
就让我痛苦了很久
也是这里的大大帮忙我才学会的
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 23:25)
※ 编辑: gsuper 来自: 140.113.239.247 (09/06 23:27)
23F：→ bmka:加油吧..这种努力的态度还是值得鼓励 09/06 23:30