convert raw counts to TPM
Hi,
I need to convert raw counts to TPMs, but I don't know how to get the gene lengths.
I'm using the following package, can you help me out?
install.packages("DGEobj.utils")
library(DGEobj.utils)
TPMs <- convertCounts(raw_counts, unit = "tpm", geneLength, log = FALSE,
normalize = "none", prior.count = NULL)
• 1,121 views
•
link
updated
by
dthorbur
•
written
by
pinheirofabiano •
•
written
by
pinheirofabiano •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Timeout error using Biomart to get gene lengths
written by AlexStar •
-
Error in bmRequest(request = request, httr_config = httr_config, verbose = verbose) : Gateway Timeo…
written by AlexStar •
-
Convert FPKM to TPM in R
written by AlexStar •
-
Warning - longer object length is not a multiple of shorter object length
written by AlexStar •
-
TPM for scRNAseq
written by little_more •
-
Error in apply(counts, 2, function(x) rpkm(x, lengths)) : dim(X) must have a positive length
written by shail.nair05 •
-
TMM normalized TPM
written by firestar
-
Why counts to tpm gives different values compared to tpm from kallisto?
written by newbie
-
What are some references of works that use a TPM expression threshold for filtering samples/genes?
written by n,n
-
Differential analysis heatmap without hclust
written by newbie
Have you read the documentation for this function? You need a vector of transcript lengths where
geneLengthis - it's unclear here if that's a placeholder using the functions variable names or an object. You can find this information in the same place you found the reference transcriptome, or you can generate it yourself with something likesamtools faidx.Regardless, please take a tour of how to ask a good question. There is plenty of information you could have provided that would help us identify where you're stuck.
Note that, if I'm not mistaken, the name
geneLengthis misleading since you probably want to the transcript length. Possibly use the length of the longest transcript within each gene or maybe something more elaborate than that as long as you exclude introns.Good catch. Edited my comment.