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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1077.7171	2314.9495	0.7427	0.7427
GTRI	64	-1079.4865	2318.4882	0.1266	0.8693
GTRIG	65	-1078.0148	2318.6532	0.1166	0.9859
HKYG	60	-1088.3885	2324.0904	0.0077	0.9936
HKYI	60	-1088.9148	2325.1430	0.0045	0.9981
HKYIG	61	-1088.3444	2327.0184	0.0018	0.9999
SYMG	61	-1092.5200	2335.3695	0.0000	1.0000
SYMI	61	-1093.1841	2336.6978	0.0000	1.0000
K2PG	57	-1099.1894	2336.7773	0.0000	1.0000
K2PI	57	-1099.8090	2338.0165	0.0000	1.0000
SYMIG	62	-1092.4272	2338.2229	0.0000	1.0000
K2PIG	58	-1099.0739	2339.4959	0.0000	1.0000
GTR	63	-1093.8623	2344.1547	0.0000	1.0000
HKY	59	-1107.2300	2358.7797	0.0000	1.0000
SYM	60	-1109.7528	2366.8189	0.0000	1.0000
K2P	56	-1118.8830	2373.2365	0.0000	1.0000
F81G	59	-1114.8747	2374.0692	0.0000	1.0000
F81I	59	-1115.7741	2375.8679	0.0000	1.0000
F81IG	60	-1114.6253	2376.5641	0.0000	1.0000
JC69G	56	-1124.1301	2383.7308	0.0000	1.0000
JC69I	56	-1125.0905	2385.6517	0.0000	1.0000
JC69IG	57	-1123.8137	2386.0259	0.0000	1.0000
F81	58	-1134.2728	2409.8938	0.0000	1.0000
JC69	55	-1143.9533	2420.4707	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1077.7171	2283.4343	0.6915	0.6915
GTRIG	65	-1078.0148	2286.0296	0.1889	0.8804
GTRI	64	-1079.4865	2286.9730	0.1179	0.9983
HKYG	60	-1088.3885	2296.7769	0.0009	0.9991
HKYI	60	-1088.9148	2297.8296	0.0005	0.9997
HKYIG	61	-1088.3444	2298.6889	0.0003	1.0000
SYMG	61	-1092.5200	2307.0399	0.0000	1.0000
SYMI	61	-1093.1841	2308.3682	0.0000	1.0000
SYMIG	62	-1092.4272	2308.8544	0.0000	1.0000
K2PG	57	-1099.1894	2312.3788	0.0000	1.0000
K2PI	57	-1099.8090	2313.6179	0.0000	1.0000
GTR	63	-1093.8623	2313.7245	0.0000	1.0000
K2PIG	58	-1099.0739	2314.1477	0.0000	1.0000
HKY	59	-1107.2300	2332.4600	0.0000	1.0000
SYM	60	-1109.7528	2339.5055	0.0000	1.0000
F81G	59	-1114.8747	2347.7495	0.0000	1.0000
F81IG	60	-1114.6253	2349.2506	0.0000	1.0000
F81I	59	-1115.7741	2349.5482	0.0000	1.0000
K2P	56	-1118.8830	2349.7660	0.0000	1.0000
JC69G	56	-1124.1301	2360.2602	0.0000	1.0000
JC69IG	57	-1123.8137	2361.6274	0.0000	1.0000
JC69I	56	-1125.0905	2362.1811	0.0000	1.0000
F81	58	-1134.2728	2384.5457	0.0000	1.0000
JC69	55	-1143.9533	2397.9066	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1088.3885	2524.5404	0.4298	0.4298
HKYI	60	-1088.9148	2525.5930	0.2539	0.6837
GTRG	64	-1077.7171	2526.3820	0.1711	0.8548
K2PG	57	-1099.1894	2528.7541	0.0523	0.9071
GTRI	64	-1079.4865	2529.9207	0.0292	0.9363
K2PI	57	-1099.8090	2529.9932	0.0281	0.9644
HKYIG	61	-1088.3444	2530.2484	0.0248	0.9892
GTRIG	65	-1078.0148	2532.7734	0.0070	0.9962
K2PIG	58	-1099.0739	2534.3191	0.0032	0.9994
SYMG	61	-1092.5200	2538.5995	0.0004	0.9998
SYMI	61	-1093.1841	2539.9278	0.0002	1.0000
SYMIG	62	-1092.4272	2544.2100	0.0000	1.0000
GTR	63	-1093.8623	2552.8762	0.0000	1.0000
HKY	59	-1107.2300	2556.4274	0.0000	1.0000
K2P	56	-1118.8830	2562.3452	0.0000	1.0000
SYM	60	-1109.7528	2567.2690	0.0000	1.0000
F81G	59	-1114.8747	2571.7169	0.0000	1.0000
JC69G	56	-1124.1301	2572.8395	0.0000	1.0000
F81I	59	-1115.7741	2573.5156	0.0000	1.0000
JC69I	56	-1125.0905	2574.7603	0.0000	1.0000
F81IG	60	-1114.6253	2577.0141	0.0000	1.0000
JC69IG	57	-1123.8137	2578.0027	0.0000	1.0000
F81	58	-1134.2728	2604.7170	0.0000	1.0000
JC69	55	-1143.9533	2606.6898	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 261-aggressive-phylip/7803_hg18.phylip = 1 - 329;
 Partition Dummy = 1:261-aggressive-phylip/7803_hg18.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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7803_hg18.phylip = 1 - 329;
 Partition Dummy = 1:261-aggressive-phylip/7803_hg18.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 BIC model (HKYG)]

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

