Consider size argument in simulate_data function to result in multi-indexed output

[AlexanderFengler]

The size argument essentially operates as number of synthetic datasets. To bring this to it's logical conclusion, maybe it should correspond to an index dimension (dataset-n) in a multi-indexed return dataframe?

Alternatively we could simply add a column, dataset-replica or something that accumulates.