
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:49 02/17/2011

Input data from file "183-phylip/chr3_6115.phylip" (29 taxa, 333 characters)

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	57	-1332.4973	2803.0383	0.8133	0.8133
HKYG	60	-1330.1843	2807.2805	0.0975	0.9108
HKYIG	61	-1329.0881	2808.0876	0.0651	0.9759
HKYI	60	-1332.2355	2811.3828	0.0125	0.9885
SYMG	61	-1331.3815	2812.6743	0.0066	0.9950
SYMIG	62	-1330.4736	2813.8806	0.0036	0.9986
GTRG	64	-1328.9854	2817.0156	0.0008	0.9994
K2PIG	58	-1338.4000	2817.7780	0.0005	0.9999
K2PI	57	-1341.4374	2820.9185	0.0001	1.0000
SYMI	61	-1339.5596	2829.0306	0.0000	1.0000
GTRIG	65	-1335.2006	2832.5361	0.0000	1.0000
GTRI	64	-1336.8983	2832.8413	0.0000	1.0000
JC69G	56	-1370.9732	2877.0769	0.0000	1.0000
F81G	59	-1369.4508	2882.8357	0.0000	1.0000
F81I	59	-1370.1257	2884.1855	0.0000	1.0000
JC69IG	57	-1374.2233	2886.4903	0.0000	1.0000
JC69I	56	-1375.7962	2886.7228	0.0000	1.0000
F81IG	60	-1372.6041	2892.1200	0.0000	1.0000
K2P	56	-1387.6766	2910.4836	0.0000	1.0000
HKY	59	-1385.6073	2915.1487	0.0000	1.0000
SYM	60	-1386.7797	2920.4712	0.0000	1.0000
GTR	63	-1384.6391	2925.2560	0.0000	1.0000
JC69	55	-1419.4787	2971.1957	0.0000	1.0000
F81	58	-1417.8437	2976.6655	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	57	-1332.4973	2778.9946	0.4424	0.4424
HKYIG	61	-1329.0881	2780.1762	0.2451	0.6875
HKYG	60	-1330.1843	2780.3687	0.2226	0.9100
HKYI	60	-1332.2355	2784.4711	0.0286	0.9387
SYMG	61	-1331.3815	2784.7629	0.0247	0.9634
SYMIG	62	-1330.4736	2784.9473	0.0226	0.9860
GTRG	64	-1328.9854	2785.9708	0.0135	0.9995
K2PIG	58	-1338.4000	2792.7999	0.0004	0.9999
K2PI	57	-1341.4374	2796.8748	0.0001	1.0000
GTRIG	65	-1335.2006	2800.4013	0.0000	1.0000
SYMI	61	-1339.5596	2801.1191	0.0000	1.0000
GTRI	64	-1336.8983	2801.7965	0.0000	1.0000
JC69G	56	-1370.9732	2853.9464	0.0000	1.0000
F81G	59	-1369.4508	2856.9016	0.0000	1.0000
F81I	59	-1370.1257	2858.2514	0.0000	1.0000
JC69IG	57	-1374.2233	2862.4466	0.0000	1.0000
JC69I	56	-1375.7962	2863.5923	0.0000	1.0000
F81IG	60	-1372.6041	2865.2083	0.0000	1.0000
K2P	56	-1387.6766	2887.3532	0.0000	1.0000
HKY	59	-1385.6073	2889.2147	0.0000	1.0000
SYM	60	-1386.7797	2893.5594	0.0000	1.0000
GTR	63	-1384.6391	2895.2783	0.0000	1.0000
JC69	55	-1419.4787	2948.9575	0.0000	1.0000
F81	58	-1417.8437	2951.6874	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1332.4973	2996.0587	0.9975	0.9975
HKYG	60	-1330.1843	3008.8572	0.0017	0.9992
HKYIG	61	-1329.0881	3012.4729	0.0003	0.9995
HKYI	60	-1332.2355	3012.9596	0.0002	0.9997
K2PIG	58	-1338.4000	3013.6722	0.0001	0.9998
K2PI	57	-1341.4374	3013.9390	0.0001	1.0000
SYMG	61	-1331.3815	3017.0596	0.0000	1.0000
SYMIG	62	-1330.4736	3021.0521	0.0000	1.0000
GTRG	64	-1328.9854	3029.6920	0.0000	1.0000
SYMI	61	-1339.5596	3033.4158	0.0000	1.0000
GTRI	64	-1336.8983	3045.5177	0.0000	1.0000
GTRIG	65	-1335.2006	3047.9305	0.0000	1.0000
JC69G	56	-1370.9732	3067.2024	0.0000	1.0000
JC69I	56	-1375.7962	3076.8483	0.0000	1.0000
JC69IG	57	-1374.2233	3079.5108	0.0000	1.0000
F81G	59	-1369.4508	3081.5820	0.0000	1.0000
F81I	59	-1370.1257	3082.9318	0.0000	1.0000
F81IG	60	-1372.6041	3093.6968	0.0000	1.0000
K2P	56	-1387.6766	3100.6092	0.0000	1.0000
HKY	59	-1385.6073	3113.8951	0.0000	1.0000
SYM	60	-1386.7797	3122.0480	0.0000	1.0000
GTR	63	-1384.6391	3135.1913	0.0000	1.0000
JC69	55	-1419.4787	3158.4053	0.0000	1.0000
F81	58	-1417.8437	3172.5597	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_6115.phylip = 1 - 333;
 Partition Dummy = 1:183-phylip/chr3_6115.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_6115.phylip = 1 - 333;
 Partition Dummy = 1:183-phylip/chr3_6115.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_6115.phylip = 1 - 333;
 Partition Dummy = 1:183-phylip/chr3_6115.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

