Using dgelist function in EdgeR to calculate RPKM

Question

[bart]

Hi,

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?

Thanks!

[Simon Ahn]

I think so. Refer to this website: https://support.bioconductor.org/p/64585/