Commit 2af95767 authored by Holger Brandl's avatar Holger Brandl

changed left to inner join when adding gene length for fpkm/tpm normalization

parent e698d45e
......@@ -55,7 +55,10 @@ geneLength = transmute(geneInfo, external_gene_name, ensembl_gene_id, gene_lengt
cLong = gather(counts, replicate, num_alignments, - ensembl_gene_id)
cLong %<>% left_join(geneLength)
## why left join? genes without annotation can not normalized due to missing length. Resulting tables with have missing data
# cLong %<>% left_join(geneLength)
cLong %<>% inner_join(geneLength)
cLong %>%
distinct_all(ensembl_gene_id) %>%
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment