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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3838.5622	7717.1244	0.4379	0.4379
GTRI	20	-3838.6935	7717.3870	0.3840	0.8218
GTRIG	21	-3838.4623	7718.9247	0.1780	0.9998
HKYI	16	-3851.2801	7734.5603	0.0001	0.9999
HKYG	16	-3851.3000	7734.6000	0.0001	1.0000
HKYIG	17	-3851.1341	7736.2681	0.0000	1.0000
GTR	19	-3852.1262	7742.2525	0.0000	1.0000
HKY	15	-3865.5261	7761.0523	0.0000	1.0000
SYMG	17	-3871.1671	7776.3343	0.0000	1.0000
SYMI	17	-3871.2334	7776.4669	0.0000	1.0000
SYMIG	18	-3871.0466	7778.0933	0.0000	1.0000
K2PI	13	-3879.4867	7784.9734	0.0000	1.0000
K2PG	13	-3879.4957	7784.9914	0.0000	1.0000
K2PIG	14	-3879.3355	7786.6710	0.0000	1.0000
SYM	16	-3884.7146	7801.4292	0.0000	1.0000
K2P	12	-3893.8829	7811.7658	0.0000	1.0000
F81I	15	-3970.6241	7971.2482	0.0000	1.0000
F81G	15	-3971.0987	7972.1974	0.0000	1.0000
F81IG	16	-3970.6225	7973.2450	0.0000	1.0000
F81	14	-3982.2380	7992.4760	0.0000	1.0000
JC69I	12	-3994.9872	8013.9744	0.0000	1.0000
JC69G	12	-3995.4305	8014.8610	0.0000	1.0000
JC69IG	13	-3994.9794	8015.9588	0.0000	1.0000
JC69	11	-4007.1318	8036.2637	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3851.2801	7815.0257	0.3977	0.3977
HKYG	16	-3851.3000	7815.0654	0.3898	0.7875
GTRG	20	-3838.5622	7817.7062	0.1041	0.8916
GTRI	20	-3838.6935	7817.9688	0.0913	0.9829
HKYIG	17	-3851.1341	7821.7626	0.0137	0.9966
GTRIG	21	-3838.4623	7824.5355	0.0034	1.0000
HKY	15	-3865.5261	7836.4886	0.0000	1.0000
GTR	19	-3852.1262	7837.8051	0.0000	1.0000
K2PI	13	-3879.4867	7850.3516	0.0000	1.0000
K2PG	13	-3879.4957	7850.3695	0.0000	1.0000
K2PIG	14	-3879.3355	7857.0782	0.0000	1.0000
SYMG	17	-3871.1671	7861.8288	0.0000	1.0000
SYMI	17	-3871.2334	7861.9614	0.0000	1.0000
SYMIG	18	-3871.0466	7868.6168	0.0000	1.0000
K2P	12	-3893.8829	7872.1149	0.0000	1.0000
SYM	16	-3884.7146	7881.8946	0.0000	1.0000
F81I	15	-3970.6241	8046.6845	0.0000	1.0000
F81G	15	-3971.0987	8047.6337	0.0000	1.0000
F81IG	16	-3970.6225	8053.7104	0.0000	1.0000
F81	14	-3982.2380	8062.8832	0.0000	1.0000
JC69I	12	-3994.9872	8074.3234	0.0000	1.0000
JC69G	12	-3995.4305	8075.2101	0.0000	1.0000
JC69IG	13	-3994.9794	8081.3370	0.0000	1.0000
JC69	11	-4007.1318	8091.5836	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_21409.phylip = 1 - 1129;
 Partition Dummy = 1:683-tomodel-phylip/chr2_21409.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_21409.phylip = 1 - 1129;
 Partition Dummy = 1:683-tomodel-phylip/chr2_21409.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 (HKYI)]

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

