Hello all, I'm trying to use DEXseq to analyze usage of regions across different sample. My data is single cell RNA-seq, hence I sum all counts per sample and treat them as bulk RNA-seq data. It's due to DEXseq cannot perform on scRNA data due to memory problem. I additionally add 1 pseudo-count to count matrix (after summation per sample) since size factor estimation uses logarithm of product of counts.
After creating successfully DEXSeqDataSet and estimated size factor, I continue with dispersion estimation. However, It prompted the following error message and I don't know the reason for that
dxd.sct <- estimateDispersions(dxd.sct)
Error in `$<-.data.frame`(`*tmp*`, "dispersion", value = NA) :
replacement has 1 row, data has 0
I'm using R 4.1.3 and DEXseq 1.40.0. Do you guys have any idea about it?
Thanks for any help.
1 answer
Log in to answer this question.
-
DESeq2 error
written by sooni •
-
Interprete plot of dispersion estimation by DESeq2
written by tien •
-
Help requiring DEXSeq in R
written by achanda •
-
Estimating power in RNA-seq differential expression
written by ams686 •
-
Normalization of ChIP-Seq count matrix ??
written by Ankit
-
prepDE.py read length choice for paired-end sequencing reads
written by 1234gingko •
-
RNA-Seq Analysis with DEXSeq
written by justinkablan225 •
-
Normalization scaling factors: formula for applying them to raw counts
written by sovrappensiero •
-
Gene level RNA seq expression measure
written by Bnf83
-
Fitting Distribution Of The Gene Expression Levels Of Rna Seq Data
written by k.nirmalraman
