T-test简介:

对分组变量的差异显著性检验,T-test是最为常用的检验方法,t-test要求数据符合正态分布。不符合用(Mann-Whitney U)

正态性检验  shapiro.test(x),  注:p值大于0.05 都不能拒绝原假设,数据是正态的

分析目的:

根据NMR数据找case control 生物标记物

第一步需要T-test进行初筛,数据格式如下:

想要自动进行分组检验,可以用R的which 命令进行数据筛选:

Tdata[2,which(Tdata[1,]=="0")]  #筛选第二行 分组样本为“0”标签的 Glucose的分组数值

最终输出文件:

整体应用代码如下:

Pvalue

log2_FC

for(i in 2:nrow(Tdata)){

if(sd(Tdata[i,which(Tdata[1,]=="0")])==0 && sd(Tdata[i,which(Tdata[1,]=="1")])){

Pvalue[i-1]

log2_FC[i-1]

}

else{

y=t.test(as.numeric(Tdata[i,which(Tdata[1,]=="0")]),as.numeric(Tdata[i,which(Tdata[1,]=="1")]))

Pvalue[i-1]

log2_FC[i-1]

}

}

# 对p value进行FDR校正

fdr=p.adjust(Pvalue, "BH")

# 在原文件后面加入log2FC,p value和FDR,共3列;

out

write.table(out,file="ttest.out.xls",quote=FALSE,sep="\t",row.names=FALSE)

python 相关性fdr矫正_T-test 、 pvalue FDR矫正 R语言实现 付代码实例相关推荐

  1. r语言软件GDINA_finTech MSc代做、代写Python程序语言、代写MSc program、代做Python设计帮做C/C++编程|代写R语言...

    finTech MSc代做.代写Python程序语言.代写MSc program.代做Python设计帮做C/C++编程|代写R语言Strathclyde Business School, finTe ...

  2. python监听文件最后修改人_Python持续监听文件变化代码实例

    在日常的工作中,有时候会有这样的需求,需要一个常驻任务,持续的监听一个目录下文件的变化,对此作出回应. pyinotify就是这样的一个python包,使用方式如下: 一旦src.txt有新的内容,程 ...

  3. python实现一个简单的加法计算器_Python tkinter实现简单加法计算器代码实例

    tkinter 是 Python 的标准 GUI 库.Python 使用 tkinter 可以快速的创建 GUI 应用程序.由于 tkinter 是内置到 python 的安装包中.只要安装好 Pyt ...

  4. python pop check mail_python初学者,用python3实现基本的学生管理系统代码实例

    本篇文章主要分享python学生管理系统的使用,文章非常详细地介绍了通过示例代码实现的学生管理系统,该系统对每个人的研究或工作都有一定的参考学习价值,希望你能在其中有所收获. 这个是用python实现 ...

  5. python批量删除_Python批量删除只保留最近几天table的代码实例

    Python批量删除table,只保留最近几天的table 代码如下: #!/usr/bin/python3 """ 批量删除table,只保留最近几天的table &q ...

  6. python图表并生成html_Python3 pyecharts生成Html文件柱状图及折线图代码实例

    代码如下 pip3 install pyechatrs #! /usr/bin/python3 from pyecharts.charts import Bar, Page from pyechart ...

  7. python爬取微博评论点赞数_Python selenium爬取微博数据代码实例

    爬取某人的微博数据,把某人所有时间段的微博数据都爬下来. 具体思路: 创建driver-–get网页--找到并提取信息-–保存csv--翻页--get网页(开始循环)-----没有"下一页& ...

  8. 用python计算圆柱体的体积和表面积_Java圆柱体表面积和体积计算代码实例

    题目描述 Java编程实现圆柱体表面积和体积的计算,要求从键盘上输入底面圆的半径和圆柱体的高,然后输出结果. 题目分析 首先定义Circle类,然后在main方法中创建Circle类的对象c,利用对象 ...

  9. python编的俄罗斯方块游戏_Python使用pygame模块编写俄罗斯方块游戏的代码实例

    文章先介绍了关于俄罗斯方块游戏的几个术语. 边框--由10*20个空格组成,方块就落在这里面. 盒子--组成方块的其中小方块,是组成方块的基本单元. 方块--从边框顶掉下的东西,游戏者可以翻转和改变位 ...

最新文章

  1. 转载文章,感觉真的很心酸
  2. Ubuntu设置squid代理
  3. code block怎样实现图形界面_微服务入门:Openresty实现API网关
  4. SQL关联查询————LEFT JOIN关键字的使用
  5. 图像入门二之视频操作
  6. “10%时间”:优点和缺点——敏捷海滩会议上Elizabeth Pope的报告
  7. 程序员的那点事(转自java老师李明志)
  8. 冒泡排序(普通,加强对比)
  9. c语言计算器括号怎么解决,C语言计算器,该如何解决
  10. 【批处理bat】注释一段文本/字符串
  11. PXE配置-Tftpd64
  12. 模拟退火算法求解超市布局问题
  13. JUCE-VST-001 JUCE的下载与安装
  14. oracle 拼音首字母查询,用Oracle的NLSSORT获取拼音首字母
  15. GC overhead limit exceeded问题
  16. 骁龙8+参数 骁龙8+什么水平 骁龙8+处理器怎么样
  17. 科幻迷福利-黑客帝国4明年1.14
  18. Appfuse创建网站访问慢的问题
  19. 如何快速深度写论文?
  20. 在html表格中进行计算,如何在word表格中实现计算功能?超级实用,值得细细品味...

热门文章

  1. 梁启超《最苦与最乐》
  2. 信息学奥赛一本通——1012:计算多项式的值
  3. 抖音霸屏软件为什么会成为爆款?
  4. 负数在计算机中的存储——补码
  5. 软件测试的两张藏宝图和三个层次
  6. egret引擎中使用tiled运行在微信小游戏中
  7. 数学分析教程史济怀练习9.1
  8. 公司邮箱如何申请注册?企业邮箱怎么购买登陆?
  9. Your task is to find for a given phone number any of its divisions into groups of two or three digit
  10. 如何有效进行配置三库管理(配置管理)?