mvrlm_coef_ref1 <- list(mrps51 = structure(list(coef = c(259.03890173607, 6.68154143619772, 
23.9866059982592, 10.6150999341518, -13.1592307881691, 2.75824176909601, 
-4.4875559425223, -20.678112855744), se = c(2.2816468792565, 
1.42503510822531, 2.06828746436643, 2.52861277942036, 5.02521090016233, 
5.17216772870593, 9.00670970004303, 8.76937346644863), t = c(113.531547800456, 
4.68868549106743, 11.5973269729249, 4.19799346920373, -2.61864249075477, 
0.533285445053833, -0.498245873573661, -2.35799204297295), dof = c(42.829508150211, 
52.5940440614573, 67.8238573048457, 64.7283647725202, 12.1012868544905, 
21.2584602568237, 24.5184968140515, 10.8662204573521), `Pr(>|t|)` = c(0, 
1.9895753584187e-05, 0, 8.38739354385165e-05, 0.0223094186142752, 
0.599367742342138, 0.622750819116767, 0.0382105140913368)), row.names = c("(Intercept)", 
"dsexFemale", "m072801B *", "m072801C", "m072801D", "m072801E", 
"m072801Omitted", "m072801Not Reached"), class = "data.frame"), 
    mrps22 = structure(list(coef = c(247.935794113977, 1.88046014566052, 
    32.801512606725, 14.6130426942565, -8.06400341685136, 5.48361564516263, 
    -5.88650246936558, -48.0796618295349), se = c(3.04597949712119, 
    2.14117856700418, 3.19097876218864, 3.4320927991799, 7.27180355926605, 
    7.01239056397551, 10.5900323482414, 25.9085640074414), t = c(81.3977225875308, 
    0.878236021338268, 10.2794518708194, 4.25776444557336, -1.10894131712013, 
    0.78198947921318, -0.555853115060892, -1.85574398549165), 
        dof = c(24.3776367361704, 32.1881924749892, 34.797025669306, 
        27.6564106476932, 12.4467345541906, 15.0290083311421, 
        21.9247781460483, 3.62323768362101), `Pr(>|t|)` = c(0, 
        0.386324471922378, 4.38671321489892e-12, 0.000214184359123237, 
        0.288423026803619, 0.44636662205117, 0.583939770340537, 
        0.144473374512257)), row.names = c("(Intercept)", "dsexFemale", 
    "m072801B *", "m072801C", "m072801D", "m072801E", "m072801Omitted", 
    "m072801Not Reached"), class = "data.frame"))


mvrlm_resid_cov1 <- structure(c(4011562.87426289, 3987484.87970119, 3987484.87970119, 
6844469.82478946), dim = c(2L, 2L), dimnames = list(c("mrps51", 
"mrps22"), c("mrps51", "mrps22")))

mvrlm_coef_ref2 <- list(composite = structure(list(coef = c(255.691490082219, 4.72169682680796, 
25.8816446607727, 10.9443482381984, -11.998428885648, 1.71037227388718, 
-6.12840564532562, -44.3173851163579), se = c(2.36687451613302, 
1.61951337241073, 2.21448992130905, 3.07604760806762, 6.30547962509871, 
5.44451072679053, 8.45062633868229, 28.6379333320317), t = c(108.029170257816, 
2.91550345137284, 11.6874068433209, 3.55792550462951, -1.90285745082559, 
0.314146184976923, -0.725201351913191, -1.54750640007911), dof = c(38.4313229420013, 
36.6531150167032, 58.5413906063462, 51.5507603100821, 12.6185974933964, 
22.3400034235907, 26.4220838567071, 2.49076445460829), `Pr(>|t|)` = c(0, 
0.00602508669235746, 0, 0.000812455180321692, 0.0801221625675621, 
0.756324868698339, 0.4747028035334, 0.237325370493719)), row.names = c("(Intercept)", 
"dsexFemale", "m072801B *", "m072801C", "m072801D", "m072801E", 
"m072801Omitted", "m072801Not Reached"), class = "data.frame"), 
    mrps22 = structure(list(coef = c(247.935794113977, 1.88046014566052, 
    32.801512606725, 14.6130426942565, -8.06400341685136, 5.48361564516263, 
    -5.88650246936558, -48.0796618295349), se = c(3.04597949712119, 
    2.14117856700418, 3.19097876218864, 3.4320927991799, 7.27180355926605, 
    7.01239056397551, 10.5900323482414, 25.9085640074414), t = c(81.3977225875308, 
    0.878236021338268, 10.2794518708194, 4.25776444557336, -1.10894131712013, 
    0.78198947921318, -0.555853115060892, -1.85574398549165), 
        dof = c(24.3776367361704, 32.1881924749892, 34.797025669306, 
        27.6564106476932, 12.4467345541906, 15.0290083311421, 
        21.9247781460483, 3.62323768362101), `Pr(>|t|)` = c(0, 
        0.386324471922378, 4.38671321489892e-12, 0.000214184359123237, 
        0.288423026803619, 0.44636662205117, 0.583939770340537, 
        0.144473374512257)), row.names = c("(Intercept)", "dsexFemale", 
    "m072801B *", "m072801C", "m072801D", "m072801E", "m072801Omitted", 
    "m072801Not Reached"), class = "data.frame"))

