Raw data 计算TPM

tpm3 <- function(counts,len) {

x <- counts/len

return(t(t(x)*1e6/colSums(x)))

}