表达量矩阵分组很复杂也可以使用limma的3大策略
时间:2022-07-26
本文章向大家介绍表达量矩阵分组很复杂也可以使用limma的3大策略,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
论一个很复杂的分组到底该怎么办
第一次看到这么多分组头都大了。首先要考虑如何分组得到grouplist,其次考虑如何在limma包中分组分析。
听说limma包的官方文档中对这些特殊的情况描述的很细致,于是我找到了这张图,觉得和我目前所面临的情况十分相似
首先下载数据,
rm(list = ls()) ## 魔幻操作,一键清空~
options(stringsAsFactors = F)
library(AnnoProbe)
library(GEOquery)
gse=geoChina('GSE51401')
eSet=gse[[1]]
probes_expr <- exprs(eSet);dim(probes_expr)
head(probes_expr[,1:4])
boxplot(probes_expr,las=2)
#probes_expr=log2(probes_expr+1)
phenoDat <- pData(eSet)
head(phenoDat[,1:4])
由于我不会交叉着分组...所以直接把网页上的分组信息复制粘贴存为了TXT格式的GSE51401文件,然后使用R语言读取
a = read.table(file ='GSE51401')
# 分组
index1=grep('TEC',a$V2)
b=a[index1,]
b$Subject = 'TEC'
index2=grep('NEC',a$V2)
d=a[index2,]
d$Subject = 'NEC'
e = a[c(36,38,41,48,49,54,57,62),]
e$Subject = 'TC'
index4=grep('NTC',a$V2)
f=a[index4,]
f$Subject = 'NTC'
h = rbind(b,d,e,f)
g = h$Subject
table(g)
于是g变成了这样的
并且很好的各种都排列在一起
获取平台信息
eSet@annotation
# Illumina HumanHT-12 V4.0 expression beadchip
GPL=eSet@annotation
###这里是GPL570
##对应的找注释平台和包的网页在http://www.bio-info-trainee.com/1399.html
if(F){
if(!require(hgu133plus2.db))BiocManager::install("hgu133plus2.db")
library(hgu133plus2.db)
ls("package:hgu133plus2.db")
ids <- toTable(hgu133plus2SYMBOL)
head(ids)
}else if(F){
getGEO(GPL)
ids = data.table::fread(paste0(GPL,".soft"),header = T,skip = "ID",data.table = F)
ids = ids[c("ID","Gene Symbol"),]
colnames(ids) = c("probe_id")
}else if(T){
ids = idmap(GPL,type = "bioc")
}
probes_anno=ids
head(probes_anno)
genes_expr <- filterEM(probes_expr,probes_anno )
head(genes_expr)
开始分析
n = genes_expr
# 首先是每个组都和第一个组比较,比如时间或者浓度梯度的药物处理
library(limma)
design=model.matrix(~factor(g))
design
fit=lmFit(n,design)
fit=eBayes(fit)
DEG=topTable(fit,coef=2,n=Inf)
head(DEG)
deg_heatmap(DEG, n, g, topn = 20)
DEG=topTable(fit,coef=3,n=Inf)
head(DEG)
deg_heatmap(DEG, n, g, topn = 20)
DEG=topTable(fit,coef=4,n=Inf)
head(DEG)
deg_heatmap(DEG, n, g, topn = 20)
我查了下coef的含义,昨天之前我还真没注意这个小东西的作用。之前都是直接处理两个分组,或者从多个分组中选取两个分组分析,昨天处理的数据全是乱七八糟的分组..刚开始直接就做了,3个分组的limma分析也直接用2个分组的套路分析,然后后来某一刻顿悟...发现哦这样不行,于是想起来之前小洁老师讲过一篇GSE474的例子,于是可以完美的借用呀。
## 然后是每个组都和其它所有样本比较,比如不同的组织
g
str(g)
table(g)
comp1toOther=do.call(cbind,lapply(c('NEC', 'NTC','TC', 'TEC'),function(x) as.numeric( g %in% x)))
library(limma)
comp1toOther
fit=lmFit(n,comp1toOther)
fit=eBayes(fit)
DEG=topTable(fit,coef=1,n=Inf)
head(DEG)
deg_heatmap(DEG, n, g, topn = 20)
DEG=topTable(fit,coef=3,n=Inf)
head(DEG)
deg_heatmap(DEG, n, g, topn = 20)
代码都是Jimmy老师给的,呃以前就是走流程来分析,现在好歹也多了解了一点点细节,哦统计计算方面的问题不予评价。
## 最后是任意选取两个或者多个分组,组合比较
kp=g %in% c('NEC', 'NTC','TC', 'TEC')
genes_expr=n[,kp]
group_list=g[kp]
library(limma)
design=model.matrix(~factor(group_list))
design
fit=lmFit(genes_expr,design)
fit=eBayes(fit)
###可以调节分组
DEG=topTable(fit,coef=2,n=Inf)
head(DEG)
deg_heatmap(DEG, n, g, topn = 20)
火山图
need_deg=data.frame(symbols=rownames(DEG), logFC=DEG$logFC, p=DEG$P.Value)
deg_volcano(need_deg,1)
deg_volcano(need_deg,2)
deg_volcano(need_deg, style = 1, p_thred = 0.05, logFC_thred = 1)
- Selenium2+python自动化25-js处理日历控件
- 转负二进制(个人模版)
- Selenium2+python自动化24-js处理富文本
- 【干货】对抗自编码器PyTorch手把手实战系列——PyTorch实现对抗自编码器
- Selenium2+python自动化23-富文本(自动发帖)
- 2-Sat+输出可行解(个人模版)
- 协同过滤原理及Python实现
- 每周学点大数据 | No.25二叉搜索树回顾(二)
- RBF神经网络及Python实现(附源码)
- 【干货】计算机视觉实战系列03——用Python做图像处理
- Adaboost从原理到实现(Python)
- Selenium+python自动化22-发送各种类型附件邮件
- Selenium2+python自动化38-显示等待(WebDriverWait)
- 逆元(个人模版)
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法