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

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

Minimum AIC  model: GTRG
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	-1071.1363	2299.4301	0.7101	0.7101
HKYIG	61	-1071.0394	2302.6706	0.1405	0.8506
GTRG	64	-1066.5085	2304.1257	0.0679	0.9184
HKYI	60	-1073.5493	2304.2560	0.0636	0.9820
GTRIG	65	-1066.3596	2307.4066	0.0132	0.9952
GTRI	64	-1069.1543	2309.4174	0.0048	1.0000
K2PG	57	-1087.3160	2321.6921	0.0000	1.0000
SYMG	61	-1082.1963	2324.9845	0.0000	1.0000
K2PIG	58	-1087.3266	2325.0452	0.0000	1.0000
HKY	59	-1086.1126	2325.9828	0.0000	1.0000
K2PI	57	-1089.7078	2326.4756	0.0000	1.0000
SYMIG	62	-1082.1192	2328.3000	0.0000	1.0000
SYMI	61	-1084.9792	2330.5503	0.0000	1.0000
GTR	63	-1082.0633	2331.6937	0.0000	1.0000
K2P	56	-1102.6411	2349.0433	0.0000	1.0000
SYM	60	-1097.2741	2351.7056	0.0000	1.0000
F81G	59	-1102.0801	2357.9179	0.0000	1.0000
F81IG	60	-1101.9334	2361.0241	0.0000	1.0000
F81I	59	-1104.0486	2361.8547	0.0000	1.0000
JC69G	56	-1114.8794	2373.5200	0.0000	1.0000
JC69IG	57	-1115.2906	2377.6413	0.0000	1.0000
JC69I	56	-1117.5353	2378.8317	0.0000	1.0000
F81	58	-1116.9962	2384.3843	0.0000	1.0000
JC69	55	-1130.2533	2401.0016	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1066.5085	2261.0169	0.4356	0.4356
HKYG	60	-1071.1363	2262.2727	0.2325	0.6681
GTRIG	65	-1066.3596	2262.7191	0.1860	0.8540
HKYIG	61	-1071.0394	2264.0788	0.0942	0.9483
GTRI	64	-1069.1543	2266.3086	0.0309	0.9792
HKYI	60	-1073.5493	2267.0987	0.0208	1.0000
SYMG	61	-1082.1963	2286.3927	0.0000	1.0000
SYMIG	62	-1082.1192	2288.2385	0.0000	1.0000
K2PG	57	-1087.3160	2288.6321	0.0000	1.0000
GTR	63	-1082.0633	2290.1267	0.0000	1.0000
HKY	59	-1086.1126	2290.2252	0.0000	1.0000
K2PIG	58	-1087.3266	2290.6532	0.0000	1.0000
SYMI	61	-1084.9792	2291.9584	0.0000	1.0000
K2PI	57	-1089.7078	2293.4156	0.0000	1.0000
SYM	60	-1097.2741	2314.5482	0.0000	1.0000
K2P	56	-1102.6411	2317.2821	0.0000	1.0000
F81G	59	-1102.0801	2322.1603	0.0000	1.0000
F81IG	60	-1101.9334	2323.8668	0.0000	1.0000
F81I	59	-1104.0486	2326.0971	0.0000	1.0000
JC69G	56	-1114.8794	2341.7588	0.0000	1.0000
JC69IG	57	-1115.2906	2344.5813	0.0000	1.0000
JC69I	56	-1117.5353	2347.0705	0.0000	1.0000
F81	58	-1116.9962	2349.9924	0.0000	1.0000
JC69	55	-1130.2533	2370.5065	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1071.1363	2475.4503	0.8618	0.8618
HKYI	60	-1073.5493	2480.2762	0.0772	0.9390
HKYIG	61	-1071.0394	2480.8093	0.0591	0.9981
GTRG	64	-1066.5085	2488.4063	0.0013	0.9994
K2PG	57	-1087.3160	2491.1508	0.0003	0.9998
GTRIG	65	-1066.3596	2493.6615	0.0001	0.9998
GTRI	64	-1069.1543	2493.6980	0.0001	0.9999
K2PI	57	-1089.7078	2495.9343	0.0000	1.0000
K2PIG	58	-1087.3266	2496.7249	0.0000	1.0000
HKY	59	-1086.1126	2499.8498	0.0000	1.0000
SYMG	61	-1082.1963	2503.1232	0.0000	1.0000
SYMIG	62	-1082.1192	2508.5220	0.0000	1.0000
SYMI	61	-1084.9792	2508.6890	0.0000	1.0000
GTR	63	-1082.0633	2513.9632	0.0000	1.0000
K2P	56	-1102.6411	2516.2479	0.0000	1.0000
SYM	60	-1097.2741	2527.7258	0.0000	1.0000
F81G	59	-1102.0801	2531.7849	0.0000	1.0000
F81I	59	-1104.0486	2535.7218	0.0000	1.0000
F81IG	60	-1101.9334	2537.0444	0.0000	1.0000
JC69G	56	-1114.8794	2540.7245	0.0000	1.0000
JC69I	56	-1117.5353	2546.0363	0.0000	1.0000
JC69IG	57	-1115.2906	2547.1000	0.0000	1.0000
F81	58	-1116.9962	2556.0640	0.0000	1.0000
JC69	55	-1130.2533	2565.9193	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/1265_hg18.phylip = 1 - 258;
 Partition Dummy = 1:261-aggressive-phylip/1265_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 (HKYG)]

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


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

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

