欲哭无泪的p-value = 0.051 | 做几次重复能得到较低的p-value
已有 1469 次阅读 2017-12-15 14:12 |个人分类:RNA-seq|系统分类:科普集锦|关键词:RNA-seq| RNA-seq
本文转载自嘉因微信公众号，已获得授权。查看最新文章，敬请关注嘉因，微信ID：rainbow-genome
作者：小哈 来源：嘉因
60分万岁，多1分浪费
p-value = 0.051。。。。。。
场景一：做RNA-seq，做几次重复？应该做几次？发paper时认可几次重复？
场景二：RNA-seq做了3次重复，用p-value < 0.05筛出的差异基因太少，只用2次重复来筛，筛出了好多，好开森 ~ ~ ~
场景三：基因KO组跟对照组比，计算药物处理后表型是否有显著差异。做了3次实验，p-value > 0.05；继续做，做到5次重复，p-value < 0.05，好开森 ~ ~ ~
减少重复次数？增加重复次数？怎样是合理的？
来看看statQuest视频讲Power in statistics，需要几次重复能获得较低的p-value。
倒霉蛋儿不只你一个，增加重复次数会发生什么？
多做了重复后，30%的情况p-value降到0.05以下，make sense？NO！那叫假阳性
不可以为了获得好的p-value而增加重复，那会增加假阳性。
正确的做法是：
- 做实验之前，算一下到底需要多少样本量；
- 如果实验前没评估样本量，就评估之后重新开始做实验。
（你们这些做统计的不懂实验狗的忧伤）
怎样评估需要做几次重复呢？
实验前评估需要多少次重复，power calculation，受4个因素影响：
第4个是统计方法，最广泛应用的t-test有最强的power。重点考虑前三个因素：
最好是有前期数据或已发表的数据，有第六感也成：
如果没有前期数据和第六感，还有办法：
对于RNA-seq数据：
具体怎么算？
交给电脑算，点击左下角“阅读原文”直达G*Power网站。刚更新，17 July 2017 - Release 3.1.9.3，有Windows版本和Mac版本。用法简单：
你可能需要一些基础知识，搬来简洁易懂的statQuest视频。
one or two tail t-test的另一个视频
回到开头的场景，解决方案：
才不管怎么算呢，We always use 3 patients...
原文地址：https://www.cnblogs.com/wangprince2017/p/9813367.html