Is there a way to get edgeR to take a matrix of count data containing say 10 subjects and 10 genes then perform regression for each gene ONLY using subjects with non-zero counts for that gene? So some genes may have 10 subjects, some 9 some 8 etc. This is for the purpose of investigating how 0 counts are impacting my results.
Barring that does anyone know how I would program a simple R loop to use lm() to do the same thing? Im merely intereted in comparing the genes' t-statistics with 0s in, and with 0s out.
•
written
by
RNAseqer
1 answer
No, you cannot subset to different columns for different genes using edgeR.
However you can do that with limma. Just convert the counts to logCPMs, then set the values for zero counts to NA.
Log in to answer this question.
-
identifying 3'utrs for multiple transcripts of genes from a summarized experiment.
written by RNAseqer
-
LASSO regression for gene selection and clinical variables
written by fabian •
-
Filtering genes in scRNA-seq
written by learner-MD •
-
filter out gene with zero counts in differential gene analysis
written by tommy •
-
Calculating power using ssizeRNA
written by anamaria
-
Graphing timecourse data from DESeq2
written by swbarnes2
-
RNA Seq Analysis:Feature Counts are zero
written by Genomics •
-
Correlation analysis for low expressed genes
written by kmsh410 •
-
RNA-seq edgeR normalization
written by GuriGuri1565 •
-
Finding driver genes from a cancer dataset
written by Gene_MMP8
