Re: st: MATA all combinations / pairs of a row

[Nick Cox <njcoxstata@gmail.com>]

One idea only from me. Look at -select()- and think of rowsums being 2.

Nick

On Wed, Oct 5, 2011 at 10:08 PM, Sebastian Eppner <eppner@uni-potsdam.de> wrote:
> Hi,
> I have N elements. In this example N will be 3. But in my application
> it will be usually around 7, going up to 20 in some cases…
>
> Each element can be either 1 or 0. My basic matrix has all 2^N, here
> 2^3 permutations of the elements:
>
> BASIC = ( 1,1,1) \ (1,1,0) \ (1,0,1) \ (0,1,1) \ (0,0,1) \ (0,1,0) \
> (1,0,0) \(0,0,0)
>
> The elements have 2 “properties” , PROP1 and PROP2  each, stored in 2
> rows of another matrix, for example:
>
> PROPERTIES = (1.4,4.6,8.1) \ (0.2, 0.4, 0.4)
>
> Every row in BASIC resembles a unique combination of the elements.
> What I want to do, is to calculate some indicators for every
> combination, using the Properties of the elements.
>
> An easy example would be to calculate the total sum of PROP1 for each
> combination. I just need to multiply the “diagonalized” first row of
> the PROPERTIES Matrix with the BASIC Matrix and then calculate the
> rowsum of this resulting matrix.
>
>
> What I want to do now is more complex and I am really stuck. The
> indicator I want to calculate now is the following:
>
> For each row/combination in BASIC  I need to  find all PAIRS of Ones
> that are contained in the combination. For any pair, I need to
> multiply:
>
> PROP1 of Element1 * PROP1 of Element2  * abs(PROP2 of Element1-PROP2
> of ELEMENT2)
>
> It is the sum of all PAIRS that I need to know for every row….
>
>
> Of course, not every row has PAIRS in it… all rows with only 1 One
> have no pairs, so the result for these rows should be 0 (or missing, I
> don’t care).
>
> Also, many rows (those with only 2 Ones) have only a single Pair…
>
> In my Example with 3 elements, only the first row sums up more than
> one (3) Pair(s)….
>
>
> I am looking for a way to solve this problem by not using loops… and
> as much matrix algorithms as possible. Since I have a lot of BASIC
> Matrices with more elements than 3… I guess the whole thing would be
> very very time consuming if I started smth with loops… Maybe there is
> no way to get rid of any loop… any idea would be very welcome (even
> some, with a little looping).

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/