查看感兴趣基因的表达量及其相关关系

问题来源：
探寻：NCR2 表达对GBM肿瘤芯片数据中cell-cycle 和 cytokines 通路基因的相关关系

通过KEGG.db 来下载KEGG相关通路的基因

if (!requireNamespace("BiocManager", quietly = TRUE))install.packages("BiocManager")
BiocManager::install("KEGG.db", version = "3.8")library(KEGG.db)

去官网查看细胞周期相关的通路：
http://www.genome.jp/kegg-bin/show_pathway?hsa04110
细胞因子相关的通路：
http://www.genome.jp/kegg-bin/show_pathway?hsa04060
用KEGG来下载相应的基因列表：

cellcycle_genes=KEGGPATHID2EXTID[['hsa04110']]
cytokine_genes=KEGGPATHID2EXTID[['hsa04060']]

下载这些基因在GBM芯片数据中的表达量

library(cgdsr)
#library(DT)
mycgds = CGDS("http://www.cbioportal.org/")
all_TCGA_studies <- getCancerStudies(mycgds)  #cbioportal 中相关的研究
all_TCGA_studies[grepl('gbm',all_TCGA_studies$cancer_study_id),]  #找出与gbm研究相关的研究
all_tables <- getCaseLists(mycgds, "gbm_tcga")
all_dataset <- getGeneticProfiles(mycgds, "gbm_tcga")#最后确定要
my_dataset <- 'gbm_tcga_mrna_U133'
my_table <- "gbm_tcga_mrna_U133"
#下载那两个基因list表达图谱
cellcycle_expr <- getProfileData(mycgds, cellcycle_genes, my_dataset, my_table)
cytokine_expr <- getProfileData(mycgds, cytokine_genes, my_dataset, my_table)
#下载单个基因的表达向量
NCR2_expr <- getProfileData(mycgds, 'NCR2', my_dataset, my_table)

关于getProfileData 函数的说明：
mycgds : A CGDS objects
cellcycle_genes: 基因名字的一个向量
mydataset： genetic profile ID 的向量或者单个 genetic profile 字符串
前三个参数都是必须参数
my_table : A case list ID

画出散点图

#定义一个画散点图的函数
cor_plot <- function(x,y){plot(x,y, xlab = 'NCR2', ylab = 'gene')model = lm(y ~ x)summary(model) int =  model$coefficient["(Intercept)"]slope =model$coefficient["x"]abline(int, slope,lty=1, lwd=2, col="red")   r= format(cor(x,y),digits = 4)p= format(cor.test(x,y)$p.value,digits = 4)title(main = paste0('p value=',p), sub  = paste0('correlation=',r))
}#执行函数画出散点图
cor_plot(x=NCR2_expr$NCR2,y=cellcycle_expr$CCNB1)
cor_plot(x=NCR2_expr$NCR2,y=cellcycle_expr$BUB1B)
cor_plot(x=NCR2_expr$NCR2,y=cytokine_expr$CCL4)

查看感兴趣基因的表达量及其相关关系相关推荐

R语言---查看指定癌症中感兴趣基因的表达量---笔记整理
原文链接:https://mp.weixin.qq.com/s?__biz=MzAxMDkxODM1Ng==&mid=2247486845&idx=1&sn=b735a4690 ...
差异表达基因热图怎么看_学徒作业我想看为什么这几个基因的表达量相关性非常高...
不知道多少人看完了全网第一个单细胞课程(免费基础课程),免费学习地址在B站:https://www.bilibili.com/video/av38741055 而且完成结业考核20题:https:/ ...
衡量基因相对表达量的RPKM、FPKM、TPM详解
衡量基因相对表达量的RPKM和FPKM.及TPM 1.RPKM(Reads Per Kilobase per Million)和FPKM(Fragments Per Kilobase per Mill ...
荧光定量PCR：基因相对表达量计算方法
荧光定量PCR之后计算目的基因的相对表达量一般采用2-△△ct的方法.我们还是假设对照组和处理组各有三个生物学重复(即对照组3个cDNA样品cDNA1, cDNA2, cDNA3,处理组3个cDNA样 ...
2021.01.05丨根据基因名称拼接表达量与相关注释
这一步是在进行最后的数据汇总工作中用到的,将基因的count与FPKM值和基因注释的结果组合在一起,得到一个完整的数据.方便客户进行后续研究.算法与之前那篇基因ID匹配注释文本一文相似,用了两个for ...
TCGA任意基因在配对样本中的表达量展示
前几天看到了生信技能树的推文:什么情况下我们可以修改基因名字里面提到了2个函数很好用: 其实这个需求我知道在小洁老师的R包tinyarray里有函数可以实现:宝藏R包tinyarray:常用图表一键 ...
基因去重表达矩阵有重复的基因名探针名字
每次都要检测数据 >dat[1:4,1:4] library(hgu133plus2.db) ids=toTable(hgu133plus2SYMBOL) #toTable这个函数:通过看hgu ...
基因共表达网络分析java,RNA-seq数据的基因共表达网络分析
目录 1 背景 BackgroundTypes of biological networks Motivation for using co-expression networks Network i ...
写给对 ”游戏开发” 感兴趣的朋友们
我们程序员能为世界带来什么? 有些程序员做出了淘宝/Amazon,为世界带来了电子商务. 有些程序员做出了微信/Twitter,为世界带来了社交平台. 而有些程序员做出的东西,虽然无法应用于我们的 ...

查看感兴趣基因的表达量及其相关关系

查看感兴趣基因的表达量及其相关关系相关推荐

最新文章

热门文章