I am trying to calculating the fpkm values from the htseq-count result. I think I already get the gene.size values for each of the transcript, while the "dds" contains more rows than the gene.size since there are NR##### (non-coding RNAs) in the dds list.
When I was tryin tying to use
mcols(dds)$basepairs <- gene.size there is error codes: Error in `[[<-`(`*tmp*`, name, value = list(gene = 1:33398, length = c(6363L, : 33398 elements in value to replace 33420 elements
I am wondering if anybody can help with this. I am not sure if the dds and gene.size is the ordered in the same way! Many thanks!