转录组分析——差异表达分析

1.数据的类型,数据名称为“GSE124133_gene_count_matrix.txt”
在这里插入图片描述
2.整理分类数据集condition type.txt,数据如下:
在这里插入图片描述
3.利用R语言中DESeq包进行差异表达分析

#查看R语言当前所在的文件夹
getwd()
#进入到数据所在的文件夹所在的路径(这一步很重要,如果不正确,可能导致数据文件无法正常的导入)
setwd("E:/shengwuxinxi/基因组重测序数据操作步骤/作业")
#如果以前没有使用过DESeq,则需要到Bioconductor上查找你的R版本所对应的安装包方式[Bioconductor](http://www.bioconductor.org/)
#加载安装包
library(DESeq)
library(cluster)
library(DESeq2)
#载入"GSE124133_transcript_count_matrix.txt“数据
database_all <- read.table(file = "GSE124133_transcript_count_matrix.txt", sep = "\t", header = T, row.names = 1)
head(database_all)
#载入分组情况文件condition type.txt
colDataMOTSHFD<-read.table("condition type.txt",header=T,sep="\t")
head(colDataMOTSHFD)
##DESeq利用进行差异表达分析分析
ddsMOTSHFD  <-  DESeqDataSetFromMatrix( countData = database_all, colData = colDataMOTSHFD, design =~sex+pdj+status)
dds <- DESeq(ddsMOTSHFD)
res <- results(dds)
resdata<merge(as.data.frame(res),as.data.frame(counts(dds,normalized=TRUE)),by="row.names",sort=FALSE)
head(resdata)
#将计算结果文件,导出表格形式
write.csv(resdata,file = "3-tumorous_v_survive.csv")

4.结果文件为,可以根据实验的具体情况对数据进行筛选,和处理
在这里插入图片描述
说明:本文数据来源于文献《Integrated analysis of lncRNA and mRNA repertoires in Marek’s disease infected spleens identifies genes relevant to resistance》

  • 2
    点赞
  • 25
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值