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

Input data from file "683-tomodel-phylip/chr1_32461.phylip" (7 taxa, 505 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	-1586.6676	3215.0708	0.4242	0.4242
HKYG	16	-1591.7385	3216.5918	0.1983	0.6224
GTRIG	21	-1586.6475	3217.2080	0.1457	0.7681
GTRI	20	-1587.9774	3217.6903	0.1145	0.8826
HKYIG	17	-1591.7165	3218.6896	0.0695	0.9520
HKYI	16	-1593.1576	3219.4300	0.0480	1.0000
GTR	19	-1601.9439	3243.4549	0.0000	1.0000
HKY	15	-1606.6336	3244.2488	0.0000	1.0000
K2PG	13	-1632.0550	3290.8514	0.0000	1.0000
K2PIG	14	-1631.9783	3292.8137	0.0000	1.0000
K2PI	13	-1633.1536	3293.0486	0.0000	1.0000
F81G	15	-1633.3226	3297.6268	0.0000	1.0000
SYMG	17	-1631.6178	3298.4922	0.0000	1.0000
F81IG	16	-1633.3086	3299.7320	0.0000	1.0000
SYMIG	18	-1631.5767	3300.5607	0.0000	1.0000
SYMI	17	-1632.6780	3300.6126	0.0000	1.0000
F81I	15	-1634.9604	3300.9024	0.0000	1.0000
K2P	12	-1648.4978	3321.6298	0.0000	1.0000
F81	14	-1647.3510	3323.5591	0.0000	1.0000
SYM	16	-1648.0656	3329.2460	0.0000	1.0000
JC69G	12	-1666.2779	3357.1899	0.0000	1.0000
JC69IG	13	-1666.2488	3359.2389	0.0000	1.0000
JC69I	12	-1667.4706	3359.5754	0.0000	1.0000
JC69	11	-1681.5771	3385.6897	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1586.6676	3213.3352	0.4546	0.4546
GTRIG	21	-1586.6475	3215.2950	0.1706	0.6252
HKYG	16	-1591.7385	3215.4771	0.1558	0.7810
GTRI	20	-1587.9774	3215.9548	0.1227	0.9037
HKYIG	17	-1591.7165	3217.4329	0.0586	0.9623
HKYI	16	-1593.1576	3218.3152	0.0377	1.0000
GTR	19	-1601.9439	3241.8878	0.0000	1.0000
HKY	15	-1606.6336	3243.2672	0.0000	1.0000
K2PG	13	-1632.0550	3290.1101	0.0000	1.0000
K2PIG	14	-1631.9783	3291.9566	0.0000	1.0000
K2PI	13	-1633.1536	3292.3072	0.0000	1.0000
F81G	15	-1633.3226	3296.6452	0.0000	1.0000
SYMG	17	-1631.6178	3297.2355	0.0000	1.0000
F81IG	16	-1633.3086	3298.6172	0.0000	1.0000
SYMIG	18	-1631.5767	3299.1533	0.0000	1.0000
SYMI	17	-1632.6780	3299.3559	0.0000	1.0000
F81I	15	-1634.9604	3299.9208	0.0000	1.0000
K2P	12	-1648.4978	3320.9956	0.0000	1.0000
F81	14	-1647.3510	3322.7019	0.0000	1.0000
SYM	16	-1648.0656	3328.1313	0.0000	1.0000
JC69G	12	-1666.2779	3356.5557	0.0000	1.0000
JC69IG	13	-1666.2488	3358.4975	0.0000	1.0000
JC69I	12	-1667.4706	3358.9412	0.0000	1.0000
JC69	11	-1681.5771	3385.1542	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1591.7385	3283.0700	0.7762	0.7762
HKYI	16	-1593.1576	3285.9082	0.1878	0.9640
HKYIG	17	-1591.7165	3289.2504	0.0353	0.9994
GTRG	20	-1586.6676	3297.8264	0.0005	0.9998
GTRI	20	-1587.9774	3300.4459	0.0001	1.0000
GTRIG	21	-1586.6475	3304.0107	0.0000	1.0000
HKY	15	-1606.6336	3306.6356	0.0000	1.0000
GTR	19	-1601.9439	3322.1545	0.0000	1.0000
K2PG	13	-1632.0550	3345.0293	0.0000	1.0000
K2PI	13	-1633.1536	3347.2265	0.0000	1.0000
K2PIG	14	-1631.9783	3351.1004	0.0000	1.0000
F81G	15	-1633.3226	3360.0136	0.0000	1.0000
F81I	15	-1634.9604	3363.2892	0.0000	1.0000
F81IG	16	-1633.3086	3366.2102	0.0000	1.0000
SYMG	17	-1631.6178	3369.0530	0.0000	1.0000
SYMI	17	-1632.6780	3371.1734	0.0000	1.0000
K2P	12	-1648.4978	3371.6903	0.0000	1.0000
SYMIG	18	-1631.5767	3375.1954	0.0000	1.0000
F81	14	-1647.3510	3381.8458	0.0000	1.0000
SYM	16	-1648.0656	3395.7242	0.0000	1.0000
JC69G	12	-1666.2779	3407.2504	0.0000	1.0000
JC69I	12	-1667.4706	3409.6359	0.0000	1.0000
JC69IG	13	-1666.2488	3413.4168	0.0000	1.0000
JC69	11	-1681.5771	3431.6243	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 683-tomodel-phylip/chr1_32461.phylip = 1 - 505;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32461.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 683-tomodel-phylip/chr1_32461.phylip = 1 - 505;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32461.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 683-tomodel-phylip/chr1_32461.phylip = 1 - 505;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32461.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

