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

Input data from file "683-tomodel-phylip/chr13_484.phylip" (7 taxa, 1039 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	-3227.9189	6495.8378	0.4172	0.4172
GTRG	20	-3228.4202	6496.8404	0.2527	0.6699
GTRIG	21	-3227.9225	6497.8449	0.1529	0.8228
HKYI	16	-3233.5786	6499.1571	0.0793	0.9022
HKYG	16	-3233.7430	6499.4860	0.0673	0.9695
HKYIG	17	-3233.5342	6501.0683	0.0305	1.0000
GTR	19	-3247.1375	6532.2750	0.0000	1.0000
HKY	15	-3252.8102	6535.6204	0.0000	1.0000
F81I	15	-3300.1975	6630.3949	0.0000	1.0000
F81G	15	-3300.3924	6630.7847	0.0000	1.0000
F81IG	16	-3300.1820	6632.3639	0.0000	1.0000
SYMG	17	-3310.7180	6655.4361	0.0000	1.0000
SYMI	17	-3311.1827	6656.3654	0.0000	1.0000
K2PG	13	-3315.6102	6657.2204	0.0000	1.0000
SYMIG	18	-3310.7419	6657.4839	0.0000	1.0000
K2PI	13	-3316.0153	6658.0306	0.0000	1.0000
K2PIG	14	-3315.6323	6659.2646	0.0000	1.0000
F81	14	-3318.5973	6665.1947	0.0000	1.0000
K2P	12	-3336.9052	6697.8104	0.0000	1.0000
SYM	16	-3333.1476	6698.2952	0.0000	1.0000
JC69G	12	-3370.4358	6764.8717	0.0000	1.0000
JC69I	12	-3370.6761	6765.3522	0.0000	1.0000
JC69IG	13	-3370.4514	6766.9028	0.0000	1.0000
JC69	11	-3390.2011	6802.4022	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3233.5786	6578.2934	0.5316	0.5316
HKYG	16	-3233.7430	6578.6222	0.4510	0.9825
HKYIG	17	-3233.5342	6585.1506	0.0172	0.9998
GTRI	20	-3227.9189	6594.7581	0.0001	0.9999
GTRG	20	-3228.4202	6595.7607	0.0001	1.0000
GTRIG	21	-3227.9225	6601.7112	0.0000	1.0000
HKY	15	-3252.8102	6609.8106	0.0000	1.0000
GTR	19	-3247.1375	6626.2493	0.0000	1.0000
F81I	15	-3300.1975	6704.5851	0.0000	1.0000
F81G	15	-3300.3924	6704.9749	0.0000	1.0000
F81IG	16	-3300.1820	6711.5001	0.0000	1.0000
K2PG	13	-3315.6102	6721.5186	0.0000	1.0000
K2PI	13	-3316.0153	6722.3288	0.0000	1.0000
K2PIG	14	-3315.6323	6728.5088	0.0000	1.0000
F81	14	-3318.5973	6734.4389	0.0000	1.0000
SYMG	17	-3310.7180	6739.5183	0.0000	1.0000
SYMI	17	-3311.1827	6740.4476	0.0000	1.0000
SYMIG	18	-3310.7419	6746.5121	0.0000	1.0000
K2P	12	-3336.9052	6757.1626	0.0000	1.0000
SYM	16	-3333.1476	6777.4314	0.0000	1.0000
JC69G	12	-3370.4358	6824.2238	0.0000	1.0000
JC69I	12	-3370.6761	6824.7043	0.0000	1.0000
JC69IG	13	-3370.4514	6831.2010	0.0000	1.0000
JC69	11	-3390.2011	6856.8084	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/chr13_484.phylip = 1 - 1039;
 Partition Dummy = 1:683-tomodel-phylip/chr13_484.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/chr13_484.phylip = 1 - 1039;
 Partition Dummy = 1:683-tomodel-phylip/chr13_484.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/chr13_484.phylip = 1 - 1039;
 Partition Dummy = 1:683-tomodel-phylip/chr13_484.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

