作者eco100 ( )
站内Statistics
标题[问题] 卡方跟K-S的比较
时间Wed Aug 14 12:06:32 2013
假设我有两组资料如下：
d1 d2 d3 d4 d5 d6 d7 d8 d9
控制 21.8 47.7 65 142 245 350.5 478.3 643.1 818.7
实验 11.4 52.1 52.8 173.3 285.1 394.2 521.3 727.3 879.7
我想比较这两组的分配是否有差异
分别用了卡方检定跟KS检定
卡方我把控制组数据当成 期望值Ei 实验组当成Oi
卡方值 = Σ(Oi-Ei)^2/Ei = 46.05 > X^2(0.95,8) p-value = 2.32232E-07
而K-S检定 则是用SPSS的两组独立样本比较
算出来的p-value是 1
我的问题是 这两个方法算出来的p-value 有点天差地远
我是觉得K-S的结果似乎比较合乎数据的样子
是卡方的计算有问题吗？
还请大家赐教，谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.109.32.9
1F：→ andrew43:你做的二个检定不是同一个目的啊. 08/14 17:04
2F：→ andrew43:你要做 Kolmogorov-Smirnov Goodness-of-Fit Test 才是. 08/14 17:19
3F：→ eco100:嗯？我做的不是适合度检定吗？ 就是比较两组数据CDF的差异 08/14 21:39
4F：→ eco100:以最大的差异当做统计量 算出p-value？ 08/14 21:40
5F：→ andrew43:你可以自己检查看看. 08/16 08:33
6F：推 anovachen:这些数值代表的是d_i出现的频率吗? (i=1~9) 08/17 00:04
7F：→ anovachen:更正: 不是频率...= = 08/17 00:06
8F：→ anovachen:不是频率也不是次数...我也不晓得该怎麽做QQ 08/17 00:07
9F：→ eco100:这些值单纯就是实验出来的数据 是连续型的 所以不适合卡方? 08/17 00:15
10F：→ anovachen:相同条件下重复九次实验得到的九笔数据吗？ 08/17 00:46
11F：→ andrew43:我说错了. 你的KS是对的, 但误用 chisq 了. 08/17 01:48
12F：→ andrew43:多谢anovachen问了我没注意的问题. 08/17 01:49
sorry 我自己都没说清楚 我再详细描述一下好了
事实上d1~d9是不同的天数
然後呢 其实每一天 都收集了repeat20次的data
而不同天的repeat 是独立的 (也就是 没有block)
post上来的 是每天的平均值
这种型态的资料 一开始的想法当然是想用 重覆的two-way anova
但是 每天的变异数会不一样 天数愈久 变异数就愈大 所以就不适合
再来从数据上看 随着天数增加 不管实验组或控制组 值都会愈来愈大
但这也是预料之中 真正有兴趣的是想判断出
实验组跟控制组 的表现"pattern"是否一样
或者说 实验组数据的成长幅度 是否跟控制组一样
因为最後就觉得不管怎样 都是直接用平均值做比较
我就没有把始末说这麽清楚 当时主要是想问为何K-S跟卡方的结果差异这麽大
所以目前看来 卡方是不适合的？
那K-S用在这个case上是OK的吗？
或者有其他更适合的方法？
谢谢！
※ 编辑: eco100 来自: 118.166.213.21 (08/17 14:18)
13F：→ yhliu:卡方 与 K-S 都不适用这资料! 08/18 12:04
14F：→ yhliu:计量值反应, 可以考虑 "重复量测" 之分析. 最好用原始数据, 08/18 12:05
15F：→ yhliu:不要直接算平均值再拿平均值分析. 08/18 12:06
16F：→ yhliu:或许也可考虑 ANCOVA? 不过, 误差变异之同质性, 及时间序列 08/18 12:07
17F：→ yhliu:之相关性是问题. 又: 如果认为实验组与控制组成立趋势不同, 08/18 12:08
18F：→ yhliu:则 ANCOVA 不适合. 08/18 12:09
19F：→ eco100:ANCOVA? 糟糕 不太熟 我会再做一下功课 谢谢老师 08/18 12:45
20F：→ eco100:那再请教一下 我有另外做t-test 其实就是每天的repeat 08/18 12:46
21F：→ eco100:分开做t-test 想说其中一天的p-value小於0.001的话 08/18 12:47
22F：→ eco100:这判断 实验组跟控制组pattern是不一样的 这样合理吗？ 08/18 12:48
23F：→ eco100:就 08/18 12:55
24F：→ yhliu:ANCOVA是在两组趋势相同时比较两组之差距. 但 "pattern" 或 08/19 19:51
25F：→ yhliu:许应指两组趋势不同? 至於各日分别比较, 这是以各日结果看两 08/19 19:52
26F：→ yhliu:组是否有差异, 并不是 "pattern" 的差异. 08/19 19:53