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

Input data from file "183-phylip/chr13_710.phylip" (29 taxa, 386 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	60	-1066.5604	2275.6438	0.4094	0.4094
HKYI	60	-1066.5886	2275.7002	0.3980	0.8074
HKYIG	61	-1066.4931	2278.3320	0.1068	0.9142
GTRI	64	-1063.2250	2280.3691	0.0386	0.9528
GTRG	64	-1063.2342	2280.3874	0.0382	0.9910
GTRIG	65	-1063.2337	2283.2800	0.0090	1.0000
K2PG	57	-1080.8607	2295.8800	0.0000	1.0000
K2PI	57	-1081.3251	2296.8088	0.0000	1.0000
K2PIG	58	-1080.3371	2297.6039	0.0000	1.0000
SYMG	61	-1078.7693	2302.8843	0.0000	1.0000
SYMI	61	-1078.7870	2302.9197	0.0000	1.0000
SYMIG	62	-1078.5458	2305.2773	0.0000	1.0000
F81G	59	-1093.3474	2326.4126	0.0000	1.0000
F81I	59	-1093.3685	2326.4548	0.0000	1.0000
F81IG	60	-1093.2329	2328.9888	0.0000	1.0000
HKY	59	-1096.7122	2333.1423	0.0000	1.0000
GTR	63	-1092.3214	2335.6863	0.0000	1.0000
JC69G	56	-1106.6573	2344.7190	0.0000	1.0000
JC69I	56	-1107.1276	2345.6595	0.0000	1.0000
JC69IG	57	-1106.8385	2347.8356	0.0000	1.0000
K2P	56	-1112.9672	2357.3387	0.0000	1.0000
SYM	60	-1109.9058	2362.3348	0.0000	1.0000
F81	58	-1123.7651	2384.4599	0.0000	1.0000
JC69	55	-1138.0277	2404.7220	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1066.5604	2253.1207	0.2795	0.2795
HKYI	60	-1066.5886	2253.1771	0.2718	0.5513
GTRI	64	-1063.2250	2254.4501	0.1438	0.6951
GTRG	64	-1063.2342	2254.4684	0.1425	0.8376
HKYIG	61	-1066.4931	2254.9863	0.1100	0.9475
GTRIG	65	-1063.2337	2256.4675	0.0524	1.0000
K2PG	57	-1080.8607	2275.7215	0.0000	1.0000
K2PI	57	-1081.3251	2276.6502	0.0000	1.0000
K2PIG	58	-1080.3371	2276.6743	0.0000	1.0000
SYMG	61	-1078.7693	2279.5386	0.0000	1.0000
SYMI	61	-1078.7870	2279.5740	0.0000	1.0000
SYMIG	62	-1078.5458	2281.0916	0.0000	1.0000
F81G	59	-1093.3474	2304.6948	0.0000	1.0000
F81I	59	-1093.3685	2304.7370	0.0000	1.0000
F81IG	60	-1093.2329	2306.4657	0.0000	1.0000
GTR	63	-1092.3214	2310.6428	0.0000	1.0000
HKY	59	-1096.7122	2311.4245	0.0000	1.0000
JC69G	56	-1106.6573	2325.3147	0.0000	1.0000
JC69I	56	-1107.1276	2326.2552	0.0000	1.0000
JC69IG	57	-1106.8385	2327.6770	0.0000	1.0000
K2P	56	-1112.9672	2337.9345	0.0000	1.0000
SYM	60	-1109.9058	2339.8117	0.0000	1.0000
F81	58	-1123.7651	2363.5302	0.0000	1.0000
JC69	55	-1138.0277	2386.0553	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1066.5604	2490.4710	0.4914	0.4914
HKYI	60	-1066.5886	2490.5274	0.4777	0.9691
HKYIG	61	-1066.4931	2496.2924	0.0268	0.9959
K2PG	57	-1080.8607	2501.2042	0.0023	0.9982
K2PI	57	-1081.3251	2502.1330	0.0014	0.9996
K2PIG	58	-1080.3371	2506.1128	0.0002	0.9998
GTRI	64	-1063.2250	2507.6237	0.0001	0.9999
GTRG	64	-1063.2342	2507.6420	0.0001	1.0000
GTRIG	65	-1063.2337	2513.5969	0.0000	1.0000
SYMG	61	-1078.7693	2520.8447	0.0000	1.0000
SYMI	61	-1078.7870	2520.8801	0.0000	1.0000
SYMIG	62	-1078.5458	2526.3535	0.0000	1.0000
F81G	59	-1093.3474	2538.0892	0.0000	1.0000
F81I	59	-1093.3685	2538.1314	0.0000	1.0000
F81IG	60	-1093.2329	2543.8160	0.0000	1.0000
HKY	59	-1096.7122	2544.8189	0.0000	1.0000
JC69G	56	-1106.6573	2546.8416	0.0000	1.0000
JC69I	56	-1107.1276	2547.7821	0.0000	1.0000
JC69IG	57	-1106.8385	2553.1598	0.0000	1.0000
K2P	56	-1112.9672	2559.4614	0.0000	1.0000
GTR	63	-1092.3214	2559.8605	0.0000	1.0000
SYM	60	-1109.9058	2577.1619	0.0000	1.0000
F81	58	-1123.7651	2592.9688	0.0000	1.0000
JC69	55	-1138.0277	2603.6264	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 183-phylip/chr13_710.phylip = 1 - 386;
 Partition Dummy = 1:183-phylip/chr13_710.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 183-phylip/chr13_710.phylip = 1 - 386;
 Partition Dummy = 1:183-phylip/chr13_710.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 183-phylip/chr13_710.phylip = 1 - 386;
 Partition Dummy = 1:183-phylip/chr13_710.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

