Question: DESeq2: before differential analysis
gravatar for limin201709
12 months 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 • 358 views
ADD COMMENTlink modified 12 months ago by Devon Ryan88k • written 12 months ago by limin20170910
gravatar for Devon Ryan
12 months ago by
Devon Ryan88k
Freiburg, Germany
Devon Ryan88k 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 12 months ago by Devon Ryan88k

sorry, it is to avoid zero counts.

ADD REPLYlink written 12 months 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: 795 users visited in the last hour