Question: DESeq2: before differential analysis
gravatar for limin201709
2.5 years ago by
limin20170910 wrote:

I was doing RNA seq data process, and I would like to normalize the raw count and +0.1 to avoid the zero log2 fold change. But if I +1 to normalize counts data, it becomes a matric and how can I continue the pipeline.

L3Ddata = FG_20141202_TF_20160315_aspFumA1163_RNA_seq[,13:18]
colnames(L3Ddata)=test$samplename# rows - genes, columns - samples
rownames(L3Dcolname)=test$samplename##adapt the name 

L3D_dds = DESeqDataSetFromMatrix(L3Ddata,colData = L3Dcolname, design=~condition)
all(rownames(L3Dcolname)%in%colnames(L3Ddata))##test if the colname of count data equals to rowname of description
L3D_dds=L3D_dds+0.1  #here it becomes a matrix
#normalize counts +0.1 to avid zero counts
L3D_dds<- DESeq(L3D_dds)# takes quite some time
rna-seq R • 637 views
ADD COMMENTlink modified 2.5 years ago by Devon Ryan96k • written 2.5 years ago by limin20170910
gravatar for Devon Ryan
2.5 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

There is absolutely no reason to add a pseudocount to your actual counts. A log2 fold change of 0 occurs when there's no change and has nothing to do with there being 0s in the data. If you're worried about extreme fold-changes due to having zeros, then that's what the lfcShrink() command is for.

ADD COMMENTlink written 2.5 years ago by Devon Ryan96k

sorry, it is to avoid zero counts.

ADD REPLYlink written 2.5 years ago by limin20170910
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 745 users visited in the last hour