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

Input data from file "683-tomodel-phylip/chr7_6238.phylip" (7 taxa, 603 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	-1649.1766	3339.7966	0.5627	0.5627
GTRI	20	-1650.0907	3341.6246	0.2256	0.7883
GTRIG	21	-1649.1261	3341.8425	0.2023	0.9906
HKYG	16	-1658.0995	3349.1273	0.0053	0.9959
HKYI	16	-1659.0230	3350.9742	0.0021	0.9980
HKYIG	17	-1658.0494	3351.1450	0.0019	0.9999
GTR	19	-1659.4837	3358.2710	0.0001	1.0000
HKY	15	-1667.9336	3366.6849	0.0000	1.0000
SYMG	17	-1668.7258	3372.4977	0.0000	1.0000
SYMI	17	-1669.4764	3373.9990	0.0000	1.0000
SYMIG	18	-1668.7134	3374.5980	0.0000	1.0000
K2PG	13	-1676.7305	3380.0789	0.0000	1.0000
K2PI	13	-1677.5593	3381.7366	0.0000	1.0000
K2PIG	14	-1676.7026	3382.1194	0.0000	1.0000
SYM	16	-1676.8799	3386.6880	0.0000	1.0000
K2P	12	-1685.3957	3395.3203	0.0000	1.0000
F81G	15	-1691.5191	3413.8559	0.0000	1.0000
F81I	15	-1692.0670	3414.9517	0.0000	1.0000
F81IG	16	-1691.5115	3415.9513	0.0000	1.0000
F81	14	-1699.6744	3428.0630	0.0000	1.0000
JC69G	12	-1706.9865	3438.5019	0.0000	1.0000
JC69I	12	-1707.5289	3439.5867	0.0000	1.0000
JC69IG	13	-1706.9822	3440.5824	0.0000	1.0000
JC69	11	-1714.6558	3451.7583	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1649.1766	3338.3533	0.5553	0.5553
GTRI	20	-1650.0907	3340.1813	0.2226	0.7779
GTRIG	21	-1649.1261	3340.2521	0.2149	0.9927
HKYG	16	-1658.0995	3348.1990	0.0040	0.9968
HKYI	16	-1659.0230	3350.0459	0.0016	0.9984
HKYIG	17	-1658.0494	3350.0988	0.0016	0.9999
GTR	19	-1659.4837	3356.9674	0.0001	1.0000
HKY	15	-1667.9336	3365.8671	0.0000	1.0000
SYMG	17	-1668.7258	3371.4515	0.0000	1.0000
SYMI	17	-1669.4764	3372.9528	0.0000	1.0000
SYMIG	18	-1668.7134	3373.4268	0.0000	1.0000
K2PG	13	-1676.7305	3379.4609	0.0000	1.0000
K2PI	13	-1677.5593	3381.1186	0.0000	1.0000
K2PIG	14	-1676.7026	3381.4051	0.0000	1.0000
SYM	16	-1676.8799	3385.7597	0.0000	1.0000
K2P	12	-1685.3957	3394.7915	0.0000	1.0000
F81G	15	-1691.5191	3413.0382	0.0000	1.0000
F81I	15	-1692.0670	3414.1340	0.0000	1.0000
F81IG	16	-1691.5115	3415.0230	0.0000	1.0000
F81	14	-1699.6744	3427.3487	0.0000	1.0000
JC69G	12	-1706.9865	3437.9731	0.0000	1.0000
JC69I	12	-1707.5289	3439.0578	0.0000	1.0000
JC69IG	13	-1706.9822	3439.9644	0.0000	1.0000
JC69	11	-1714.6558	3451.3116	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1658.0995	3418.6297	0.6797	0.6797
HKYI	16	-1659.0230	3420.4766	0.2699	0.9496
HKYIG	17	-1658.0494	3424.9314	0.0291	0.9787
GTRG	20	-1649.1766	3426.3916	0.0140	0.9927
GTRI	20	-1650.0907	3428.2196	0.0056	0.9984
HKY	15	-1667.9336	3431.8959	0.0009	0.9993
GTRIG	21	-1649.1261	3432.6924	0.0006	0.9999
K2PG	13	-1676.7305	3436.6859	0.0001	0.9999
K2PI	13	-1677.5593	3438.3435	0.0000	1.0000
GTR	19	-1659.4837	3440.6038	0.0000	1.0000
K2PIG	14	-1676.7026	3443.0320	0.0000	1.0000
SYMG	17	-1668.7258	3446.2841	0.0000	1.0000
K2P	12	-1685.3957	3447.6145	0.0000	1.0000
SYMI	17	-1669.4764	3447.7854	0.0000	1.0000
SYMIG	18	-1668.7134	3452.6613	0.0000	1.0000
SYM	16	-1676.8799	3456.1904	0.0000	1.0000
F81G	15	-1691.5191	3479.0670	0.0000	1.0000
F81I	15	-1692.0670	3480.1628	0.0000	1.0000
F81IG	16	-1691.5115	3485.4537	0.0000	1.0000
F81	14	-1699.6744	3488.9755	0.0000	1.0000
JC69G	12	-1706.9865	3490.7961	0.0000	1.0000
JC69I	12	-1707.5289	3491.8808	0.0000	1.0000
JC69IG	13	-1706.9822	3497.1894	0.0000	1.0000
JC69	11	-1714.6558	3499.7326	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/chr7_6238.phylip = 1 - 603;
 Partition Dummy = 1:683-tomodel-phylip/chr7_6238.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/chr7_6238.phylip = 1 - 603;
 Partition Dummy = 1:683-tomodel-phylip/chr7_6238.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/chr7_6238.phylip = 1 - 603;
 Partition Dummy = 1:683-tomodel-phylip/chr7_6238.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

