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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2830.2259	5692.4518	0.3820	0.3820
HKYI	16	-2830.2658	5692.5317	0.3671	0.7491
HKYIG	17	-2830.2259	5694.4518	0.1405	0.8896
GTRG	20	-2828.2919	5696.5838	0.0484	0.9380
GTRI	20	-2828.3859	5696.7719	0.0441	0.9821
GTRIG	21	-2828.2844	5698.5689	0.0179	1.0000
K2PI	13	-2856.2008	5738.4015	0.0000	1.0000
K2PG	13	-2856.2369	5738.4739	0.0000	1.0000
K2PIG	14	-2856.2369	5740.4738	0.0000	1.0000
SYMI	17	-2854.5575	5743.1150	0.0000	1.0000
SYMG	17	-2854.5927	5743.1853	0.0000	1.0000
SYMIG	18	-2854.5926	5745.1853	0.0000	1.0000
HKY	15	-2865.8942	5761.7885	0.0000	1.0000
GTR	19	-2863.5451	5765.0902	0.0000	1.0000
K2P	12	-2893.8454	5811.6908	0.0000	1.0000
F81I	15	-2892.4597	5814.9194	0.0000	1.0000
F81G	15	-2892.5116	5815.0232	0.0000	1.0000
SYM	16	-2892.3366	5816.6733	0.0000	1.0000
F81IG	16	-2892.5117	5817.0233	0.0000	1.0000
JC69I	12	-2915.8321	5855.6643	0.0000	1.0000
JC69G	12	-2915.9143	5855.8285	0.0000	1.0000
JC69IG	13	-2915.9143	5857.8286	0.0000	1.0000
F81	14	-2925.6315	5879.2630	0.0000	1.0000
JC69	11	-2950.1300	5922.2600	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2830.2259	5773.8931	0.5026	0.5026
HKYI	16	-2830.2658	5773.9729	0.4829	0.9855
HKYIG	17	-2830.2259	5780.9831	0.0145	1.0000
GTRG	20	-2828.2919	5798.3854	0.0000	1.0000
GTRI	20	-2828.3859	5798.5734	0.0000	1.0000
K2PI	13	-2856.2008	5804.5725	0.0000	1.0000
K2PG	13	-2856.2369	5804.6449	0.0000	1.0000
GTRIG	21	-2828.2844	5805.4605	0.0000	1.0000
K2PIG	14	-2856.2369	5811.7349	0.0000	1.0000
SYMI	17	-2854.5575	5829.6463	0.0000	1.0000
SYMG	17	-2854.5927	5829.7166	0.0000	1.0000
SYMIG	18	-2854.5926	5836.8066	0.0000	1.0000
HKY	15	-2865.8942	5838.1396	0.0000	1.0000
GTR	19	-2863.5451	5861.8016	0.0000	1.0000
K2P	12	-2893.8454	5872.7717	0.0000	1.0000
F81I	15	-2892.4597	5891.2706	0.0000	1.0000
F81G	15	-2892.5116	5891.3744	0.0000	1.0000
SYM	16	-2892.3366	5898.1145	0.0000	1.0000
F81IG	16	-2892.5117	5898.4645	0.0000	1.0000
JC69I	12	-2915.8321	5916.7452	0.0000	1.0000
JC69G	12	-2915.9143	5916.9094	0.0000	1.0000
JC69IG	13	-2915.9143	5923.9996	0.0000	1.0000
F81	14	-2925.6315	5950.5240	0.0000	1.0000
JC69	11	-2950.1300	5978.2508	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_14888.phylip = 1 - 1200;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14888.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 AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_14888.phylip = 1 - 1200;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14888.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 591-tomodel-phylip/chr1_14888.phylip = 1 - 1200;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14888.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

