In reply to Darius Baer.

Read in the ID variable for the second group. Then ID[idx] contains the IDs. You might need to reshape the IDs. For the example in the blog post, you would use

ID = T("A":"D");
M = shape(ID[idx], 0, k);
print M;