Perl Monk, Perl Meditation 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?
 [stevieb]: I don't know what you mean. I write the code in C/C++ that accesses the ICs per the datasheet, wrap it for Perl, then use them. My code ensures that the functions/methods from Perl do the right thing for the chips/sensors etc [stevieb]: At least that's what I'm working on right now. I've also been updating my automated test software Test::BrewBuild so it handles this work properly. It's not really easy meshing high level languages with low-level hardware :)

How do I use this? | Other CB clients
Other Users?
Others chilling in the Monastery: (6)
As of 2017-06-25 23:08 GMT
Sections?
Information?
Find Nodes?
Leftovers?
Voting Booth?
How many monitors do you use while coding?

Results (572 votes). Check out past polls.