mvrlm_residpv2 <- list(structure(c(56.3414319138114, 67.6811760321951, 31.6214319138114, 
31.3714319138114, 16.8214319138114, 39.9814319138114, 52.2402771197727, 
40.1865338230753, 37.1902771197727, 18.7002771197727, 28.2802771197727, 
42.2702771197727, 42.8348609592637, 46.0938239385802, 29.4148609592636, 
27.1048609592636, 35.9148609592636, 43.9548609592637, 66.3440065007161, 
51.5220073016552, 32.1340065007162, 14.9740065007161, 29.0940065007161, 
48.9940065007162, 55.7152656574385, 45.7407851895365, 22.9152656574385, 
22.2052656574385, 30.0852656574385, 43.3852656574385), dim = 6:5, dimnames = list(
    NULL, c("mrpcm1", "mrpcm2", "mrpcm3", "mrpcm4", "mrpcm5"))), 
    structure(c(69.7922331336375, 51.3626932792981, 58.6822331336375, 
    36.8522331336375, 22.5522331336375, 44.4422331336375, 69.7922331336375, 
    51.3626932792981, 58.6822331336375, 36.8522331336375, 22.5522331336375, 
    44.4422331336375, 69.7922331336375, 51.3626932792981, 58.6822331336375, 
    36.8522331336375, 22.5522331336375, 44.4422331336375, 69.7922331336375, 
    51.3626932792981, 58.6822331336375, 36.8522331336375, 22.5522331336375, 
    44.4422331336375, 69.7922331336375, 51.3626932792981, 58.6822331336375, 
    36.8522331336375, 22.5522331336375, 44.4422331336375), dim = 6:5, dimnames = list(
        NULL, c("mrps22", "mrps22", "mrps22", "mrps22", "mrps22"
        ))))

mvrlm_resid_cov2 <- structure(c(3637215.64625793, 4395877.60747972, 4395877.60747972, 
6844469.82478946), dim = c(2L, 2L), dimnames = list(c("composite", 
"mrps22"), c("composite", "mrps22")))

mvrlm_coef_ref3 <- list(algebra = structure(list(coef = c(258.600210017001, 6.49222474005324, 
24.7391227972011, 11.680971281142, -12.8871534885896, 1.98741028174384, 
-5.3110754662388, -33.4928541526612), se = c(2.38014319215483, 
1.59643140022273, 2.19055814051215, 3.04655538457315, 6.58913052399054, 
5.24356370859732, 9.28127572620601, 23.4862810605838), t = c(108.649013584297, 
4.06671075196057, 11.2935248508935, 3.83415687772852, -1.95582003447472, 
0.379019001616266, -0.572235501122199, -1.42606034843341), dof = c(41.7238003891765, 
44.4749535721814, 58.8980282326851, 63.7672250667681, 13.5860296321186, 
22.8602630567283, 23.4053209178299, 2.73505264533492), `Pr(>|t|)` = c(0, 
0.000191951299627124, 2.22044604925031e-16, 0.000291059122298876, 
0.0713696694435905, 0.708171355984658, 0.572618435850891, 0.257359883043741
)), row.names = c("(Intercept)", "dsexFemale", "m072801B *", 
"m072801C", "m072801D", "m072801E", "m072801Omitted", "m072801Not Reached"
), class = "data.frame"), geometry = structure(list(coef = c(255.501195547217, 
5.15869213348664, 22.3457824869964, 8.80889910534128, -9.26139075681579, 
-0.174882805590771, -4.71382186560159, -31.7661493980583), se = c(2.44292791904898, 
1.64604974410557, 2.31468229061241, 3.64791490343731, 5.44020835677012, 
5.98021963077922, 7.37914196558784, 24.8394292430355), t = c(104.588102479373, 
3.13398313262384, 9.65393072631332, 2.41477647876077, -1.70239633290706, 
-0.029243542275718, -0.638803520461348, -1.27885987585503), dof = c(35.3688352916959, 
36.7297442383663, 53.9657891394272, 48.9980631621269, 13.5180224951499, 
26.0264650878676, 23.1254044077314, 1.69566886861829), `Pr(>|t|)` = c(0, 
0.00338260070611218, 2.35811370430383e-13, 0.0195230234936867, 
0.111552374437042, 0.976893448896208, 0.52922752340554, 0.348275076610693
)), row.names = c("(Intercept)", "dsexFemale", "m072801B *", 
"m072801C", "m072801D", "m072801E", "m072801Omitted", "m072801Not Reached"
), class = "data.frame"))

