Databases
StackOverflow2013
Postsdbo
POHow can I quickly run a function over an array of every possible array of length L with given possible elements?
 

Note that there are some explanatory texts on larger screens.

plurals
  1. POHow can I quickly run a function over an array of every possible array of length L with given possible elements?
    primarykey
    data
    text
    <p>I have a function calc_dG that, for any array corresponding to a short DNA sequence (3 to 15 bases or so), gives me the binding energy of that sequence. Actually, it's just an array lookup. nndG is an array of binding energies for adjacent pairs of bases, and thus the binding energies can be calculated with <code>nndG[4*S[:-1]+S[1:]]</code> when using an a,g,c,t -> 0,1,2,3 way of denoting sequences numerically: this means that arrays of many sequences can be calculated at once very quickly in numpy. </p> <p>I need to find, for a length L, every sequence that both fits some template and results in a binding energy value in a certain range. </p> <p>This is very easy to do with iterators: just iterate through every possible array input, calculate the binding energy, and then record the arrays that are in the range. This, however, is far too slow when implemented in Python (for length 15 with 4 possible values for each element there are 4**15 possible arrays, etc etc). I could use Weave or some other method of implementing it in C, but I'd prefer to find an array-based solution that is simple and fast.</p> <p>For example, if every element has the same possible values (eg, [0,1,2,3]), then generating an array of every possible length L 1D array with those values can be done with <code>lambda x: indices(repeat([4],L)).reshape((L,-1)).transpose()</code>; then I can just do <code>calc_dG( result )</code>, and use result[results that are in the desired range] to get the arrays that I want as a final result. This is much faster than using Python iterators, and likely almost as fast, if not faster, than using C iterators. Unfortunately, it doesn't work for arbitrary templates, and for longer sequences, will run out of memory, as it has to store every possible array in memory before calculating values.</p> <p>Is there some way to do all of this without resorting to C?</p>
    singulars
    1. PO
      singulars
      1. PTAnswer
    1. This table or related slice is empty.
    1. PTQuestion
    1. USSteve Tjoa
    1. UScge
    plurals
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. PO
      singulars
      1. PTAnswer
    1. VO
      singulars
      1. POHow can I quickly run a function over an array of every possible array of length L with given possible elements?
      1. This table or related slice is empty.
      1. VTUpMod
    2. VO
      singulars
      1. POHow can I quickly run a function over an array of every possible array of length L with given possible elements?
      1. This table or related slice is empty.
      1. VTUpMod
    1. COThe question is a bit hard to follow. Regardless, have you discovered Biopython yet? It's great.
      singulars
      1. POHow can I quickly run a function over an array of every possible array of length L with given possible elements?
      1. USSteve Tjoa
    2. COWhat is S? What is indices? Try to reformulate the question for people out of the bio domain.
      singulars
      1. POHow can I quickly run a function over an array of every possible array of length L with given possible elements?
      1. UScyborg
 

Querying!

 
Guidance
A row detail

Detail views are divided into sections. All the information in the data section comes from columns in the selected row. The other sections display data from other, related rows.

Related data can be related in a to-one or a to-many fashion. Captions of data related in a to-many fashion link to a list view showing a filtered view of the table.

Try moving around until you find a non-empty to-many entry and click on the label to get to one. You can move back to the root by clicking on the database name in the header.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload