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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-9491.6274	19116.9456	0.7524	0.7524
GTRIG	65	-9491.6481	19119.1688	0.2476	1.0000
HKYG	60	-9514.8898	19154.7729	0.0000	1.0000
HKYIG	61	-9514.9099	19156.9830	0.0000	1.0000
SYMG	61	-9595.1265	19317.4161	0.0000	1.0000
SYMIG	62	-9595.1455	19319.6271	0.0000	1.0000
K2PG	57	-9613.6826	19345.8662	0.0000	1.0000
K2PIG	58	-9613.7028	19348.0677	0.0000	1.0000
GTRI	64	-9634.2523	19402.1955	0.0000	1.0000
HKYI	60	-9660.5291	19446.0513	0.0000	1.0000
SYMI	61	-9744.7588	19616.6808	0.0000	1.0000
K2PI	57	-9766.2795	19651.0600	0.0000	1.0000
GTR	63	-9813.2107	19757.9334	0.0000	1.0000
F81G	59	-9836.2298	19795.2857	0.0000	1.0000
F81IG	60	-9836.2501	19797.4934	0.0000	1.0000
HKY	59	-9846.4352	19815.6966	0.0000	1.0000
JC69G	56	-9899.2406	19914.8240	0.0000	1.0000
JC69IG	57	-9899.2609	19917.0228	0.0000	1.0000
SYM	60	-9930.3770	19985.7473	0.0000	1.0000
K2P	56	-9958.2110	20032.7649	0.0000	1.0000
F81I	59	-9982.2220	20087.2702	0.0000	1.0000
JC69I	56	-10050.6153	20217.5735	0.0000	1.0000
F81	58	-10162.2315	20445.1251	0.0000	1.0000
JC69	55	-10235.2873	20584.7622	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-9491.6274	19111.2548	0.7351	0.7351
GTRIG	65	-9491.6481	19113.2962	0.2649	1.0000
HKYG	60	-9514.8898	19149.7797	0.0000	1.0000
HKYIG	61	-9514.9099	19151.8199	0.0000	1.0000
SYMG	61	-9595.1265	19312.2529	0.0000	1.0000
SYMIG	62	-9595.1455	19314.2910	0.0000	1.0000
K2PG	57	-9613.6826	19341.3652	0.0000	1.0000
K2PIG	58	-9613.7028	19343.4056	0.0000	1.0000
GTRI	64	-9634.2523	19396.5046	0.0000	1.0000
HKYI	60	-9660.5291	19441.0582	0.0000	1.0000
SYMI	61	-9744.7588	19611.5177	0.0000	1.0000
K2PI	57	-9766.2795	19646.5590	0.0000	1.0000
GTR	63	-9813.2107	19752.4215	0.0000	1.0000
F81G	59	-9836.2298	19790.4595	0.0000	1.0000
F81IG	60	-9836.2501	19792.5002	0.0000	1.0000
HKY	59	-9846.4352	19810.8705	0.0000	1.0000
JC69G	56	-9899.2406	19910.4811	0.0000	1.0000
JC69IG	57	-9899.2609	19912.5218	0.0000	1.0000
SYM	60	-9930.3770	19980.7541	0.0000	1.0000
K2P	56	-9958.2110	20028.4221	0.0000	1.0000
F81I	59	-9982.2220	20082.4440	0.0000	1.0000
JC69I	56	-10050.6153	20213.2306	0.0000	1.0000
F81	58	-10162.2315	20440.4630	0.0000	1.0000
JC69	55	-10235.2873	20580.5746	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-9491.6274	19452.4427	0.9754	0.9754
GTRIG	65	-9491.6481	19459.8151	0.0244	0.9998
HKYG	60	-9514.8898	19469.6433	0.0002	1.0000
HKYIG	61	-9514.9099	19477.0146	0.0000	1.0000
SYMG	61	-9595.1265	19637.4476	0.0000	1.0000
SYMIG	62	-9595.1455	19644.8168	0.0000	1.0000
K2PG	57	-9613.6826	19645.2356	0.0000	1.0000
K2PIG	58	-9613.7028	19652.6071	0.0000	1.0000
GTRI	64	-9634.2523	19737.6925	0.0000	1.0000
HKYI	60	-9660.5291	19760.9218	0.0000	1.0000
SYMI	61	-9744.7588	19936.7124	0.0000	1.0000
K2PI	57	-9766.2795	19950.4294	0.0000	1.0000
GTR	63	-9813.2107	20088.2783	0.0000	1.0000
F81G	59	-9836.2298	20104.9921	0.0000	1.0000
F81IG	60	-9836.2501	20112.3639	0.0000	1.0000
HKY	59	-9846.4352	20125.4030	0.0000	1.0000
JC69G	56	-9899.2406	20209.0205	0.0000	1.0000
JC69IG	57	-9899.2609	20216.3922	0.0000	1.0000
SYM	60	-9930.3770	20300.6177	0.0000	1.0000
K2P	56	-9958.2110	20326.9614	0.0000	1.0000
F81I	59	-9982.2220	20396.9765	0.0000	1.0000
JC69I	56	-10050.6153	20511.7700	0.0000	1.0000
F81	58	-10162.2315	20749.6645	0.0000	1.0000
JC69	55	-10235.2873	20873.7829	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/10292_hg18.phylip = 1 - 1527;
 Partition Dummy = 1:261-aggressive-phylip/10292_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/10292_hg18.phylip = 1 - 1527;
 Partition Dummy = 1:261-aggressive-phylip/10292_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 (GTRG)]

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


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

