摘要 |
A method of estimating transcript abundances includes: (a) obtaining transcript fragment sequencing data from a potential mixture of transcripts of a genetic locus of interest; (b) assigning this data to genetic coordinates of the locus of interest to obtain a data set of fragment genetic coordinate coverage and a coverage envelope curve; (c) setting a number of transcripts of the mixture; (d) pre-setting a probability distribution function of modelled genetic coverage for each transcript i composed of the product of a weight factor αi and the sum of at least 2 probability subfunctions j independently weighted by a weight factor βi,j; (e) adding the probability distribution functions for each transcript to obtain a sum function; (f) fitting the sum function to the coverage envelope curve to optimize αi and βi,j to increase the fit; and (g) repeating steps (e) and (f) until a pre-set convergence criterion has been fulfilled. |