
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:05 08/01/2011

Input data from file "261-aggressive-phylip/9201_hg18.phylip" (29 taxa, 395 characters)

Minimum AIC  model: SYMIG
Minimum AICc model: SYMIG
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
SYMIG	62	-1169.6378	2486.8057	0.3977	0.3977
SYMG	61	-1171.3639	2487.4426	0.2892	0.6869
SYMI	61	-1172.4849	2489.6845	0.0943	0.7811
K2PIG	58	-1177.2063	2490.7817	0.0545	0.8356
K2PG	57	-1178.5821	2490.7843	0.0544	0.8900
GTRIG	65	-1167.9128	2491.9047	0.0311	0.9211
GTRI	64	-1169.5529	2492.3179	0.0253	0.9463
K2PI	57	-1179.5258	2492.6719	0.0212	0.9675
GTRG	64	-1169.7357	2492.6835	0.0210	0.9885
HKYI	60	-1177.0396	2495.9953	0.0040	0.9926
HKYG	60	-1177.1110	2496.1382	0.0037	0.9963
HKYIG	61	-1175.7214	2496.1575	0.0037	1.0000
SYM	60	-1187.2990	2516.5141	0.0000	1.0000
K2P	56	-1194.3339	2519.5554	0.0000	1.0000
GTR	63	-1186.0407	2522.4440	0.0000	1.0000
HKY	59	-1193.0203	2525.1749	0.0000	1.0000
JC69G	56	-1198.0708	2527.0292	0.0000	1.0000
JC69I	56	-1198.9640	2528.8155	0.0000	1.0000
JC69IG	57	-1197.7410	2529.1021	0.0000	1.0000
F81G	59	-1196.4801	2532.0945	0.0000	1.0000
F81IG	60	-1196.1436	2534.2033	0.0000	1.0000
F81I	59	-1197.5883	2534.3110	0.0000	1.0000
JC69	55	-1212.7995	2553.7701	0.0000	1.0000
F81	58	-1211.3710	2559.1111	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMIG	62	-1169.6378	2463.2755	0.4423	0.4423
SYMG	61	-1171.3639	2464.7279	0.2140	0.6563
GTRIG	65	-1167.9128	2465.8256	0.1236	0.7799
SYMI	61	-1172.4849	2466.9697	0.0697	0.8496
GTRI	64	-1169.5529	2467.1058	0.0652	0.9148
GTRG	64	-1169.7357	2467.4714	0.0543	0.9690
K2PIG	58	-1177.2063	2470.4127	0.0125	0.9815
K2PG	57	-1178.5821	2471.1642	0.0086	0.9901
K2PI	57	-1179.5258	2473.0517	0.0033	0.9934
HKYIG	61	-1175.7214	2473.4427	0.0027	0.9961
HKYI	60	-1177.0396	2474.0791	0.0020	0.9981
HKYG	60	-1177.1110	2474.2221	0.0019	1.0000
SYM	60	-1187.2990	2494.5979	0.0000	1.0000
GTR	63	-1186.0407	2498.0815	0.0000	1.0000
K2P	56	-1194.3339	2500.6678	0.0000	1.0000
HKY	59	-1193.0203	2504.0406	0.0000	1.0000
JC69G	56	-1198.0708	2508.1416	0.0000	1.0000
JC69IG	57	-1197.7410	2509.4819	0.0000	1.0000
JC69I	56	-1198.9640	2509.9279	0.0000	1.0000
F81G	59	-1196.4801	2510.9602	0.0000	1.0000
F81IG	60	-1196.1436	2512.2871	0.0000	1.0000
F81I	59	-1197.5883	2513.1767	0.0000	1.0000
JC69	55	-1212.7995	2535.5990	0.0000	1.0000
F81	58	-1211.3710	2538.7421	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1178.5821	2697.9606	0.6236	0.6236
K2PI	57	-1179.5258	2699.8482	0.2427	0.8662
K2PIG	58	-1177.2063	2701.1880	0.1242	0.9904
SYMG	61	-1171.3639	2707.4399	0.0055	0.9959
SYMI	61	-1172.4849	2709.6818	0.0018	0.9977
SYMIG	62	-1169.6378	2709.9665	0.0015	0.9992
HKYI	60	-1177.0396	2712.8123	0.0004	0.9996
HKYG	60	-1177.1110	2712.9552	0.0003	0.9999
HKYIG	61	-1175.7214	2716.1548	0.0001	1.0000
GTRI	64	-1169.5529	2721.7545	0.0000	1.0000
GTRG	64	-1169.7357	2722.1200	0.0000	1.0000
K2P	56	-1194.3339	2723.4854	0.0000	1.0000
GTRIG	65	-1167.9128	2724.4532	0.0000	1.0000
JC69G	56	-1198.0708	2730.9592	0.0000	1.0000
JC69I	56	-1198.9640	2732.7455	0.0000	1.0000
SYM	60	-1187.2990	2733.3311	0.0000	1.0000
JC69IG	57	-1197.7410	2736.2784	0.0000	1.0000
HKY	59	-1193.0203	2738.7949	0.0000	1.0000
F81G	59	-1196.4801	2745.7145	0.0000	1.0000
F81I	59	-1197.5883	2747.9310	0.0000	1.0000
GTR	63	-1186.0407	2748.7513	0.0000	1.0000
F81IG	60	-1196.1436	2751.0202	0.0000	1.0000
JC69	55	-1212.7995	2754.4378	0.0000	1.0000
F81	58	-1211.3710	2769.5174	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9201_hg18.phylip = 1 - 395;
 Partition Dummy = 1:261-aggressive-phylip/9201_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9201_hg18.phylip = 1 - 395;
 Partition Dummy = 1:261-aggressive-phylip/9201_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9201_hg18.phylip = 1 - 395;
 Partition Dummy = 1:261-aggressive-phylip/9201_hg18.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

