add function for creating subsets

[bmacGTPM]

and computing mean, variance, sd.

And possibly plotting, possibly with type = 'subset', although that may be the same as scatter