
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:16 04/04/2011

Input data from file "683-tomodel-phylip/chr6_3144.phylip" (7 taxa, 736 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


(Output sorted by AICc. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AICc		wAICc	cwAICc
GTRI	20	-1622.5450	3286.2649	0.4151	0.4151
GTRG	20	-1622.5887	3286.3522	0.3974	0.8126
GTRIG	21	-1622.5887	3288.4716	0.1377	0.9503
HKYI	16	-1629.6983	3292.1533	0.0219	0.9722
HKYG	16	-1629.7585	3292.2737	0.0206	0.9927
HKYIG	17	-1629.7586	3294.3695	0.0072	0.9999
SYMI	17	-1636.1252	3307.1027	0.0000	1.0000
SYMG	17	-1636.1636	3307.1796	0.0000	1.0000
K2PI	13	-1640.5807	3307.6655	0.0000	1.0000
K2PG	13	-1640.6183	3307.7409	0.0000	1.0000
SYMIG	18	-1636.1636	3309.2813	0.0000	1.0000
K2PIG	14	-1640.6184	3309.8193	0.0000	1.0000
GTR	19	-1637.2553	3313.5720	0.0000	1.0000
HKY	15	-1644.3301	3319.3268	0.0000	1.0000
F81I	15	-1649.1677	3329.0021	0.0000	1.0000
F81G	15	-1649.2292	3329.1251	0.0000	1.0000
F81IG	16	-1649.2293	3331.2152	0.0000	1.0000
SYM	16	-1650.9958	3334.7482	0.0000	1.0000
K2P	12	-1655.7019	3335.8354	0.0000	1.0000
JC69I	12	-1659.8731	3344.1777	0.0000	1.0000
JC69G	12	-1659.9204	3344.2724	0.0000	1.0000
JC69IG	13	-1659.9205	3346.3452	0.0000	1.0000
F81	14	-1663.6505	3355.8835	0.0000	1.0000
JC69	11	-1674.6602	3371.6850	0.0000	1.0000



(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1622.5450	3285.0901	0.4154	0.4154
GTRG	20	-1622.5887	3285.1774	0.3977	0.8131
GTRIG	21	-1622.5887	3287.1774	0.1463	0.9594
HKYI	16	-1629.6983	3291.3967	0.0177	0.9771
HKYG	16	-1629.7585	3291.5171	0.0167	0.9938
HKYIG	17	-1629.7586	3293.5171	0.0061	1.0000
SYMI	17	-1636.1252	3306.2503	0.0000	1.0000
SYMG	17	-1636.1636	3306.3273	0.0000	1.0000
K2PI	13	-1640.5807	3307.1614	0.0000	1.0000
K2PG	13	-1640.6183	3307.2367	0.0000	1.0000
SYMIG	18	-1636.1636	3308.3273	0.0000	1.0000
K2PIG	14	-1640.6184	3309.2368	0.0000	1.0000
GTR	19	-1637.2553	3312.5105	0.0000	1.0000
HKY	15	-1644.3301	3318.6602	0.0000	1.0000
F81I	15	-1649.1677	3328.3354	0.0000	1.0000
F81G	15	-1649.2292	3328.4585	0.0000	1.0000
F81IG	16	-1649.2293	3330.4586	0.0000	1.0000
SYM	16	-1650.9958	3333.9916	0.0000	1.0000
K2P	12	-1655.7019	3335.4038	0.0000	1.0000
JC69I	12	-1659.8731	3343.7462	0.0000	1.0000
JC69G	12	-1659.9204	3343.8408	0.0000	1.0000
JC69IG	13	-1659.9205	3345.8410	0.0000	1.0000
F81	14	-1663.6505	3355.3010	0.0000	1.0000
JC69	11	-1674.6602	3371.3204	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1629.6983	3365.0163	0.3662	0.3662
HKYG	16	-1629.7585	3365.1367	0.3448	0.7110
K2PI	13	-1640.5807	3366.9774	0.1374	0.8484
K2PG	13	-1640.6183	3367.0527	0.1323	0.9806
HKYIG	17	-1629.7586	3371.7380	0.0127	0.9934
K2PIG	14	-1640.6184	3373.6540	0.0049	0.9982
GTRI	20	-1622.5450	3377.1147	0.0009	0.9991
GTRG	20	-1622.5887	3377.2020	0.0008	0.9999
GTRIG	21	-1622.5887	3383.8033	0.0000	1.0000
SYMI	17	-1636.1252	3384.4712	0.0000	1.0000
SYMG	17	-1636.1636	3384.5482	0.0000	1.0000
HKY	15	-1644.3301	3387.6786	0.0000	1.0000
K2P	12	-1655.7019	3390.6186	0.0000	1.0000
SYMIG	18	-1636.1636	3391.1494	0.0000	1.0000
F81I	15	-1649.1677	3397.3539	0.0000	1.0000
F81G	15	-1649.2292	3397.4769	0.0000	1.0000
JC69I	12	-1659.8731	3398.9609	0.0000	1.0000
JC69G	12	-1659.9204	3399.0556	0.0000	1.0000
GTR	19	-1637.2553	3399.9339	0.0000	1.0000
F81IG	16	-1649.2293	3404.0783	0.0000	1.0000
JC69IG	13	-1659.9205	3405.6570	0.0000	1.0000
SYM	16	-1650.9958	3407.6113	0.0000	1.0000
F81	14	-1663.6505	3419.7182	0.0000	1.0000
JC69	11	-1674.6602	3421.9339	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_3144.phylip = 1 - 736;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3144.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_3144.phylip = 1 - 736;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3144.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_3144.phylip = 1 - 736;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3144.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

