Add cluster variable to `ampute()` for multilevel amputation

[hanneoberman]

Is your feature request related to a problem? Please describe. No, but it would help simulators investigating multilevel missing data problems

Describe the solution you'd like An argument in the ampute function to specify a clustering variable

Describe alternatives you've considered Splitting up the data by the clustering variable and applying ampute() on each