Key: Difference between revisions
(Problems) |
|||
Line 67: | Line 67: | ||
</source> | </source> | ||
However, while this code works, it is inefficient in that the unique is computed twice. This can be avoided by letting Key return the unique and using that: | However, while this code works, it is inefficient in that the unique is computed twice. This can be avoided by letting Key return the unique and using that: | ||
<source> | <source lang=apl> | ||
(keys counts)←,⌿{⍺,≢⍵}⌸'TCCGCGGTGGCG' | (keys counts)←,⌿{⍺,≢⍵}⌸'TCCGCGGTGGCG' | ||
keys⌷⍨⊃⍒counts | keys⌷⍨⊃⍒counts | ||
Line 73: | Line 73: | ||
</source> | </source> | ||
Unfortunately, this can introduce a different inefficiency, in that the result of Key's operand can end up being a [[heterogeneous array]] (containing multiple [[datatype]]s), and these are stored as pointer arrays, consuming memory for one pointer per element, and forcing "pointer chasing" when addressing the data. A possible work-around is to collect the unique keys separately from the result of counts: | Unfortunately, this can introduce a different inefficiency, in that the result of Key's operand can end up being a [[heterogeneous array]] (containing multiple [[datatype]]s), and these are stored as pointer arrays, consuming memory for one pointer per element, and forcing "pointer chasing" when addressing the data. A possible work-around is to collect the unique keys separately from the result of counts: | ||
<source> | <source lang=apl> | ||
data←'TCCGCGGTGGCG' | data←'TCCGCGGTGGCG' | ||
keys←0⌿data | keys←0⌿data |
Revision as of 09:28, 16 January 2022
⌸
|
Key (⌸
) is a primitive monadic operator which takes a dyadic function operand where specified keys group the indices or major cells of an argument. It was introduced in Dyalog APL version 14.0 and is commonly compared to SQL's GROUP BY statement.
Description
Monadically, Key will group identical major cells together and applies the function operand once for each unique major cell. The function is applied with the unique major cell as left argument, while the right argument is the indices of major cells that match it:
{⍺⍵}⌸'Mississippi' ┌─┬────────┐ │M│1 │ ├─┼────────┤ │i│2 5 8 11│ ├─┼────────┤ │s│3 4 6 7 │ ├─┼────────┤ │p│9 10 │ └─┴────────┘
In the dyadic case, Key applies the function to collections of major cells from the right argument corresponding to unique elements of the left argument:
'Mississippi'{⍺⍵}⌸'ABCDEFGHIJK' ┌─┬────┐ │M│A │ ├─┼────┤ │i│BEHK│ ├─┼────┤ │s│CDFG│ ├─┼────┤ │p│IJ │ └─┴────┘
The monadic case, f⌸Y
is equivalent to Y f⌸ ⍳≢Y
.
Problems
Vocabulary
A common problem with Key is the inability to control the order of the result (as Key will use the order of appearance) and the "vocabulary" (as Key will never include information for a major cell that doesn't occur). For example, here we want to count occurrences of the letters A, C, G, T:
{⍺,≢⍵}⌸'TCCGCGGTGGCG' T 2 C 4 G 6
Since A is entirely missing in the argument, it isn't mentioned in the result either. Likewise, the result is mis-ordered due to G and T appearing before the first C. A common solution is to inject the vocabulary before the actual data, and then decrement from the counts:
{⍺,¯1+≢⍵}⌸'ACGT','TCCGCGGTGGCG' A 0 C 4 G 6 T 2
Now that the meaning of each count is known, the operand's left argument can be ignored, and the decrementing can be factored out from the operand:
¯1+{≢⍵}⌸'ACGT','TCCGCGGTGGCG' 0 4 6 2
Computing the unique
Key computes the set of unique major cells. Often, this collection is needed separately from the occurrence information, but can be hard to extract. For example, to get the most frequently occurring letter:
⊃⍒{≢⍵}⌸'TCCGCGGTGGCG' 3
Notice that 3 is the index in the unique set of letters, and so it is tempting to write:
{(⊃⍒{≢⍵}⌸⍵)⌷∪⍵}'TCCGCGGTGGCG' G
However, while this code works, it is inefficient in that the unique is computed twice. This can be avoided by letting Key return the unique and using that:
(keys counts)←,⌿{⍺,≢⍵}⌸'TCCGCGGTGGCG' keys⌷⍨⊃⍒counts G
Unfortunately, this can introduce a different inefficiency, in that the result of Key's operand can end up being a heterogeneous array (containing multiple datatypes), and these are stored as pointer arrays, consuming memory for one pointer per element, and forcing "pointer chasing" when addressing the data. A possible work-around is to collect the unique keys separately from the result of counts:
data←'TCCGCGGTGGCG' keys←0⌿data counts←{keys⍪←⍺ ⋄ ≢⍵}⌸data keys⌷⍨⊃⍒counts G
If there are a large number of unique values, the repeated updating of the accumulating keys
variable can be an issue in itself.
External links
Lessons
Documentation