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

Input data from file "591-tomodel-phylip/chr5_1813.phylip" (7 taxa, 374 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-683.7161	1409.8119	0.3186	0.3186
GTRI	20	-683.7895	1409.9586	0.2961	0.6147
GTR	19	-685.0692	1410.2853	0.2514	0.8661
GTRIG	21	-683.5907	1411.8063	0.1175	0.9836
HKYI	16	-692.1104	1417.7446	0.0060	0.9897
HKYG	16	-692.1136	1417.7509	0.0060	0.9957
HKYIG	17	-692.1099	1419.9390	0.0020	0.9977
HKY	15	-695.1070	1421.5547	0.0009	0.9986
F81I	15	-695.5803	1422.5014	0.0006	0.9992
F81G	15	-695.5816	1422.5040	0.0006	0.9997
F81IG	16	-695.5798	1424.6833	0.0002	0.9999
F81	14	-698.5359	1426.2417	0.0001	1.0000
SYMG	17	-698.9188	1433.5567	0.0000	1.0000
SYMI	17	-698.9366	1433.5923	0.0000	1.0000
SYMIG	18	-698.6799	1435.2866	0.0000	1.0000
SYM	16	-701.9201	1437.3640	0.0000	1.0000
K2PI	13	-705.6575	1438.3261	0.0000	1.0000
K2PG	13	-705.7167	1438.4444	0.0000	1.0000
K2PIG	14	-705.6568	1440.4834	0.0000	1.0000
JC69I	12	-708.8018	1442.4678	0.0000	1.0000
JC69G	12	-708.8249	1442.5140	0.0000	1.0000
K2P	12	-709.0874	1443.0391	0.0000	1.0000
JC69IG	13	-708.8007	1444.6126	0.0000	1.0000
JC69	11	-712.0192	1446.7678	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-683.7161	1407.4323	0.3244	0.3244
GTRI	20	-683.7895	1407.5790	0.3015	0.6259
GTR	19	-685.0692	1408.1384	0.2279	0.8538
GTRIG	21	-683.5907	1409.1813	0.1353	0.9891
HKYI	16	-692.1104	1416.2208	0.0040	0.9931
HKYG	16	-692.1136	1416.2271	0.0040	0.9971
HKYIG	17	-692.1099	1418.2199	0.0015	0.9986
HKY	15	-695.1070	1420.2140	0.0005	0.9991
F81I	15	-695.5803	1421.1606	0.0003	0.9995
F81G	15	-695.5816	1421.1632	0.0003	0.9998
F81IG	16	-695.5798	1423.1595	0.0001	0.9999
F81	14	-698.5359	1425.0718	0.0000	1.0000
SYMG	17	-698.9188	1431.8376	0.0000	1.0000
SYMI	17	-698.9366	1431.8732	0.0000	1.0000
SYMIG	18	-698.6799	1433.3598	0.0000	1.0000
SYM	16	-701.9201	1435.8402	0.0000	1.0000
K2PI	13	-705.6575	1437.3150	0.0000	1.0000
K2PG	13	-705.7167	1437.4333	0.0000	1.0000
K2PIG	14	-705.6568	1439.3135	0.0000	1.0000
JC69I	12	-708.8018	1441.6035	0.0000	1.0000
JC69G	12	-708.8249	1441.6498	0.0000	1.0000
K2P	12	-709.0874	1442.1748	0.0000	1.0000
JC69IG	13	-708.8007	1443.6015	0.0000	1.0000
JC69	11	-712.0192	1446.0385	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-692.1104	1479.0089	0.1952	0.1952
HKYG	16	-692.1136	1479.0152	0.1946	0.3898
HKY	15	-695.1070	1479.0778	0.1886	0.5784
F81	14	-698.5359	1480.0114	0.1183	0.6967
F81I	15	-695.5803	1480.0245	0.1175	0.8142
F81G	15	-695.5816	1480.0270	0.1173	0.9315
GTR	19	-685.0692	1482.6993	0.0308	0.9624
HKYIG	17	-692.1099	1484.9322	0.0101	0.9725
GTRG	20	-683.7161	1485.9174	0.0062	0.9786
F81IG	16	-695.5798	1485.9476	0.0061	0.9847
GTRI	20	-683.7895	1486.0641	0.0057	0.9904
K2PI	13	-705.6575	1488.3303	0.0018	0.9923
K2PG	13	-705.7167	1488.4486	0.0017	0.9940
JC69I	12	-708.8018	1488.6946	0.0015	0.9956
JC69G	12	-708.8249	1488.7408	0.0015	0.9971
JC69	11	-712.0192	1489.2053	0.0012	0.9983
K2P	12	-709.0874	1489.2659	0.0012	0.9994
GTRIG	21	-683.5907	1491.5907	0.0004	0.9998
K2PIG	14	-705.6568	1494.2531	0.0001	0.9999
JC69IG	13	-708.8007	1494.6168	0.0001	1.0000
SYMG	17	-698.9188	1498.5499	0.0000	1.0000
SYMI	17	-698.9366	1498.5855	0.0000	1.0000
SYM	16	-701.9201	1498.6283	0.0000	1.0000
SYMIG	18	-698.6799	1503.9964	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_1813.phylip = 1 - 374;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1813.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_1813.phylip = 1 - 374;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1813.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_1813.phylip = 1 - 374;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1813.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

