Problems? Is your data what you think it is? PerlMonks

### Re: how to make this code more efficient? Need to Compute the Min and Max of a Temperature using an equation

 on Apr 20, 2010 at 18:25 UTC ( #835851=note: print w/replies, xml ) Need Help??

Treat it less as a brute force challenge and think about your data. You have a sting
NNYDBAVDVHVHNGGNR

What parameters will increase the Tm? - C and G residues
What parameters will decrease the Tm? - A and T residues

For each residue, find if IUPAC can give you a C or a G for the max temperature at that residue, then a T or an A for the min temperature.

If you really need to calculate over large numbers of primers then divide the primers up into percentiles, eg 90% GC, 80% GC and so on. Its an approximation but good enough for what I think you are doing.

If you need more accuracy yet, look at your hash values - this shows you the residue pairs that give you most and least values, so see waht pairs of residues are likely for the extremes of these values and bias the calculations to those primers

Just out of curiosity where are you getting the enthalpy and entrophy values from? St. Lucia?

One other thought - if you want accurate temps, you're going to want to calculate on the monovalent and divalent ion concentrations. I suggest looking at Primer 3 for ways to calculate these values. Or simply call Primer3 rather than reinventing this in your app - it gives this information anyway.

Good luck

yet another biologist hacking perl....

• Comment on Re: how to make this code more efficient? Need to Compute the Min and Max of a Temperature using an equation

Create A New User
Node Status?
node history
Node Type: note [id://835851]
help
Chatterbox?
 LanX will keep his pope tshirt at home [Eily]: what, not a Bishop t-shirt ? [holli]: a pope t-shirt? are the words "fuck the" printed on it? [Eily]: ... [Eily]: I just realized that there's an onion in your eye

How do I use this? | Other CB clients
Other Users?
Others meditating upon the Monastery: (13)
As of 2017-12-14 16:06 GMT
Sections?
Information?
Find Nodes?
Leftovers?
Voting Booth?
What programming language do you hate the most?

Results (397 votes). Check out past polls.

Notices?