mvrlm_residpv3 <- list(structure(c(60.4029508294734, 68.6844922656711, 28.4829508294734, 
37.8529508294734, 12.7129508294734, 40.2129508294734, 63.5562350488423, 
42.0818606246983, 30.8262350488422, 34.8962350488422, 17.9762350488422, 
41.7062350488422, 43.2150308481005, 38.4519050568938, 20.6150308481004, 
35.3650308481004, 24.5150308481004, 36.5950308481005, 71.9341112495783, 
45.6512171425833, 26.9541112495783, 22.9841112495782, 33.8641112495782, 
48.7541112495783, 54.7038842527286, 44.7438608391427, 15.9038842527286, 
28.2438842527287, 35.2438842527287, 45.8338842527286), dim = 6:5, dimnames = list(
    NULL, c("mrps51", "mrps52", "mrps53", "mrps54", "mrps55"))), 
    structure(c(70.4928099389219, 48.9632469644581, 20.0528099389219, 
    22.3228099389219, 7.64280993892186, 37.9028099389218, 37.7053737798508, 
    24.8521324685974, 29.6753737798508, 17.5453737798508, 35.5653737798509, 
    42.8853737798509, 48.9743370184677, 43.2189366239021, 31.1443370184677, 
    23.5443370184677, 32.6443370184677, 39.1643370184677, 72.0158610772791, 
    55.7385275118879, 39.775861077279, 6.14586107727905, 15.525861077279, 
    41.505861077279, 81.4232673469813, 40.8522662600885, 23.1032673469813, 
    13.4032673469814, 27.6932673469813, 35.9132673469813), dim = 6:5, dimnames = list(
        NULL, c("mrps31", "mrps32", "mrps33", "mrps34", "mrps35"
        ))))

mvrlm_resid_cov3 <- structure(c(3567321.0079866, 2989766.25340655, 2989766.25340655, 
3099735.84516655), dim = c(2L, 2L), dimnames = list(c("algebra", 
"geometry"), c("algebra", "geometry")))

mvrlm_coef_ref4 <- list(algebra = structure(list(coef = c(258.600210017001, 6.49222474005324, 
24.7391227972011, 11.680971281142, -12.8871534885896, 1.98741028174384, 
-5.3110754662388, -33.4928541526612), se = c(2.38014319215483, 
1.59643140022273, 2.19055814051215, 3.04655538457315, 6.58913052399054, 
5.24356370859732, 9.28127572620601, 23.4862810605838), t = c(108.649013584297, 
4.06671075196057, 11.2935248508935, 3.83415687772852, -1.95582003447472, 
0.379019001616266, -0.572235501122199, -1.42606034843341), dof = c(41.7238003891765, 
44.4749535721814, 58.8980282326851, 63.7672250667681, 13.5860296321186, 
22.8602630567283, 23.4053209178299, 2.73505264533492), `Pr(>|t|)` = c(0, 
0.000191951299627124, 2.22044604925031e-16, 0.000291059122298876, 
0.0713696694435905, 0.708171355984658, 0.572618435850891, 0.257359883043741
)), row.names = c("(Intercept)", "dsexFemale", "m072801B *", 
"m072801C", "m072801D", "m072801E", "m072801Omitted", "m072801Not Reached"
), class = "data.frame"), geometry = structure(list(coef = c(255.501195547217, 
5.15869213348664, 22.3457824869964, 8.80889910534128, -9.26139075681579, 
-0.174882805590771, -4.71382186560159, -31.7661493980583), se = c(2.44292791904898, 
1.64604974410557, 2.31468229061241, 3.64791490343731, 5.44020835677012, 
5.98021963077922, 7.37914196558784, 24.8394292430355), t = c(104.588102479373, 
3.13398313262384, 9.65393072631332, 2.41477647876077, -1.70239633290706, 
-0.029243542275718, -0.638803520461348, -1.27885987585503), dof = c(35.3688352916959, 
36.7297442383663, 53.9657891394272, 48.9980631621269, 13.5180224951499, 
26.0264650878676, 23.1254044077314, 1.69566886861829), `Pr(>|t|)` = c(0, 
0.00338260070611218, 2.35811370430383e-13, 0.0195230234936867, 
0.111552374437042, 0.976893448896208, 0.52922752340554, 0.348275076610693
)), row.names = c("(Intercept)", "dsexFemale", "m072801B *", 
"m072801C", "m072801D", "m072801E", "m072801Omitted", "m072801Not Reached"
), class = "data.frame"), measurement = structure(list(coef = c(248.576438252752, 
2.52127622080741, 32.3417305093413, 12.3962876016519, -14.1713411898605, 
3.12917121690047, -5.25885207512995, -37.1845446310284), se = c(3.03889863244681, 
2.29688961181709, 3.01882011632236, 4.31525292916813, 9.26991690622827, 
7.3908251358158, 10.1585009416688, 31.2308992356816), t = c(81.7982000447997, 
1.0976915076092, 10.7133678931295, 2.87266767559825, -1.52874522320034, 
0.423385908798811, -0.517679931844946, -1.19063317230855), dof = c(35.5443984857029, 
32.3078726249878, 39.178675851679, 37.0424297422304, 15.3035448270071, 
17.3614807133915, 23.5615918607897, 2.71962293797103), `Pr(>|t|)` = c(0, 
0.280453169025543, 3.31290550548147e-13, 0.00669895295156353, 
0.146731923709766, 0.677212228473132, 0.609508356158105, 0.327281552297339
)), row.names = c("(Intercept)", "dsexFemale", "m072801B *", 
"m072801C", "m072801D", "m072801E", "m072801Omitted", "m072801Not Reached"
), class = "data.frame"))

