
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:13 04/04/2011

Input data from file "683-tomodel-phylip/chr2_26639.phylip" (7 taxa, 1020 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3221.4447	6482.8895	0.5896	0.5896
GTRIG	21	-3221.4446	6484.8892	0.2169	0.8065
GTRI	20	-3223.1334	6486.2667	0.1089	0.9154
HKYG	16	-3227.8635	6487.7270	0.0525	0.9679
HKYIG	17	-3227.8737	6489.7474	0.0191	0.9870
HKYI	16	-3229.2607	6490.5214	0.0130	1.0000
GTR	19	-3241.6169	6521.2337	0.0000	1.0000
HKY	15	-3247.5494	6525.0989	0.0000	1.0000
K2PG	13	-3307.0584	6640.1169	0.0000	1.0000
K2PIG	14	-3307.1098	6642.2197	0.0000	1.0000
K2PI	13	-3308.2369	6642.4738	0.0000	1.0000
SYMG	17	-3306.1264	6646.2527	0.0000	1.0000
SYMIG	18	-3306.1429	6648.2858	0.0000	1.0000
SYMI	17	-3307.3742	6648.7484	0.0000	1.0000
F81G	15	-3311.3294	6652.6587	0.0000	1.0000
F81IG	16	-3311.3395	6654.6789	0.0000	1.0000
F81I	15	-3312.3838	6654.7676	0.0000	1.0000
K2P	12	-3329.1744	6682.3488	0.0000	1.0000
F81	14	-3329.8986	6687.7971	0.0000	1.0000
SYM	16	-3328.1843	6688.3687	0.0000	1.0000
JC69G	12	-3376.3963	6776.7927	0.0000	1.0000
JC69IG	13	-3376.4024	6778.8047	0.0000	1.0000
JC69I	12	-3377.4469	6778.8937	0.0000	1.0000
JC69	11	-3396.7083	6815.4166	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3227.8635	6566.5679	0.7819	0.7819
HKYI	16	-3229.2607	6569.3624	0.1933	0.9752
HKYIG	17	-3227.8737	6573.5159	0.0242	0.9994
GTRG	20	-3221.4447	6581.4406	0.0005	0.9999
GTRI	20	-3223.1334	6584.8179	0.0001	1.0000
GTRIG	21	-3221.4446	6588.3680	0.0000	1.0000
HKY	15	-3247.5494	6599.0123	0.0000	1.0000
GTR	19	-3241.6169	6614.8573	0.0000	1.0000
K2PG	13	-3307.0584	6704.1752	0.0000	1.0000
K2PI	13	-3308.2369	6706.5321	0.0000	1.0000
K2PIG	14	-3307.1098	6711.2055	0.0000	1.0000
F81G	15	-3311.3294	6726.5721	0.0000	1.0000
F81I	15	-3312.3838	6728.6810	0.0000	1.0000
SYMG	17	-3306.1264	6730.0212	0.0000	1.0000
SYMI	17	-3307.3742	6732.5169	0.0000	1.0000
F81IG	16	-3311.3395	6733.5198	0.0000	1.0000
SYMIG	18	-3306.1429	6736.9819	0.0000	1.0000
K2P	12	-3329.1744	6741.4795	0.0000	1.0000
F81	14	-3329.8986	6756.7830	0.0000	1.0000
SYM	16	-3328.1843	6767.2096	0.0000	1.0000
JC69G	12	-3376.3963	6835.9234	0.0000	1.0000
JC69I	12	-3377.4469	6838.0244	0.0000	1.0000
JC69IG	13	-3376.4024	6842.8630	0.0000	1.0000
JC69	11	-3396.7083	6869.6197	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 683-tomodel-phylip/chr2_26639.phylip = 1 - 1020;
 Partition Dummy = 1:683-tomodel-phylip/chr2_26639.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 683-tomodel-phylip/chr2_26639.phylip = 1 - 1020;
 Partition Dummy = 1:683-tomodel-phylip/chr2_26639.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 (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_26639.phylip = 1 - 1020;
 Partition Dummy = 1:683-tomodel-phylip/chr2_26639.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

