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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1395.2152	2823.2132	0.3890	0.3890
HKYI	16	-1395.3050	2823.3927	0.3556	0.7447
HKYIG	17	-1395.1576	2825.1970	0.1443	0.8890
GTRG	20	-1393.2729	2827.7614	0.0400	0.9290
GTRI	20	-1393.3602	2827.9361	0.0367	0.9657
HKY	15	-1399.3701	2829.4300	0.0174	0.9830
GTRIG	21	-1393.2128	2829.7648	0.0147	0.9977
GTR	19	-1397.5041	2834.1064	0.0017	0.9994
F81I	15	-1403.6198	2837.9292	0.0002	0.9997
F81G	15	-1403.6869	2838.0635	0.0002	0.9999
F81IG	16	-1403.6870	2840.1567	0.0001	1.0000
F81	14	-1407.4973	2843.5972	0.0000	1.0000
K2PG	13	-1432.0967	2890.7149	0.0000	1.0000
K2PI	13	-1432.1567	2890.8349	0.0000	1.0000
K2PIG	14	-1432.0755	2892.7535	0.0000	1.0000
K2P	12	-1435.9638	2896.3740	0.0000	1.0000
SYMG	17	-1431.2679	2897.4177	0.0000	1.0000
SYMI	17	-1431.3068	2897.4954	0.0000	1.0000
SYMIG	18	-1431.2561	2899.4992	0.0000	1.0000
JC69I	12	-1438.4792	2901.4048	0.0000	1.0000
JC69G	12	-1438.5593	2901.5649	0.0000	1.0000
SYM	16	-1435.0967	2902.9762	0.0000	1.0000
JC69IG	13	-1438.5592	2903.6400	0.0000	1.0000
JC69	11	-1442.2529	2906.8830	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1395.2152	2822.4304	0.3777	0.3777
HKYI	16	-1395.3050	2822.6100	0.3453	0.7230
HKYIG	17	-1395.1576	2824.3152	0.1472	0.8701
GTRG	20	-1393.2729	2826.5458	0.0482	0.9184
GTRI	20	-1393.3602	2826.7205	0.0442	0.9626
GTRIG	21	-1393.2128	2828.4256	0.0188	0.9815
HKY	15	-1399.3701	2828.7403	0.0161	0.9976
GTR	19	-1397.5041	2833.0081	0.0019	0.9995
F81I	15	-1403.6198	2837.2396	0.0002	0.9997
F81G	15	-1403.6869	2837.3739	0.0002	0.9999
F81IG	16	-1403.6870	2839.3739	0.0001	1.0000
F81	14	-1407.4973	2842.9946	0.0000	1.0000
K2PG	13	-1432.0967	2890.1934	0.0000	1.0000
K2PI	13	-1432.1567	2890.3134	0.0000	1.0000
K2PIG	14	-1432.0755	2892.1509	0.0000	1.0000
K2P	12	-1435.9638	2895.9277	0.0000	1.0000
SYMG	17	-1431.2679	2896.5359	0.0000	1.0000
SYMI	17	-1431.3068	2896.6135	0.0000	1.0000
SYMIG	18	-1431.2561	2898.5122	0.0000	1.0000
JC69I	12	-1438.4792	2900.9585	0.0000	1.0000
JC69G	12	-1438.5593	2901.1186	0.0000	1.0000
SYM	16	-1435.0967	2902.1935	0.0000	1.0000
JC69IG	13	-1438.5592	2903.1185	0.0000	1.0000
JC69	11	-1442.2529	2906.5058	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1395.2152	2895.5197	0.4188	0.4188
HKYI	16	-1395.3050	2895.6992	0.3829	0.8017
HKY	15	-1399.3701	2897.2615	0.1753	0.9770
HKYIG	17	-1395.1576	2901.9725	0.0166	0.9937
F81I	15	-1403.6198	2905.7607	0.0025	0.9962
F81G	15	-1403.6869	2905.8950	0.0023	0.9985
F81	14	-1407.4973	2906.9477	0.0014	0.9999
F81IG	16	-1403.6870	2912.4632	0.0001	1.0000
GTRG	20	-1393.2729	2917.9074	0.0000	1.0000
GTRI	20	-1393.3602	2918.0821	0.0000	1.0000
GTR	19	-1397.5041	2919.8016	0.0000	1.0000
GTRIG	21	-1393.2128	2924.3553	0.0000	1.0000
K2PG	13	-1432.0967	2949.5785	0.0000	1.0000
K2PI	13	-1432.1567	2949.6984	0.0000	1.0000
K2P	12	-1435.9638	2950.7446	0.0000	1.0000
JC69I	12	-1438.4792	2955.7754	0.0000	1.0000
JC69G	12	-1438.5593	2955.9355	0.0000	1.0000
K2PIG	14	-1432.0755	2956.1040	0.0000	1.0000
JC69	11	-1442.2529	2956.7547	0.0000	1.0000
JC69IG	13	-1438.5592	2962.5035	0.0000	1.0000
SYMG	17	-1431.2679	2974.1932	0.0000	1.0000
SYMI	17	-1431.3068	2974.2708	0.0000	1.0000
SYM	16	-1435.0967	2975.2827	0.0000	1.0000
SYMIG	18	-1431.2561	2980.7376	0.0000	1.0000



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

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

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


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

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


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

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

