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

Input data from file "683-tomodel-phylip/chr3_15081.phylip" (7 taxa, 957 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	-2755.1077	5550.2154	0.4672	0.4672
GTRG	20	-2755.3753	5550.7506	0.3575	0.8248
GTRIG	21	-2755.0949	5552.1899	0.1741	0.9989
HKYI	16	-2765.9639	5563.9279	0.0005	0.9994
HKYG	16	-2766.1014	5564.2027	0.0004	0.9998
HKYIG	17	-2765.9268	5565.8536	0.0002	1.0000
GTR	19	-2766.2920	5570.5840	0.0000	1.0000
HKY	15	-2776.2510	5582.5021	0.0000	1.0000
F81I	15	-2816.5808	5663.1615	0.0000	1.0000
F81G	15	-2816.8731	5663.7462	0.0000	1.0000
F81IG	16	-2816.5838	5665.1677	0.0000	1.0000
F81	14	-2825.2519	5678.5037	0.0000	1.0000
SYMG	17	-2879.1810	5792.3619	0.0000	1.0000
SYMI	17	-2879.3133	5792.6266	0.0000	1.0000
SYMIG	18	-2879.1414	5794.2827	0.0000	1.0000
K2PI	13	-2888.4285	5802.8570	0.0000	1.0000
K2PG	13	-2888.4763	5802.9526	0.0000	1.0000
K2PIG	14	-2888.3647	5804.7294	0.0000	1.0000
SYM	16	-2890.7577	5813.5154	0.0000	1.0000
K2P	12	-2899.3797	5822.7594	0.0000	1.0000
JC69I	12	-2923.4918	5870.9835	0.0000	1.0000
JC69G	12	-2923.6870	5871.3741	0.0000	1.0000
JC69IG	13	-2923.4815	5872.9629	0.0000	1.0000
JC69	11	-2933.0656	5888.1311	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2765.9639	5641.7487	0.4980	0.4980
HKYG	16	-2766.1014	5642.0236	0.4341	0.9321
GTRI	20	-2755.1077	5647.4915	0.0282	0.9603
GTRG	20	-2755.3753	5648.0266	0.0216	0.9818
HKYIG	17	-2765.9268	5648.5382	0.0167	0.9985
GTRIG	21	-2755.0949	5654.3298	0.0009	0.9995
HKY	15	-2776.2510	5655.4592	0.0005	1.0000
GTR	19	-2766.2920	5662.9963	0.0000	1.0000
F81I	15	-2816.5808	5736.1186	0.0000	1.0000
F81G	15	-2816.8731	5736.7033	0.0000	1.0000
F81IG	16	-2816.5838	5742.9886	0.0000	1.0000
F81	14	-2825.2519	5746.5970	0.0000	1.0000
K2PI	13	-2888.4285	5866.0864	0.0000	1.0000
K2PG	13	-2888.4763	5866.1820	0.0000	1.0000
K2PIG	14	-2888.3647	5872.8226	0.0000	1.0000
SYMG	17	-2879.1810	5875.0466	0.0000	1.0000
SYMI	17	-2879.3133	5875.3112	0.0000	1.0000
K2P	12	-2899.3797	5881.1250	0.0000	1.0000
SYMIG	18	-2879.1414	5881.8312	0.0000	1.0000
SYM	16	-2890.7577	5891.3363	0.0000	1.0000
JC69I	12	-2923.4918	5929.3491	0.0000	1.0000
JC69G	12	-2923.6870	5929.7397	0.0000	1.0000
JC69IG	13	-2923.4815	5936.1924	0.0000	1.0000
JC69	11	-2933.0656	5941.6330	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 683-tomodel-phylip/chr3_15081.phylip = 1 - 957;
 Partition Dummy = 1:683-tomodel-phylip/chr3_15081.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 683-tomodel-phylip/chr3_15081.phylip = 1 - 957;
 Partition Dummy = 1:683-tomodel-phylip/chr3_15081.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 683-tomodel-phylip/chr3_15081.phylip = 1 - 957;
 Partition Dummy = 1:683-tomodel-phylip/chr3_15081.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

