Entering edit mode
5.4 years ago
kousi31
▴
100
Hi all,
I converted the FPKM read counts obtained from Stringtie into TPM by following this link In R : setting function fpkm -
fpkmToTpm <- function(fpkm)
{
exp(log(fpkm) - log(sum(fpkm)) + log(1e6))
}
setwd ("/DATA/TPM/")
In=read.table("file_t_data_FPKM", stringsAsFactors = F, row.names = 1, header = T)
Data <- data.frame(fpkm=In[,2])
In$TPM <- with(Data, fpkmToTpm(fpkm))
write.table(In, "file_TPM", append = FALSE, sep = "\t", row.names = TRUE, col.names = TRUE)
When I visualized the reads in IGV I found all samples had reads aligned even when the read counts are around zero.
Example loci with 9 replicates (early(E), mid(M) and late(L) lactation stages with 3 biological replicates in each stage).
TID E1 E2 E3 M1 M2 M3 L1 L2 L3
MSTRG.16476.18 0.0034 0.10 2.28 0.0040 0.0046 0.034 18.89 0.64 7.08
In this loci I could visualize read counts in all sample. When the expression in the Mid stage is almost nil.