搜索
查看: 2364|回复: 1

[mRNA-seq] TCGA数据提取后用DESeq2作分析,结果得到log2foldchange全是1

[复制链接]

4

主题

13

帖子

141

积分

注册会员

Rank: 2

积分
141
发表于 2018-12-29 20:25:17 | 显示全部楼层 |阅读模式
TCGA数据提取后用DESeq2作分析,结果得到log2foldchange全是1,原因可能是什么?
代码如下:
[mw_shl_code=python,true]# 用DESeq2进行差异表达分析
source("https://bioconductor.org/biocLite.R")
biocLite("DESeq")
biocLite("pasilla")
biocLite("BioParallel")
library(limma)
library(DESeq2)
library(pasilla)

foldChange = 1
padj = 0.05
setwd("/Users/tea/Bioinformation/TCGA/work_space/DataTransformation")
rt = read.table("geneName_Matrix_Processed.csv",sep=",",header=T,check.names=F)
rt = as.matrix(rt)
dimnames=list(rownames(exp),colnames(exp))
countData = matrix(as.numeric(as.matrix(exp)),nrow=nrow(exp),dimnames=dimnames)
condition = factor(c(rep("normal",44),rep("tumor",502)))
colData = data.frame(row.names=colnames(countData), condition)
dds = DESeqDataSetFromMatrix(countData, DataFrame(condition), design= ~ condition )
head(dds)  
dds = DESeq(dds, parallel = T)   

library("BiocParallel")
register(MulticoreParam(16))

res = results(dds, parallel = T)
write.csv(res,file= "DESeq2_results_out.csv", row.names = T)[/mw_shl_code]


csv文件截图-2

csv文件截图-2

csv文件截图-1

csv文件截图-1

打开csv文件后,log2foldchange一列全是1...不知道问题可能出在哪里?谢谢大家!



上一篇:clusterProfiler跑enrichkegg总是得不到结果
下一篇:转录组期刊类型
回复

使用道具 举报

2

主题

15

帖子

699

积分

高级会员

Rank: 4

积分
699
发表于 2019-2-28 10:14:26 | 显示全部楼层
[mw_shl_code=python,true]rt = read.table("geneName_Matrix_Processed.csv",sep=",",header=T,check.names=F)
rt = as.matrix(rt)
dimnames=list(rownames(exp),colnames(exp))
countData = matrix(as.numeric(as.matrix(exp)),nrow=nrow(exp),dimnames=dimnames)[/mw_shl_code]
感觉这个程序跑不通,变量exp是哪里来的。
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树 ( 粤ICP备15016384号  

GMT+8, 2020-4-7 02:02 , Processed in 0.027802 second(s), 27 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.