How to normalize the counts when you are looking at 3'UTR?
I want to know whether I have elevation at 3'UTR between two conditions. What I did was took all the transcripts and extracted the coverage using bedtools and I computed the RPKM for both conditions. Then I used WT over condition to calculate the log2foldchange value and I got this. To get this plot I grouped the values by transcript id and I sliced maximum 10 of exon rank. I assumed the max would the 3'UTR. so the 1 is 3'UTR. But my guide is saying is that I could be the library size. How to I normalize properly?
slice_max(rank_1, n = 10) %>% mutate(seqnum = 1:length(rank))
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
DESeq2 time-series - full model matrix is less than full rank
written by gogeni5529 •
-
How to extract list of top 10 markers from each cluster?
written by Bioinfo •
-
Geom_bar plot with several x-axis variables
written by sking •
-
How to compute Spearman correlation between numeric and non-numeric data variables ?
written by prithvi.mastermind •
-
survival analysis with several conditions
written by Learner
-
gnomAD reports variants with "NA" allele count and allele frequency
written by j.lunger18 •
-
Differential Expression Analysis of Two Groups Time Series Data with Limma
written by hkarakurt
-
specific RNA-SEQ GLM design
written by gero007 •
-
The best test for the two conditional problem.
written by unique379
-
How To Calculate The Rpkm For The Count Tables Of Rna Seq Data
written by narges •
