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

Input data from file "591-tomodel-phylip/chr5_1800.phylip" (7 taxa, 642 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	-1392.8599	2827.0724	0.4406	0.4406
GTRI	20	-1393.0359	2827.4244	0.3695	0.8100
GTRIG	21	-1392.7046	2828.8995	0.1767	0.9867
HKYG	16	-1401.4484	2835.7671	0.0057	0.9924
HKYI	16	-1401.6415	2836.1534	0.0047	0.9971
HKYIG	17	-1401.1032	2837.1872	0.0028	0.9999
GTR	19	-1402.8264	2844.8746	0.0001	1.0000
HKY	15	-1412.1797	2855.1262	0.0000	1.0000
F81G	15	-1413.2895	2857.3457	0.0000	1.0000
F81I	15	-1413.5029	2857.7726	0.0000	1.0000
F81IG	16	-1412.9012	2858.6728	0.0000	1.0000
SYMG	17	-1418.8488	2872.6783	0.0000	1.0000
SYMI	17	-1419.0602	2873.1013	0.0000	1.0000
SYMIG	18	-1418.6410	2874.3799	0.0000	1.0000
F81	14	-1423.5193	2875.7085	0.0000	1.0000
K2PG	13	-1425.3681	2877.3159	0.0000	1.0000
K2PI	13	-1425.5902	2877.7601	0.0000	1.0000
K2PIG	14	-1425.0122	2878.6943	0.0000	1.0000
SYM	16	-1428.1427	2889.1559	0.0000	1.0000
K2P	12	-1435.5406	2895.5773	0.0000	1.0000
JC69G	12	-1435.6634	2895.8229	0.0000	1.0000
JC69I	12	-1435.8775	2896.2511	0.0000	1.0000
JC69IG	13	-1435.3105	2897.2007	0.0000	1.0000
JC69	11	-1445.3998	2913.2186	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1392.8599	2825.7197	0.4363	0.4363
GTRI	20	-1393.0359	2826.0718	0.3658	0.8021
GTRIG	21	-1392.7046	2827.4091	0.1875	0.9895
HKYG	16	-1401.4484	2834.8967	0.0044	0.9940
HKYI	16	-1401.6415	2835.2830	0.0037	0.9976
HKYIG	17	-1401.1032	2836.2065	0.0023	0.9999
GTR	19	-1402.8264	2843.6527	0.0001	1.0000
HKY	15	-1412.1797	2854.3594	0.0000	1.0000
F81G	15	-1413.2895	2856.5789	0.0000	1.0000
F81I	15	-1413.5029	2857.0059	0.0000	1.0000
F81IG	16	-1412.9012	2857.8024	0.0000	1.0000
SYMG	17	-1418.8488	2871.6975	0.0000	1.0000
SYMI	17	-1419.0602	2872.1205	0.0000	1.0000
SYMIG	18	-1418.6410	2873.2819	0.0000	1.0000
F81	14	-1423.5193	2875.0386	0.0000	1.0000
K2PG	13	-1425.3681	2876.7363	0.0000	1.0000
K2PI	13	-1425.5902	2877.1804	0.0000	1.0000
K2PIG	14	-1425.0122	2878.0245	0.0000	1.0000
SYM	16	-1428.1427	2888.2855	0.0000	1.0000
K2P	12	-1435.5406	2895.0812	0.0000	1.0000
JC69G	12	-1435.6634	2895.3268	0.0000	1.0000
JC69I	12	-1435.8775	2895.7550	0.0000	1.0000
JC69IG	13	-1435.3105	2896.6211	0.0000	1.0000
JC69	11	-1445.3998	2912.7996	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1401.4484	2906.3302	0.5248	0.5248
HKYI	16	-1401.6415	2906.7165	0.4326	0.9574
HKYIG	17	-1401.1032	2912.1045	0.0292	0.9866
GTRG	20	-1392.8599	2915.0115	0.0068	0.9935
GTRI	20	-1393.0359	2915.3635	0.0057	0.9992
GTRIG	21	-1392.7046	2921.1655	0.0003	0.9995
HKY	15	-1412.1797	2921.3282	0.0003	0.9998
F81G	15	-1413.2895	2923.5477	0.0001	0.9999
F81I	15	-1413.5029	2923.9747	0.0001	1.0000
GTR	19	-1402.8264	2928.4799	0.0000	1.0000
F81IG	16	-1412.9012	2929.2358	0.0000	1.0000
K2PG	13	-1425.3681	2934.7759	0.0000	1.0000
K2PI	13	-1425.5902	2935.2201	0.0000	1.0000
F81	14	-1423.5193	2937.5428	0.0000	1.0000
K2PIG	14	-1425.0122	2940.5287	0.0000	1.0000
SYMG	17	-1418.8488	2947.5955	0.0000	1.0000
SYMI	17	-1419.0602	2948.0185	0.0000	1.0000
K2P	12	-1435.5406	2948.6563	0.0000	1.0000
JC69G	12	-1435.6634	2948.9019	0.0000	1.0000
JC69I	12	-1435.8775	2949.3301	0.0000	1.0000
SYMIG	18	-1418.6410	2953.6445	0.0000	1.0000
JC69IG	13	-1435.3105	2954.6607	0.0000	1.0000
SYM	16	-1428.1427	2959.7189	0.0000	1.0000
JC69	11	-1445.3998	2961.9100	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_1800.phylip = 1 - 642;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1800.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_1800.phylip = 1 - 642;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1800.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/chr5_1800.phylip = 1 - 642;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1800.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

