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

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

Minimum AIC  model: GTRG
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	-1125.9626	2393.4753	0.3344	0.3344
HKYG	60	-1121.4960	2393.8781	0.2734	0.6077
K2PI	57	-1126.9907	2395.5314	0.1196	0.7274
HKYI	60	-1122.6828	2396.2517	0.0834	0.8108
K2PIG	58	-1125.9252	2396.4863	0.0742	0.8850
HKYIG	61	-1121.4630	2396.9768	0.0581	0.9430
GTRG	64	-1117.3566	2398.4214	0.0282	0.9712
SYMG	61	-1123.0578	2400.1664	0.0118	0.9830
GTRI	64	-1118.3704	2400.4490	0.0102	0.9933
SYMI	61	-1123.8385	2401.7279	0.0054	0.9987
GTRIG	65	-1119.0628	2405.1084	0.0010	0.9996
SYMIG	62	-1124.9773	2407.1972	0.0004	1.0000
HKY	59	-1137.0262	2421.8003	0.0000	1.0000
K2P	56	-1141.8781	2422.2459	0.0000	1.0000
GTR	63	-1132.5873	2425.6362	0.0000	1.0000
SYM	60	-1138.3682	2427.6224	0.0000	1.0000
JC69G	56	-1146.6213	2431.7322	0.0000	1.0000
F81G	59	-1142.0108	2431.7696	0.0000	1.0000
F81I	59	-1143.2890	2434.3258	0.0000	1.0000
JC69I	56	-1147.9743	2434.4382	0.0000	1.0000
JC69IG	57	-1146.5316	2434.6132	0.0000	1.0000
F81IG	60	-1141.9354	2434.7568	0.0000	1.0000
JC69	55	-1161.0693	2457.5932	0.0000	1.0000
F81	58	-1156.6829	2458.0018	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1117.3566	2362.7132	0.2988	0.2988
HKYG	60	-1121.4960	2362.9920	0.2600	0.5588
GTRI	64	-1118.3704	2364.7408	0.1084	0.6672
HKYIG	61	-1121.4630	2364.9260	0.0988	0.7661
HKYI	60	-1122.6828	2365.3656	0.0793	0.8454
K2PG	57	-1125.9626	2365.9253	0.0600	0.9054
K2PIG	58	-1125.9252	2367.8503	0.0229	0.9283
K2PI	57	-1126.9907	2367.9814	0.0215	0.9497
SYMG	61	-1123.0578	2368.1156	0.0201	0.9698
GTRIG	65	-1119.0628	2368.1256	0.0200	0.9897
SYMI	61	-1123.8385	2369.6771	0.0092	0.9989
SYMIG	62	-1124.9773	2373.9547	0.0011	1.0000
GTR	63	-1132.5873	2391.1747	0.0000	1.0000
HKY	59	-1137.0262	2392.0524	0.0000	1.0000
K2P	56	-1141.8781	2395.7562	0.0000	1.0000
SYM	60	-1138.3682	2396.7363	0.0000	1.0000
F81G	59	-1142.0108	2402.0217	0.0000	1.0000
F81IG	60	-1141.9354	2403.8708	0.0000	1.0000
F81I	59	-1143.2890	2404.5779	0.0000	1.0000
JC69G	56	-1146.6213	2405.2426	0.0000	1.0000
JC69IG	57	-1146.5316	2407.0632	0.0000	1.0000
JC69I	56	-1147.9743	2407.9486	0.0000	1.0000
F81	58	-1156.6829	2429.3659	0.0000	1.0000
JC69	55	-1161.0693	2432.1387	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1125.9626	2576.6596	0.6938	0.6938
K2PI	57	-1126.9907	2578.7157	0.2482	0.9420
K2PIG	58	-1125.9252	2582.2818	0.0417	0.9837
HKYG	60	-1121.4960	2584.8176	0.0117	0.9955
HKYI	60	-1122.6828	2587.1912	0.0036	0.9991
HKYIG	61	-1121.4630	2590.4487	0.0007	0.9998
SYMG	61	-1123.0578	2593.6383	0.0001	0.9999
SYMI	61	-1123.8385	2595.1998	0.0001	1.0000
GTRG	64	-1117.3566	2599.3272	0.0000	1.0000
GTRI	64	-1118.3704	2601.3548	0.0000	1.0000
K2P	56	-1141.8781	2602.7935	0.0000	1.0000
SYMIG	62	-1124.9773	2603.1745	0.0000	1.0000
GTRIG	65	-1119.0628	2608.4367	0.0000	1.0000
HKY	59	-1137.0262	2610.1809	0.0000	1.0000
JC69G	56	-1146.6213	2612.2799	0.0000	1.0000
JC69I	56	-1147.9743	2614.9858	0.0000	1.0000
JC69IG	57	-1146.5316	2617.7976	0.0000	1.0000
SYM	60	-1138.3682	2618.5619	0.0000	1.0000
F81G	59	-1142.0108	2620.1502	0.0000	1.0000
F81I	59	-1143.2890	2622.7064	0.0000	1.0000
GTR	63	-1132.5873	2624.0916	0.0000	1.0000
F81IG	60	-1141.9354	2625.6964	0.0000	1.0000
JC69	55	-1161.0693	2635.4788	0.0000	1.0000
F81	58	-1156.6829	2643.7973	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/chr5_1649.phylip = 1 - 298;
 Partition Dummy = 1:183-phylip/chr5_1649.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 (K2PG)]

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

