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

Input data from file "683-tomodel-phylip/chr1_14672.phylip" (7 taxa, 1113 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	-3860.8805	7761.7610	0.3919	0.3919
HKYI	16	-3865.6581	7763.3163	0.1801	0.5720
GTRG	20	-3861.8362	7763.6725	0.1507	0.7228
GTRIG	21	-3860.9003	7763.8006	0.1414	0.8641
HKYG	16	-3866.5915	7765.1830	0.0708	0.9349
HKYIG	17	-3865.6781	7765.3563	0.0649	0.9999
GTR	19	-3870.3520	7778.7039	0.0001	1.0000
HKY	15	-3874.9852	7779.9704	0.0000	1.0000
K2PI	13	-3894.3281	7814.6562	0.0000	1.0000
SYMI	17	-3890.6232	7815.2464	0.0000	1.0000
K2PG	13	-3895.1217	7816.2433	0.0000	1.0000
SYMG	17	-3891.3100	7816.6200	0.0000	1.0000
K2PIG	14	-3894.3459	7816.6918	0.0000	1.0000
SYMIG	18	-3890.6372	7817.2745	0.0000	1.0000
SYM	16	-3898.3727	7828.7453	0.0000	1.0000
K2P	12	-3902.6152	7829.2303	0.0000	1.0000
F81I	15	-3968.1311	7966.2621	0.0000	1.0000
F81IG	16	-3968.1593	7968.3186	0.0000	1.0000
F81G	15	-3969.1717	7968.3433	0.0000	1.0000
F81	14	-3974.9971	7977.9942	0.0000	1.0000
JC69I	12	-3989.7217	8003.4434	0.0000	1.0000
JC69G	12	-3990.6956	8005.3913	0.0000	1.0000
JC69IG	13	-3989.7483	8005.4965	0.0000	1.0000
JC69	11	-3996.3949	8014.7899	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3865.6581	7843.5533	0.7014	0.7014
HKYG	16	-3866.5915	7845.4200	0.2758	0.9772
HKYIG	17	-3865.6781	7850.6081	0.0206	0.9978
HKY	15	-3874.9852	7855.1926	0.0021	0.9999
GTRI	20	-3860.8805	7862.0573	0.0001	1.0000
GTRG	20	-3861.8362	7863.9687	0.0000	1.0000
GTRIG	21	-3860.9003	7869.1117	0.0000	1.0000
GTR	19	-3870.3520	7873.9854	0.0000	1.0000
K2PI	13	-3894.3281	7879.8487	0.0000	1.0000
K2PG	13	-3895.1217	7881.4359	0.0000	1.0000
K2PIG	14	-3894.3459	7886.8992	0.0000	1.0000
K2P	12	-3902.6152	7889.4081	0.0000	1.0000
SYMI	17	-3890.6232	7900.4982	0.0000	1.0000
SYMG	17	-3891.3100	7901.8719	0.0000	1.0000
SYMIG	18	-3890.6372	7907.5411	0.0000	1.0000
SYM	16	-3898.3727	7908.9823	0.0000	1.0000
F81I	15	-3968.1311	8041.4843	0.0000	1.0000
F81G	15	-3969.1717	8043.5655	0.0000	1.0000
F81	14	-3974.9971	8048.2016	0.0000	1.0000
F81IG	16	-3968.1593	8048.5556	0.0000	1.0000
JC69I	12	-3989.7217	8063.6212	0.0000	1.0000
JC69G	12	-3990.6956	8065.5690	0.0000	1.0000
JC69	11	-3996.3949	8069.9528	0.0000	1.0000
JC69IG	13	-3989.7483	8070.6891	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/chr1_14672.phylip = 1 - 1113;
 Partition Dummy = 1:683-tomodel-phylip/chr1_14672.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/chr1_14672.phylip = 1 - 1113;
 Partition Dummy = 1:683-tomodel-phylip/chr1_14672.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/chr1_14672.phylip = 1 - 1113;
 Partition Dummy = 1:683-tomodel-phylip/chr1_14672.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

