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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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	-1631.5624	3304.1657	0.4386	0.4386
GTRI	20	-1631.7599	3304.5606	0.3600	0.7987
GTRIG	21	-1631.3934	3305.9332	0.1813	0.9799
GTR	19	-1636.7091	3312.3589	0.0073	0.9872
HKYG	16	-1640.1099	3312.8905	0.0056	0.9928
HKYI	16	-1640.3113	3313.2934	0.0046	0.9974
HKYIG	17	-1639.9341	3314.6237	0.0024	0.9998
HKY	15	-1645.1185	3320.8282	0.0001	0.9999
F81G	15	-1645.6451	3321.8814	0.0001	0.9999
F81I	15	-1645.8251	3322.2413	0.0001	1.0000
F81IG	16	-1645.5185	3323.7078	0.0000	1.0000
F81	14	-1650.3851	3329.2869	0.0000	1.0000
SYMG	17	-1659.2496	3353.2548	0.0000	1.0000
SYMI	17	-1659.4462	3353.6480	0.0000	1.0000
SYMIG	18	-1659.0682	3354.9819	0.0000	1.0000
SYM	16	-1664.3709	3361.4126	0.0000	1.0000
K2PG	13	-1668.1317	3362.7105	0.0000	1.0000
K2PI	13	-1668.3186	3363.0845	0.0000	1.0000
K2PIG	14	-1667.9733	3364.4633	0.0000	1.0000
JC69G	12	-1672.7113	3369.8055	0.0000	1.0000
JC69I	12	-1672.8785	3370.1399	0.0000	1.0000
K2P	12	-1673.0779	3370.5386	0.0000	1.0000
JC69IG	13	-1672.5916	3371.6304	0.0000	1.0000
JC69	11	-1677.4599	3377.2434	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1631.5624	3303.1248	0.4354	0.4354
GTRI	20	-1631.7599	3303.5197	0.3574	0.7928
GTRIG	21	-1631.3934	3304.7868	0.1897	0.9825
GTR	19	-1636.7091	3311.4183	0.0069	0.9894
HKYG	16	-1640.1099	3312.2197	0.0046	0.9940
HKYI	16	-1640.3113	3312.6226	0.0038	0.9978
HKYIG	17	-1639.9341	3313.8682	0.0020	0.9998
HKY	15	-1645.1185	3320.2370	0.0001	0.9999
F81G	15	-1645.6451	3321.2903	0.0000	0.9999
F81I	15	-1645.8251	3321.6502	0.0000	1.0000
F81IG	16	-1645.5185	3323.0370	0.0000	1.0000
F81	14	-1650.3851	3328.7703	0.0000	1.0000
SYMG	17	-1659.2496	3352.4993	0.0000	1.0000
SYMI	17	-1659.4462	3352.8924	0.0000	1.0000
SYMIG	18	-1659.0682	3354.1364	0.0000	1.0000
SYM	16	-1664.3709	3360.7418	0.0000	1.0000
K2PG	13	-1668.1317	3362.2634	0.0000	1.0000
K2PI	13	-1668.3186	3362.6373	0.0000	1.0000
K2PIG	14	-1667.9733	3363.9467	0.0000	1.0000
JC69G	12	-1672.7113	3369.4226	0.0000	1.0000
JC69I	12	-1672.8785	3369.7571	0.0000	1.0000
K2P	12	-1673.0779	3370.1558	0.0000	1.0000
JC69IG	13	-1672.5916	3371.1833	0.0000	1.0000
JC69	11	-1677.4599	3376.9199	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1640.1099	3387.7239	0.4333	0.4333
HKYI	16	-1640.3113	3388.1268	0.3543	0.7876
HKY	15	-1645.1185	3391.0222	0.0833	0.8708
F81G	15	-1645.6451	3392.0755	0.0492	0.9200
F81I	15	-1645.8251	3392.4354	0.0411	0.9611
HKYIG	17	-1639.9341	3394.0914	0.0180	0.9791
F81	14	-1650.3851	3394.8364	0.0124	0.9914
GTRG	20	-1631.5624	3397.5050	0.0033	0.9947
GTRI	20	-1631.7599	3397.9000	0.0027	0.9974
F81IG	16	-1645.5185	3398.5412	0.0019	0.9993
GTR	19	-1636.7091	3401.0795	0.0005	0.9999
GTRIG	21	-1631.3934	3403.8861	0.0001	1.0000
K2PG	13	-1668.1317	3423.6105	0.0000	1.0000
K2PI	13	-1668.3186	3423.9845	0.0000	1.0000
JC69G	12	-1672.7113	3426.0508	0.0000	1.0000
JC69I	12	-1672.8785	3426.3853	0.0000	1.0000
K2P	12	-1673.0779	3426.7839	0.0000	1.0000
JC69	11	-1677.4599	3428.8290	0.0000	1.0000
K2PIG	14	-1667.9733	3430.0129	0.0000	1.0000
JC69IG	13	-1672.5916	3432.5304	0.0000	1.0000
SYMG	17	-1659.2496	3432.7225	0.0000	1.0000
SYMI	17	-1659.4462	3433.1156	0.0000	1.0000
SYM	16	-1664.3709	3436.2460	0.0000	1.0000
SYMIG	18	-1659.0682	3439.0787	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/chr6_7243.phylip = 1 - 828;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7243.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/chr6_7243.phylip = 1 - 828;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7243.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_7243.phylip = 1 - 828;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7243.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