mvrlm_residpv4 <-list(structure(c(60.4029508294734, 68.6844922656711, 28.4829508294734, 
37.8529508294734, 12.7129508294734, 40.2129508294734, 63.5562350488423, 
42.0818606246983, 30.8262350488422, 34.8962350488422, 17.9762350488422, 
41.7062350488422, 43.2150308481005, 38.4519050568938, 20.6150308481004, 
35.3650308481004, 24.5150308481004, 36.5950308481005, 71.9341112495783, 
45.6512171425833, 26.9541112495783, 22.9841112495782, 33.8641112495782, 
48.7541112495783, 54.7038842527286, 44.7438608391427, 15.9038842527286, 
28.2438842527287, 35.2438842527287, 45.8338842527286), dim = 6:5, dimnames = list(
    NULL, c("mrps51", "mrps52", "mrps53", "mrps54", "mrps55"))), 
    structure(c(70.4928099389219, 48.9632469644581, 20.0528099389219, 
    22.3228099389219, 7.64280993892186, 37.9028099389218, 37.7053737798508, 
    24.8521324685974, 29.6753737798508, 17.5453737798508, 35.5653737798509, 
    42.8853737798509, 48.9743370184677, 43.2189366239021, 31.1443370184677, 
    23.5443370184677, 32.6443370184677, 39.1643370184677, 72.0158610772791, 
    55.7385275118879, 39.775861077279, 6.14586107727905, 15.525861077279, 
    41.505861077279, 81.4232673469813, 40.8522662600885, 23.1032673469813, 
    13.4032673469814, 27.6932673469813, 35.9132673469813), dim = 6:5, dimnames = list(
        NULL, c("mrps31", "mrps32", "mrps33", "mrps34", "mrps35"
        ))), structure(c(52.5560805346176, 87.285128264762, 60.1760805346177, 
    57.3360805346177, 11.6560805346177, 52.4360805346176, 69.7922331336375, 
    51.3626932792981, 58.6822331336375, 36.8522331336375, 22.5522331336375, 
    44.4422331336375, 59.2679638986813, 73.3635928816054, 55.1579638986813, 
    66.1879638986813, 56.1779638986813, 66.9379638986813, 69.3687956844132, 
    65.6263564272836, 55.0687956844132, 16.6687956844132, 42.0387956844132, 
    62.7987956844132, 63.7177018341451, 61.6913853365828, 46.9277018341451, 
    32.9977018341451, 37.8377018341451, 65.5177018341451), dim = 6:5, dimnames = list(
        NULL, c("mrps21", "mrps22", "mrps23", "mrps24", "mrps25"
        ))))

mvrlm_resid_cov4 <- structure(c(3567321.0079866, 2989766.25340655, 4056151.00854307, 
2989766.25340655, 3099735.84516655, 3810011.98338393, 4056151.00854307, 
3810011.98338393, 5851003.07925973), dim = c(3L, 3L), dimnames = list(
    c("algebra", "geometry", "measurement"), c("algebra", "geometry", 
    "measurement")))