作者dfghpt (乱风)
看板Statistics
标题[问题] 用原始数据和平均後数据做t-test的准确度
时间Mon Jul 4 22:39:20 2011
不好意思 想问问大家t-test怎麽做准确度比较高
因为我不是学统计的 但做生物性实验资料都需要用统计来验证
但最近有个奇怪情况有点搞不太懂
我总共有170个样本资料
不过这170组资料是从六次的不等样本数量 在同一个实验流程和参数下做出来的
我是把六次总共170个样本数的资料直接和对照组做t-test
那都可以得到具有差异的p值(<0.05)
但学长又建议我把六次分别平均 剩下六组平均後数据
再和对照组做t-test 差异性就都不见了= = p值都变高很多
学长原本认为我的样本数大 所以error bar很小 因此轻易就能test出差异
但实质上说不定是无差异的
请问一下专业的各位哪一种方法才是比较准确的??
我对平均後差异性消失也有点百思不解....真的是因为数量的关系吗
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.32.224.66
1F：推 aeglos:已知机率的Type I error与未知机率的Type II error， 07/04 23:11
2F：→ aeglos:你喜欢面对哪一个呢？^^a 07/04 23:12
3F：推 aeglos:把原有的大样本变成极小样本是自废武功啊~~样品永远不嫌多 07/04 23:16
4F：→ gsuper:六次的不等样本数量 , 你要怎麽平均? 07/04 23:27
恩恩~就是六个平均数~
5F：→ gsuper:如果是拿六个平均数 V.S. 控制组 , 平均数具有不偏性 07/04 23:28
6F：→ gsuper:那变异数就会变很小,然後另一方面样本数也变小,当然不显着 07/04 23:29
所以还是应该以原始数据一起下去做检定比较好罗?
※ 编辑: dfghpt 来自: 114.32.224.66 (07/05 02:33)
7F：→ gsuper:6组是同一人做的吗? 时间够近吗? 你感觉实验的准度够吗? 07/05 13:45
8F：→ gsuper:如果看来看去6组都没甚麽差 (观察平均数和变异数) 07/05 13:51
9F：→ gsuper:那就像你现在这样 170 V.S. 控制组就好了 07/05 13:52
10F：→ gsuper:若差很多,那可能要考虑meta-analysis的fixed effect model 07/05 13:52
11F：→ kingbee:题外话问一下 提建议的学长 对统计了解多少 07/06 14:29
他学生物的 所以他也不太懂= =
12F：→ andrew43:重点应是你的「统计单位」是什麽. 一个单位一个重复数 07/06 15:45
统计单位是时间 在实验中每个样本有出现反应的时间 重复数通常只有一个妈@@?
13F：→ andrew43:当然有是有先平均再检定的情况, 但也是统计单位内的平均. 07/06 15:46
14F：→ andrew43:至於样本多易达显着, 这可预期, 但显着的意义为何别搞错. 07/06 15:47
15F：→ andrew43:显着与否也不一定能表达差异量, 以上小小提醒. 07/06 15:48
这边有点不太懂@@ 所以应该增加实验次数来更显示差异性?
※ 编辑: dfghpt 来自: 114.32.224.66 (07/06 20:48)
16F：→ andrew43:我猜你误会统计单位的意思. 例如调查学生身高, 则统计单 07/06 22:22
17F：→ andrew43:位为「一位学生」, 而使一位学生的身高为一个重覆数. 07/06 22:23
18F：→ andrew43:至於重覆数多能干嘛, 老生常谈就是提高检定力 (power), 07/06 22:24
19F：→ andrew43:但检定常常只指出是否有差别, 但不一定能指出差别多少. 07/06 22:25
20F：→ andrew43:所以 "有没有差" 和 "差多少" 不要一概而论即可. 07/06 22:25
21F：→ andrew43:意思就是有显着就有显着啦, 但差多少要再看看. 07/06 22:32
22F：→ gsuper:把六次样本,总和样本,控制组 的 1.样本数 2.平均 3.标准差 07/07 00:06
23F：→ gsuper:都列出来看看如何? 07/07 00:06
24F：→ yhliu:问题及资料不大了解. 07/08 19:31
25F：→ yhliu:设 X(i,j)～i.i.d. N(μ,σ^2), j=1,..,Ni, i=1,...,k. 07/08 19:32
26F：→ yhliu:要检定 H0:μ=μ0, 最佳检定策略是用所有 n=ΣNi 个观测值做 07/08 19:34
27F：→ yhliu:t 检定. 今若计算 Y(i) =Σ{X(i,j): j=1,...,Ni}/Ni, 而用 07/08 19:34
28F：→ yhliu:Y(i), i=1,...,k 当观测值做普通 t 检定, 则: 07/08 19:35
29F：→ yhliu:(1) 当 Ni 不等, 如本例, 结果是错的! 因 Y(i) 变异数不等. 07/08 19:36
30F：→ yhliu:(2) 若所有 Ni 均相同,即等样本情形, 用 Y(i) 而不用 X(i,j) 07/08 19:37
31F：→ yhliu: 是没效率的, 因误差自由度由 n-1 骤减为 k-1. 07/08 19:38