Using dgelist function in EdgeR to calculate RPKM
I’m new to RNAseq and I want to calculate rpkm values from my raw counts but I am unsure if my understanding of how to do so is correct.
From what I understand, I should do the following:
df <- DGEList(counts=counts,genes=my_genes_vector,groups=my_groups_vector)) #counts should be the dataframe with all the counts, genes should be a vector with genes and groups should be a vector to which group a sample belongs (treatment vs no treatment etc)
df <- calcNormFactors(df)
RPKM <- rpkm(df, gene.length=my_gene_length_vector) #my_gene_length_vector should be a vector with all gene lengths
Is this correct?
• 24 views
Read more here: Source link