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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
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
HKYG	60	-1040.5383	2227.4075	0.4085	0.4085
GTRG	64	-1035.3881	2229.1412	0.1717	0.5802
HKYI	60	-1041.5994	2229.5298	0.1414	0.7215
HKYIG	61	-1040.2534	2229.8137	0.1227	0.8442
GTRI	64	-1036.3079	2230.9808	0.0684	0.9126
GTRIG	65	-1035.1517	2231.7319	0.0470	0.9596
K2PG	57	-1047.9070	2233.3442	0.0210	0.9806
K2PI	57	-1049.0135	2235.5573	0.0069	0.9876
K2PIG	58	-1047.6659	2235.7747	0.0062	0.9938
SYMG	61	-1043.7365	2236.7798	0.0038	0.9976
SYMI	61	-1044.7087	2238.7243	0.0014	0.9990
SYMIG	62	-1043.5428	2239.3900	0.0010	1.0000
HKY	59	-1058.7100	2260.7963	0.0000	1.0000
K2P	56	-1065.0545	2264.7472	0.0000	1.0000
GTR	63	-1054.9302	2265.1840	0.0000	1.0000
F81G	59	-1061.9230	2267.2223	0.0000	1.0000
SYM	60	-1061.1900	2268.7110	0.0000	1.0000
F81I	59	-1062.9798	2269.3359	0.0000	1.0000
F81IG	60	-1061.7803	2269.8916	0.0000	1.0000
JC69G	56	-1068.6579	2271.9542	0.0000	1.0000
JC69I	56	-1069.5988	2273.8358	0.0000	1.0000
JC69IG	57	-1068.4748	2274.4797	0.0000	1.0000
F81	58	-1078.3312	2297.1052	0.0000	1.0000
JC69	55	-1084.6626	2301.0919	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1035.3881	2198.7762	0.4063	0.4063
GTRIG	65	-1035.1517	2200.3033	0.1893	0.5956
GTRI	64	-1036.3079	2200.6158	0.1619	0.7575
HKYG	60	-1040.5383	2201.0765	0.1286	0.8861
HKYIG	61	-1040.2534	2202.5068	0.0629	0.9490
HKYI	60	-1041.5994	2203.1988	0.0445	0.9935
SYMG	61	-1043.7365	2209.4729	0.0019	0.9955
K2PG	57	-1047.9070	2209.8140	0.0016	0.9971
SYMIG	62	-1043.5428	2211.0856	0.0009	0.9980
K2PIG	58	-1047.6659	2211.3318	0.0008	0.9987
SYMI	61	-1044.7087	2211.4175	0.0007	0.9995
K2PI	57	-1049.0135	2212.0271	0.0005	1.0000
HKY	59	-1058.7100	2235.4200	0.0000	1.0000
GTR	63	-1054.9302	2235.8604	0.0000	1.0000
F81G	59	-1061.9230	2241.8459	0.0000	1.0000
K2P	56	-1065.0545	2242.1089	0.0000	1.0000
SYM	60	-1061.1900	2242.3801	0.0000	1.0000
F81IG	60	-1061.7803	2243.5606	0.0000	1.0000
F81I	59	-1062.9798	2243.9596	0.0000	1.0000
JC69G	56	-1068.6579	2249.3159	0.0000	1.0000
JC69IG	57	-1068.4748	2250.9495	0.0000	1.0000
JC69I	56	-1069.5988	2251.1975	0.0000	1.0000
F81	58	-1078.3312	2272.6623	0.0000	1.0000
JC69	55	-1084.6626	2279.3252	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1047.9070	2427.8960	0.5677	0.5677
K2PI	57	-1049.0135	2430.1091	0.1877	0.7554
HKYG	60	-1040.5383	2430.6365	0.1442	0.8997
HKYI	60	-1041.5994	2432.7588	0.0499	0.9496
K2PIG	58	-1047.6659	2433.2398	0.0392	0.9888
HKYIG	61	-1040.2534	2435.8928	0.0104	0.9992
SYMG	61	-1043.7365	2442.8589	0.0003	0.9995
GTRG	64	-1035.3881	2443.6402	0.0002	0.9998
SYMI	61	-1044.7087	2444.8035	0.0001	0.9999
GTRI	64	-1036.3079	2445.4798	0.0001	1.0000
SYMIG	62	-1043.5428	2448.2976	0.0000	1.0000
GTRIG	65	-1035.1517	2448.9933	0.0000	1.0000
K2P	56	-1065.0545	2456.3649	0.0000	1.0000
HKY	59	-1058.7100	2461.1540	0.0000	1.0000
JC69G	56	-1068.6579	2463.5719	0.0000	1.0000
JC69I	56	-1069.5988	2465.4535	0.0000	1.0000
F81G	59	-1061.9230	2467.5799	0.0000	1.0000
JC69IG	57	-1068.4748	2469.0315	0.0000	1.0000
F81I	59	-1062.9798	2469.6936	0.0000	1.0000
SYM	60	-1061.1900	2471.9401	0.0000	1.0000
F81IG	60	-1061.7803	2473.1206	0.0000	1.0000
GTR	63	-1054.9302	2476.8984	0.0000	1.0000
JC69	55	-1084.6626	2489.7552	0.0000	1.0000
F81	58	-1078.3312	2494.5703	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 183-phylip/chr9_5218.phylip = 1 - 339;
 Partition Dummy = 1:183-phylip/chr9_5218.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr9_5218.phylip = 1 - 339;
 Partition Dummy = 1:183-phylip/chr9_5218.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr9_5218.phylip = 1 - 339;
 Partition Dummy = 1:183-phylip/chr9_5218.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

