Variable
|
SKULL
|
LIMBS
|
---|
ln L
|
AlCc
|
∆AlCc
|
ln L
|
AlCc
|
∆AlCc
|
---|
BM1
|
395.321
|
−851.309
|
100.238
|
586.680
|
−1234.026
|
117.504
|
BMM
|
437.403
|
−925.206
|
26.341
|
645.680
|
−1341.761
|
9.769
|
OU1
|
414.072
|
−878.544
|
73.003
|
611.689
|
−1273.778
|
77.752
|
OUM
|
427.652
|
−904.941
|
46.606
|
629.454
|
−1308.544
|
42.986
|
EB
|
392.314
|
−844.888
|
106.659
|
583.708
|
−1227.677
|
123.853
|
BM_EB
|
401.074
|
−862.409
|
89.138
|
586.723
|
−1233.707
|
117.823
|
EB_BM
|
400.259
|
−860.780
|
90.767
|
586.626
|
−1233.512
|
118.018
|
BM_EBi
|
447.326
|
−944.969
|
6.578
|
649.317
|
−1348.950
|
2.580
|
EB_BMi
|
450.615
|
−951.547
|
0
|
650.606
|
−1351.530
|
0
|
BM_OU
|
445.226
|
−940.852
|
10.695
|
620.337
|
−1291.075
|
60.455
|
OU_BM
|
425.734
|
−901.867
|
49.680
|
618.049
|
−1286.499
|
65.031
|
BM_OUi
|
448.231
|
−943.881
|
7.666
|
645.977
|
−1339.372
|
12.158
|
OU_BMi
|
330.765
|
−708.949
|
242.598
|
404.802
|
−857.024
|
494.506
|
EB_OU
|
408.027
|
−866.371
|
85.177
|
598.161
|
−1246.638
|
104.892
|
OU_EB
|
430.634
|
−911.585
|
39.962
|
617.650
|
−1285.623
|
65.907
|
EB_OUi
|
446.735
|
−940.855
|
10.692
|
649.460
|
−1346.304
|
5.226
|
OU_EBi
|
332.339
|
−712.063
|
239.484
|
407.160
|
−861.704
|
489.826
|
- We tested the fit of the following evolutionary models: BM1 = Brownian Motion (unique rate), BMM = Brownian Motion (multiple rates), EB = Early Burst, and 12 evolutionary models with shifts from one model to another (e.g. BM_EB = shift of a BM to EB process, EB_BM = shift of EB to BM, BM_EBi = BM_EB with independent rates, EB_BMi = EB_BM with independent rates, etc.). Analyses were performed in R using the functions mvBM(), mvOU(), mvEB() and mvSHIFT() from the R package mvMORPH (Clavel et al., [71])