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

Input data from file "261-aggressive-phylip/1175_hg18.phylip" (29 taxa, 758 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	-3416.5627	6973.1312	0.9526	0.9526
GTRIG	65	-3418.3659	6979.1307	0.0474	1.0000
SYMG	61	-3435.2503	7003.3684	0.0000	1.0000
SYMIG	62	-3435.2591	7005.7586	0.0000	1.0000
HKYG	60	-3438.1492	7006.8006	0.0000	1.0000
HKYIG	61	-3438.1581	7009.1841	0.0000	1.0000
K2PG	57	-3445.9266	7015.2988	0.0000	1.0000
K2PIG	58	-3445.9355	7017.6622	0.0000	1.0000
GTRI	64	-3484.2298	7108.4654	0.0000	1.0000
SYMI	61	-3500.5021	7133.8721	0.0000	1.0000
HKYI	60	-3504.9877	7140.4775	0.0000	1.0000
F81G	59	-3508.3826	7144.9085	0.0000	1.0000
F81IG	60	-3508.3917	7147.2855	0.0000	1.0000
K2PI	57	-3512.8498	7149.1454	0.0000	1.0000
GTR	63	-3506.8249	7151.2694	0.0000	1.0000
JC69G	56	-3515.5146	7152.1362	0.0000	1.0000
JC69IG	57	-3515.5236	7154.4930	0.0000	1.0000
SYM	60	-3522.3556	7175.2133	0.0000	1.0000
HKY	59	-3527.9100	7183.9633	0.0000	1.0000
K2P	56	-3536.5849	7194.2767	0.0000	1.0000
F81I	59	-3573.1395	7274.4223	0.0000	1.0000
JC69I	56	-3579.9927	7281.0925	0.0000	1.0000
F81	58	-3595.9886	7317.7683	0.0000	1.0000
JC69	55	-3602.7933	7324.3615	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3416.5627	6961.1254	0.9428	0.9428
GTRIG	65	-3418.3659	6966.7318	0.0571	1.0000
SYMG	61	-3435.2503	6992.5006	0.0000	1.0000
SYMIG	62	-3435.2591	6994.5183	0.0000	1.0000
HKYG	60	-3438.1492	6996.2985	0.0000	1.0000
HKYIG	61	-3438.1581	6998.3163	0.0000	1.0000
K2PG	57	-3445.9266	7005.8531	0.0000	1.0000
K2PIG	58	-3445.9355	7007.8710	0.0000	1.0000
GTRI	64	-3484.2298	7096.4596	0.0000	1.0000
SYMI	61	-3500.5021	7123.0042	0.0000	1.0000
HKYI	60	-3504.9877	7129.9754	0.0000	1.0000
F81G	59	-3508.3826	7134.7653	0.0000	1.0000
F81IG	60	-3508.3917	7136.7834	0.0000	1.0000
GTR	63	-3506.8249	7139.6498	0.0000	1.0000
K2PI	57	-3512.8498	7139.6997	0.0000	1.0000
JC69G	56	-3515.5146	7143.0292	0.0000	1.0000
JC69IG	57	-3515.5236	7145.0473	0.0000	1.0000
SYM	60	-3522.3556	7164.7111	0.0000	1.0000
HKY	59	-3527.9100	7173.8201	0.0000	1.0000
K2P	56	-3536.5849	7185.1697	0.0000	1.0000
F81I	59	-3573.1395	7264.2790	0.0000	1.0000
JC69I	56	-3579.9927	7271.9855	0.0000	1.0000
F81	58	-3595.9886	7307.9772	0.0000	1.0000
JC69	55	-3602.7933	7315.5865	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-3416.5627	7257.4891	0.9915	0.9915
GTRIG	65	-3418.3659	7267.7263	0.0059	0.9974
K2PG	57	-3445.9266	7269.8021	0.0021	0.9995
HKYG	60	-3438.1492	7274.1395	0.0002	0.9998
SYMG	61	-3435.2503	7274.9723	0.0002	0.9999
K2PIG	58	-3445.9355	7276.4507	0.0001	1.0000
HKYIG	61	-3438.1581	7280.7879	0.0000	1.0000
SYMIG	62	-3435.2591	7281.6206	0.0000	1.0000
GTRI	64	-3484.2298	7392.8234	0.0000	1.0000
JC69G	56	-3515.5146	7402.3475	0.0000	1.0000
K2PI	57	-3512.8498	7403.6486	0.0000	1.0000
SYMI	61	-3500.5021	7405.4759	0.0000	1.0000
HKYI	60	-3504.9877	7407.8164	0.0000	1.0000
F81G	59	-3508.3826	7407.9756	0.0000	1.0000
JC69IG	57	-3515.5236	7408.9962	0.0000	1.0000
F81IG	60	-3508.3917	7414.6244	0.0000	1.0000
GTR	63	-3506.8249	7431.3829	0.0000	1.0000
SYM	60	-3522.3556	7442.5521	0.0000	1.0000
K2P	56	-3536.5849	7444.4880	0.0000	1.0000
HKY	59	-3527.9100	7447.0304	0.0000	1.0000
JC69I	56	-3579.9927	7531.3038	0.0000	1.0000
F81I	59	-3573.1395	7537.4894	0.0000	1.0000
JC69	55	-3602.7933	7570.2741	0.0000	1.0000
F81	58	-3595.9886	7576.5568	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/1175_hg18.phylip = 1 - 758;
 Partition Dummy = 1:261-aggressive-phylip/1175_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/1175_hg18.phylip = 1 - 758;
 Partition Dummy = 1:261-aggressive-phylip/1175_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/1175_hg18.phylip = 1 - 758;
 Partition Dummy = 1:261-aggressive-phylip/1175_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

