![]() # print (tprb(totalfragments,creads,targetfragments), ", ", end='') # targetfragments = totalfragments*targetfraction # print (targetfraction, " ", totalfragments, " ", reads, " ", targetfragments, " ", tprb(totalfragments,reads,targetfragments), " ", tprb(totalfragments_norRNA,reads,targetfragments),predictCT(targetfragments,qpcrampregion,fragmentsize,genomesize)) #print ("target fraction, total fragments, reads, target fragments, probability of >1 read, probability of >1 read (no rRNA), predicted CT") # print (totalfragments, " ", creads, " ", targetfragments, " ", tprb(totalfragments,creads,targetfragments), " ", tprb(totalfragments_norRNA,creads,targetfragments),predictCT(targetfragments,qpcrampregion,fragmentsize,genomesize)) Here are samples of PNG output from some of the demos. Gnuplot is distributed with a large set of demonstration scripts. #for creads in range(reads,100000000,1000000): Demo scripts for gnuplot version 5 Demo scripts for gnuplot version 5 (pngcairo terminal) See also the demo output for the SVGand canvasterminals. #print ("total fragments, reads, target fragments, probability of >1 read, probability of >1 read (no rRNA), predicted CT") Totalfragments_norRNA = ((totalfragments-targetfragments)*(1-rrnafraction)) + targetfragments Targetfragments = totalfragments*targetfraction Totalfragments = (basesperng*totalng)/fragmentsize Targetfragments = fragments*(ontargetfragments/totalfragments)Ĭt = ctoffset+log10(targetfragments)*ctslope Ontargetfragments = fragmentsize-ampregion from math import combĭef predictCT(fragments,ampregion,fragmentsize,genomesize): Here’s the modeling code from the substack post.
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |