(* Simulated Data for mcmScheffeContrasts *)



dataset1a={{1, 4.806}, {1, 4.106}, {1, 4.33}, {1, 4.197}, {1, 3.653}, {1, 3.278}, {1, 

    3.54}, {1, 3.046}, {1, 2.929}, {1, 2.109}, {2, 4.948}, {2, 6.387}, {2, 

    5.582}, {2, 6.554}, {2, 6.538}, {2, 6.092}, {2, 5.775}, {2, 4.726}, {2, 

    7.797}, {2, 6.548}, {3, 1.444}, {3, 1.16}, {3, 3.227}, {3, 1.643}, {3, 

    1.77}, {3, -0.2487}, {3, 1.522}, {3, 0.6088}, {3, 2.384}, {3, 2.049}, {4, 

    8.262}, {4, 8.045}, {4, 7.827}, {4, 7.399}, {4, 7.1}, {4, 7.514}, {4, 

    8.774}, {4, 9.083}, {4, 7.889}, {4, 7.993}}