diff --git a/bin/gen_tx2gene.r b/bin/gen_tx2gene.r index 1c18ea2..441e423 100755 --- a/bin/gen_tx2gene.r +++ b/bin/gen_tx2gene.r @@ -26,6 +26,9 @@ txdf <- AnnotationDbi::select( tab <- table(txdf$GENEID) txdf$ntx <- tab[match(txdf$GENEID, names(tab))] +txdf$TXNAME <- gsub("\\.*", "", txdf$TXNAME) +txdf$GENEID <- gsub("\\.*", "", txdf$GENEID) + tx2gene <- data.frame( tx = txdf$TXNAME, gene = txdf$GENEID diff --git a/bin/tximport.r b/bin/tximport.r index 1ba0b7c..e5c5cfa 100755 --- a/bin/tximport.r +++ b/bin/tximport.r @@ -20,6 +20,7 @@ txi_gene <- tximport( type = "kallisto", tx2gene = tx2gene, ignoreTxVersion = TRUE, + ignoreAfterBar = TRUE, dropInfReps=TRUE ) @@ -28,6 +29,7 @@ txi_tx <- tximport( type = "kallisto", txOut = TRUE, ignoreTxVersion = TRUE, + ignoreAfterBar = TRUE, dropInfReps=TRUE ) @@ -37,6 +39,7 @@ txi_tx_scaled <- tximport( txOut = TRUE, countsFromAbundance = "scaledTPM", ignoreTxVersion = TRUE, + ignoreAfterBar = TRUE, dropInfReps=TRUE )