
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:37 04/04/2011

Input data from file "591-tomodel-phylip/chr2_21324.phylip" (7 taxa, 1089 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-3314.7343	6669.4687	0.5489	0.5489
GTRG	20	-3315.7094	6671.4188	0.2070	0.7559
GTRIG	21	-3314.7362	6671.4724	0.2016	0.9575
HKYI	16	-3321.8819	6675.7639	0.0236	0.9811
HKYG	16	-3322.7143	6677.4286	0.0103	0.9913
HKYIG	17	-3321.8837	6677.7675	0.0087	1.0000
SYMI	17	-3349.1713	6732.3425	0.0000	1.0000
SYMG	17	-3350.0824	6734.1649	0.0000	1.0000
K2PI	13	-3354.1566	6734.3131	0.0000	1.0000
K2PG	13	-3355.0489	6736.0979	0.0000	1.0000
K2PIG	14	-3354.1605	6736.3209	0.0000	1.0000
SYMIG	18	-3350.1654	6736.3307	0.0000	1.0000
GTR	19	-3356.1376	6750.2752	0.0000	1.0000
HKY	15	-3361.2501	6752.5002	0.0000	1.0000
F81I	15	-3383.0327	6796.0655	0.0000	1.0000
F81IG	16	-3383.0664	6798.1327	0.0000	1.0000
F81G	15	-3384.7889	6799.5778	0.0000	1.0000
SYM	16	-3385.7870	6803.5739	0.0000	1.0000
K2P	12	-3390.5957	6805.1914	0.0000	1.0000
JC69I	12	-3408.7551	6841.5102	0.0000	1.0000
JC69IG	13	-3408.7734	6843.5469	0.0000	1.0000
JC69G	12	-3410.1130	6844.2261	0.0000	1.0000
F81	14	-3416.5474	6861.0948	0.0000	1.0000
JC69	11	-3441.9136	6905.8272	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3321.8819	6755.6521	0.6818	0.6818
HKYG	16	-3322.7143	6757.3168	0.2966	0.9783
HKYIG	17	-3321.8837	6762.6488	0.0206	0.9990
GTRI	20	-3314.7343	6769.3290	0.0007	0.9997
GTRG	20	-3315.7094	6771.2791	0.0003	1.0000
GTRIG	21	-3314.7362	6776.3257	0.0000	1.0000
K2PI	13	-3354.1566	6799.2223	0.0000	1.0000
K2PG	13	-3355.0489	6801.0071	0.0000	1.0000
K2PIG	14	-3354.1605	6806.2232	0.0000	1.0000
SYMI	17	-3349.1713	6817.2238	0.0000	1.0000
SYMG	17	-3350.0824	6819.0461	0.0000	1.0000
SYMIG	18	-3350.1654	6826.2050	0.0000	1.0000
HKY	15	-3361.2501	6827.3954	0.0000	1.0000
GTR	19	-3356.1376	6845.1425	0.0000	1.0000
K2P	12	-3390.5957	6865.1076	0.0000	1.0000
F81I	15	-3383.0327	6870.9607	0.0000	1.0000
F81G	15	-3384.7889	6874.4730	0.0000	1.0000
F81IG	16	-3383.0664	6878.0210	0.0000	1.0000
SYM	16	-3385.7870	6883.4622	0.0000	1.0000
JC69I	12	-3408.7551	6901.4264	0.0000	1.0000
JC69G	12	-3410.1130	6904.1423	0.0000	1.0000
JC69IG	13	-3408.7734	6908.4561	0.0000	1.0000
F81	14	-3416.5474	6930.9971	0.0000	1.0000
JC69	11	-3441.9136	6960.7504	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 591-tomodel-phylip/chr2_21324.phylip = 1 - 1089;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21324.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 591-tomodel-phylip/chr2_21324.phylip = 1 - 1089;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21324.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_21324.phylip = 1 - 1089;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21324.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

