r语言如何计算t分布临界值_[统计]从p-value到q-value的计算（附代码）

看到群里有老师在问“如何从p-value计算获得q-value”。今天简单进行下讨论，并附送一份可在命令行下独立运行的代码。首先，关于“P值”我们昨天讨论了“计算个“P值”，咱也不太懂，咱也不敢信”。大家或许还对下面的这张PPT还有印象：对于来自相同分布的两组样本，如果很多次采集数据并进行假设检验，“本不应显著的”两组样本也会出现P值很小的“显著结果”（“P-hacking”）。

​上图中，我们对来自相同正态分布n(5,5)的x和y进行仿真，每次模拟三个重复，并进行t检验；假设我们有这么一个“劳模”进行了1000次（x轴）独立实验，并各自计算出p值（y轴）；其中，红线为p=0.05的阈值线。从图上看出，1000次实验中，有很多次已经出现p<0.05的情况了。实际上，这种情况下“P-value”服从[0,1]区间的“均匀分布”；在随机情况下，1000次实验出现也能出现大约50次“P<0.