Sorry, there is no functionality like that for mask selection, but there are some alternatives. First you could parse the output of what you've already done with 1dcat and column selectors (that can skip columns). The skipping of sub-bricks in 3D+t datasets and columns or rows in 1D files is described in the help for 3dcalc and 1dcat.

3dmaskave can take a mask range, but it still has to be consecutive values and that would also require juggling of the mask values. You could keep all your masks in separate sub-bricks and then combine selected values into a single sub-brick volume.

I would probably stick with a 3dcalc combination of the ROIs you want in a temporary dataset. This could be done in an in-line 3dcalc command inside the 3dROIstats command. It's a bit more difficult to read, but it should work. This method is given as an example in the Misc. Utilities (Jewel Box) presentation.

afni.nimh.nih.gov/pub/dist/edu/latest/afni12_miscutil