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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	64	-1444.4702	3031.2114	0.4080	0.4080
GTRIG	65	-1443.5471	3031.8365	0.2985	0.7065
GTRG	64	-1444.7995	3031.8701	0.2935	1.0000
HKYG	60	-1464.8900	3062.2502	0.0000	1.0000
HKYI	60	-1470.1141	3072.6984	0.0000	1.0000
HKYIG	61	-1469.6493	3074.2064	0.0000	1.0000
GTR	63	-1471.7785	3083.3652	0.0000	1.0000
HKY	59	-1491.4404	3112.9216	0.0000	1.0000
F81G	59	-1493.4698	3116.9804	0.0000	1.0000
F81I	59	-1496.1193	3122.2794	0.0000	1.0000
F81IG	60	-1495.3411	3123.1524	0.0000	1.0000
SYMI	61	-1505.6322	3146.1723	0.0000	1.0000
SYMG	61	-1507.1145	3149.1368	0.0000	1.0000
K2PI	57	-1512.8068	3150.8204	0.0000	1.0000
SYMIG	62	-1506.9885	3151.3308	0.0000	1.0000
K2PIG	58	-1514.8988	3157.4172	0.0000	1.0000
K2PG	57	-1516.4091	3158.0249	0.0000	1.0000
F81	58	-1519.5507	3166.7212	0.0000	1.0000
SYM	60	-1532.5526	3197.5753	0.0000	1.0000
JC69I	56	-1537.7172	3198.2365	0.0000	1.0000
JC69IG	57	-1537.3177	3199.8422	0.0000	1.0000
JC69G	56	-1538.8259	3200.4539	0.0000	1.0000
K2P	56	-1543.2098	3209.2217	0.0000	1.0000
JC69	55	-1565.2852	3250.9758	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-1444.4702	3016.9403	0.3780	0.3780
GTRIG	65	-1443.5471	3017.0943	0.3500	0.7281
GTRG	64	-1444.7995	3017.5991	0.2719	1.0000
HKYG	60	-1464.8900	3049.7800	0.0000	1.0000
HKYI	60	-1470.1141	3060.2282	0.0000	1.0000
HKYIG	61	-1469.6493	3061.2986	0.0000	1.0000
GTR	63	-1471.7785	3069.5570	0.0000	1.0000
HKY	59	-1491.4404	3100.8808	0.0000	1.0000
F81G	59	-1493.4698	3104.9396	0.0000	1.0000
F81I	59	-1496.1193	3110.2385	0.0000	1.0000
F81IG	60	-1495.3411	3110.6822	0.0000	1.0000
SYMI	61	-1505.6322	3133.2645	0.0000	1.0000
SYMG	61	-1507.1145	3136.2290	0.0000	1.0000
SYMIG	62	-1506.9885	3137.9770	0.0000	1.0000
K2PI	57	-1512.8068	3139.6136	0.0000	1.0000
K2PIG	58	-1514.8988	3145.7975	0.0000	1.0000
K2PG	57	-1516.4091	3146.8181	0.0000	1.0000
F81	58	-1519.5507	3155.1015	0.0000	1.0000
SYM	60	-1532.5526	3185.1052	0.0000	1.0000
JC69I	56	-1537.7172	3187.4345	0.0000	1.0000
JC69IG	57	-1537.3177	3188.6354	0.0000	1.0000
JC69G	56	-1538.8259	3189.6518	0.0000	1.0000
K2P	56	-1543.2098	3198.4197	0.0000	1.0000
JC69	55	-1565.2852	3240.5704	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	64	-1444.4702	3303.2693	0.5498	0.5498
GTRG	64	-1444.7995	3303.9281	0.3955	0.9453
GTRIG	65	-1443.5471	3307.8972	0.0544	0.9997
HKYG	60	-1464.8900	3318.2134	0.0003	1.0000
HKYI	60	-1470.1141	3328.6616	0.0000	1.0000
HKYIG	61	-1469.6493	3334.2059	0.0000	1.0000
GTR	63	-1471.7785	3351.4121	0.0000	1.0000
HKY	59	-1491.4404	3364.8404	0.0000	1.0000
F81G	59	-1493.4698	3368.8992	0.0000	1.0000
F81I	59	-1496.1193	3374.1981	0.0000	1.0000
F81IG	60	-1495.3411	3379.1157	0.0000	1.0000
K2PI	57	-1512.8068	3394.6254	0.0000	1.0000
K2PG	57	-1516.4091	3401.8299	0.0000	1.0000
K2PIG	58	-1514.8988	3405.2832	0.0000	1.0000
SYMI	61	-1505.6322	3406.1718	0.0000	1.0000
SYMG	61	-1507.1145	3409.1363	0.0000	1.0000
F81	58	-1519.5507	3414.5871	0.0000	1.0000
SYMIG	62	-1506.9885	3415.3582	0.0000	1.0000
JC69I	56	-1537.7172	3437.9723	0.0000	1.0000
JC69G	56	-1538.8259	3440.1897	0.0000	1.0000
JC69IG	57	-1537.3177	3443.6472	0.0000	1.0000
K2P	56	-1543.2098	3448.9575	0.0000	1.0000
SYM	60	-1532.5526	3453.5386	0.0000	1.0000
JC69	55	-1565.2852	3486.6344	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6197_hg18.phylip = 1 - 648;
 Partition Dummy = 1:261-aggressive-phylip/6197_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6197_hg18.phylip = 1 - 648;
 Partition Dummy = 1:261-aggressive-phylip/6197_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6197_hg18.phylip = 1 - 648;
 Partition Dummy = 1:261-aggressive-phylip/6197_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

