I think it should be updated so you can get a different vector when running with replacement, but I think we should keep the requirement that ``n <=`` the total count, since the name of our method is subsample. It is easy enough for someone to call scipy's multinomial function for other use cases.