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?