Simulate the data consiting of three sided heterogeniety

Simulate.data.threeSided(model_test)