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

**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