r语言如何绘制蛋白质组和转录组相关性图
发表于:2024-09-23 作者:千家信息网编辑
千家信息网最后更新 2024年09月23日,这篇"r语言如何绘制蛋白质组和转录组相关性图"文章的知识点大部分人都不太理解,所以小编给大家总结了以下内容,内容详细,步骤清晰,具有一定的借鉴价值,希望大家阅读完这篇文章能有所收获,下面我们一起来看看
千家信息网最后更新 2024年09月23日r语言如何绘制蛋白质组和转录组相关性图
这篇"r语言如何绘制蛋白质组和转录组相关性图"文章的知识点大部分人都不太理解,所以小编给大家总结了以下内容,内容详细,步骤清晰,具有一定的借鉴价值,希望大家阅读完这篇文章能有所收获,下面我们一起来看看这篇"r语言如何绘制蛋白质组和转录组相关性图"文章吧。
目前研究蛋白质组的项目越来越多,经常有蛋白质和转录组相关的分析。
deg_file <- 'S-1_S-2_S-3_vs_R-1_R-2_R-3.all'dep_file <- 'S-1_S-2_S-3_vs_R-1_R-2_R-3.DEP.all.xls'# 设置筛选参数dep_fc <- 1.5dep_fdr <- 0.05deg_fc <-1.5deg_fdr <- 0.05# 读取文件deg <- read_delim(deg_file,delim = '\t')dep <- read_delim(dep_file,delim = '\t')dim(dep)# 蛋白ID转换extr_id <- function(x){ return(substr(x, start = 0, stop = 16))}dep$protein_id <- as.character(llply(dep$protein, extr_id))gene_ids <- deg$`#ID`protein_id <- dep$protein_id# 共有IDcommon_id <- intersect(gene_ids,protein_id)# 提取数据列表dep_common = dep[dep$protein_id %in% common_id,c('protein_id','logFC','P.Value')]colnames(dep_common) <- c('protein_id','logFC_p','FDR_p')dep_common$regulate_p = as.factor(ifelse(dep_common$FDR_p < dep_fdr & abs(dep_common$logFC_p) >=log2(dep_fc), ifelse(dep_common$logFC_p>log2(dep_fc),'Up','Down'),'Normal'))table(dep_common$regulate_p)deg_common = deg[deg$`#ID` %in% common_id,c('#ID','FDR','log2FC')]colnames(deg_common) <- c('#ID','FDR_g','log2FC_g')deg_common$regulate_g = as.factor(ifelse(deg_common$FDR_g < deg_fdr & abs(deg_common$log2FC_g) >=log2(deg_fc), ifelse(deg_common$log2FC_g>log2(deg_fc),'Up','Down'),'Normal'))table(deg_common$regulate_g)# 合并转录组和蛋白质结果merge_df <- merge(dep_common, deg_common, by.x='protein_id', by.y='#ID')# 基于调控关系进行分组merge_df$type <- paste(merge_df$regulate_p, merge_df$regulate_g,sep = '_')types <- c('Down_Down','Down_Up','Up_Down','Up_Up')merge_df[!merge_df$type %in% types,'type'] <- 'Other'# 绘图mycol <- c("#e41a1c","#ff7f00","#984ea3","#4daf4a","#80b1d3")ggplot(data=merge_df, aes(x=logFC_p,y=log2FC_g,colour=type)) + geom_point() + scale_color_manual(values = c("Up_Up"=mycol[1],"Up_Down"=mycol[2],"Down_Up"=mycol[3],"Down_Down"=mycol[4],"Other"=mycol[5]))+ geom_hline(aes(yintercept=log2(deg_fc)), linetype="dashed",colour="grey11")+ geom_hline(aes(yintercept=-log2(deg_fc)), linetype="dashed",colour="grey11") + geom_vline(aes(xintercept=-log2(dep_fc)),linetype="dashed",colour="grey11") + geom_vline(aes(xintercept=log2(dep_fc)),linetype="dashed",colour="grey11") + labs(x = "Protein_logFC", y = "Gene_logFC", title = "Protein_Gene_relation") + theme_bw()+ theme( panel.grid=element_blank(), axis.text.x=element_text(colour="black"), axis.text.y=element_text(colour="black"), panel.border=element_rect(colour = "black"), legend.key = element_blank(), legend.title = element_blank())
以上就是关于"r语言如何绘制蛋白质组和转录组相关性图"这篇文章的内容,相信大家都有了一定的了解,希望小编分享的内容对大家有帮助,若想了解更多相关的知识内容,请关注行业资讯频道。
蛋白
蛋白质
内容
相关性
语言
文章
知识
篇文章
价值
参数
大部分
就是
数据
文件
更多
步骤
知识点
结果
行业
资讯
数据库的安全要保护哪些东西
数据库安全各自的含义是什么
生产安全数据库录入
数据库的安全性及管理
数据库安全策略包含哪些
海淀数据库安全审计系统
建立农村房屋安全信息数据库
易用的数据库客户端支持安全管理
连接数据库失败ssl安全错误
数据库的锁怎样保障安全
网络安全法 微信
常熟标准软件开发口碑推荐
网络安全风险评估的形式
执行定点数据库恢复时用什么备份
成都世纪网络技术
沈逸网络安全和信息化
开启数据库的指令
恩度网络技术有限公司
人民银行软件开发中心待遇
怎么上传php数据库
关于网络安全的知识小报
5台服务器
织梦数据库名称密码
商丘讯博网络技术有限公司
execel服务器缺少根元素
本地网站访问远程数据库很慢
服务器磁盘如何处理
关于网络安全意识征文
河西区信息网络技术创造辉煌
用友数据库 重装不了
网络安全 蓝队
云服务器下能安装虚拟机吗
网安系统服务器怎么接线
数据库数据按照时间筛选数据
我的世界服务器生存冒险
数据通信网络网络技术
打开服务器安全组
数据库查询每科成绩的最高分
怎样建立班级成绩管理数据库
数据库约束出错