FDR计算
FDR的计算很简单,我折腾了一上午主要是因为遇到了以下几个问题。
问题:
- FDR是什么?有什么用?怎么计算?
- 我把几个模型的P值都合并成一个表了,所以每次运算FDR时,我需挑选特定的对象。
- 我有多个模型,所以我想着要如何构建循环。
- FDR的数据类型和输出结果都是vector,所以要如何把vector的结果写回dataframe里
解决:
逐一摸索,再次巩固R的基础知识(PS:雷子要赶明天的组会汇报,所以我必须自食其力)。
- 我的FDR理论学习笔记(微博)
- 挑选特定的数值。因为FDR是多次重复测量中犯错的概率(见第一条FDR的概念),所以测量的次数就会影响校正后的结果(废话),所以一定要把每个model中的p单独挑出来运算!
又因为FDR的运算对象是向量,所以我只提取模型中的p值这一列,而不是整个数据框截出来。
p_value <-