Re: st: Mata and subsets...

[Matthew Baker <matthew.baker@hunter.cuny.edu>]

Austin --

Thanks a million and apologies for not mentioning moremata. Also, I
think it is safe to say I would not come up with that solution myself,
and it works very nicely!

Best,

Matt Baker

On Thu, Mar 1, 2012 at 5:21 PM, Austin Nichols <austinnichols@gmail.com> wrote:
> Matthew Baker <matthew.baker@hunter.cuny.edu>
> Note the Stata FAQ tells you to give the provenance as -moremata- on SSC.
>
> Here is an unpleasantly Stataish (as opposed to Mataish) solution,
> with the looping you hope to avoid:
>
> cap ssc inst moremata
> mata
> transmorphic matrix cosubsets(n,k) {
>  m1=mm_subsets(n,k)
>  c=runningsum(J(1,n,1))
> st_local("c",invtokens(strofreal(c)))
> st_local("d",invtokens(strofreal(m1'[1,.])))
> stata("loc e: list c-d")
> m2=tokens(st_local("e"))'
>  for (i=2; i<=cols(m1); i++) {
> st_local("d",invtokens(strofreal(m1'[i,.])))
> stata("loc e: list c-d")
> m2=m2,tokens(st_local("e"))'
>  }
> return(m2)
> }
> mm_subsets(5,2)
> cosubsets(5,2)
>
>
> On Thu, Mar 1, 2012 at 3:03 PM, Matthew Baker
> <matthew.baker@hunter.cuny.edu> wrote:
>> Dear Listers --
>>
>> I have a question about obtaining subsets and their complements in
>> mata. To be clear, what I am trying to do is get a listing of both the
>> subsets and their corresponding complements.
>>
>> For example,
>> mata
>> mm_subsets(4,2)
>>
>> gives me the subsets:
>>       1   2   3   4   5   6
>>    +-------------------------+
>>  1 |  1   1   1   2   2   3  |
>>  2 |  2   3   4   3   4   4  |
>>    +-------------------------+
>>
>> But the question is: is there also an (efficient, general) way to
>> obtain the corresponding complements:
>>
>>       1   2   3   4   5   6
>>    +-------------------------+
>>  1 |  3   2   2   1   1   1  |
>>  2 |  4   4   3   4   3   2  |
>>    +-------------------------+
>>
>> I'm hoping for a general method that doesn't depend upon having four
>> items or the two-element subsets but am baffled about how to do this
>> without looping, etc.
>>
>> Thanks in advance!
>>
>> Matt Baker
>
> *
> *   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/



-- 
Dr. Matthew J. Baker
Department of Economics
Hunter College and the Graduate Center, CUNY

*
*   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/