搜索
查看: 1514|回复: 1

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

[复制链接]

4

主题

13

帖子

139

积分

注册会员

Rank: 2

积分
139
发表于 2018-12-29 20:25:17 | 显示全部楼层 |阅读模式
TCGA数据提取后用DESeq2作分析,结果得到log2foldchange全是1,原因可能是什么?
代码如下:
[Python] 纯文本查看 复制代码
# 用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)




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

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x



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

使用道具 举报

2

主题

14

帖子

579

积分

高级会员

Rank: 4

积分
579
发表于 2019-2-28 10:14:26 | 显示全部楼层
[Python] 纯文本查看 复制代码
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)

感觉这个程序跑不通,变量exp是哪里来的。
回复 支持 反对

使用道具 举报

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

本版积分规则

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

GMT+8, 2019-9-22 02:13 , Processed in 0.028323 second(s), 28 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.