Take a look at the dice function as well. You can build up the whole PDL matrix with the input data, then extract the only columns you're interested into in the order you want.