
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:55 04/04/2011

Input data from file "591-tomodel-phylip/chr8_5056.phylip" (7 taxa, 901 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-3287.3324	6606.6648	0.4943	0.4943
HKYG	16	-3288.1963	6608.3925	0.2084	0.7027
HKYIG	17	-3287.2284	6608.4569	0.2018	0.9045
GTRI	20	-3285.5275	6611.0550	0.0550	0.9595
GTRIG	21	-3285.4580	6612.9161	0.0217	0.9812
GTRG	20	-3286.6022	6613.2043	0.0188	1.0000
HKY	15	-3328.4664	6686.9328	0.0000	1.0000
SYMI	17	-3327.8175	6689.6349	0.0000	1.0000
GTR	19	-3326.0602	6690.1203	0.0000	1.0000
K2PI	13	-3332.4160	6690.8321	0.0000	1.0000
SYMG	17	-3328.5139	6691.0279	0.0000	1.0000
K2PG	13	-3332.7332	6691.4664	0.0000	1.0000
SYMIG	18	-3327.7428	6691.4856	0.0000	1.0000
K2PIG	14	-3332.2389	6692.4779	0.0000	1.0000
F81I	15	-3369.3197	6768.6393	0.0000	1.0000
F81IG	16	-3369.3220	6770.6440	0.0000	1.0000
F81G	15	-3370.8921	6771.7841	0.0000	1.0000
K2P	12	-3376.4486	6776.8973	0.0000	1.0000
SYM	16	-3372.7747	6777.5494	0.0000	1.0000
JC69I	12	-3407.4664	6838.9329	0.0000	1.0000
JC69IG	13	-3407.4673	6840.9347	0.0000	1.0000
JC69G	12	-3408.8528	6841.7056	0.0000	1.0000
F81	14	-3408.4407	6844.8814	0.0000	1.0000
JC69	11	-3448.2692	6918.5384	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3287.3324	6683.5209	0.6856	0.6856
HKYG	16	-3288.1963	6685.2486	0.2890	0.9746
HKYIG	17	-3287.2284	6690.1165	0.0253	1.0000
GTRI	20	-3285.5275	6707.1251	0.0000	1.0000
GTRG	20	-3286.6022	6709.2744	0.0000	1.0000
GTRIG	21	-3285.4580	6713.7897	0.0000	1.0000
K2PI	13	-3332.4160	6753.2776	0.0000	1.0000
K2PG	13	-3332.7332	6753.9120	0.0000	1.0000
HKY	15	-3328.4664	6758.9853	0.0000	1.0000
K2PIG	14	-3332.2389	6759.7270	0.0000	1.0000
SYMI	17	-3327.8175	6771.2945	0.0000	1.0000
SYMG	17	-3328.5139	6772.6874	0.0000	1.0000
SYMIG	18	-3327.7428	6777.9487	0.0000	1.0000
GTR	19	-3326.0602	6781.3869	0.0000	1.0000
K2P	12	-3376.4486	6834.5393	0.0000	1.0000
F81I	15	-3369.3197	6840.6919	0.0000	1.0000
F81G	15	-3370.8921	6843.8367	0.0000	1.0000
F81IG	16	-3369.3220	6847.5000	0.0000	1.0000
SYM	16	-3372.7747	6854.4055	0.0000	1.0000
JC69I	12	-3407.4664	6896.5749	0.0000	1.0000
JC69G	12	-3408.8528	6899.3476	0.0000	1.0000
JC69IG	13	-3407.4673	6903.3802	0.0000	1.0000
F81	14	-3408.4407	6912.1305	0.0000	1.0000
JC69	11	-3448.2692	6971.3769	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_5056.phylip = 1 - 901;
 Partition Dummy = 1:591-tomodel-phylip/chr8_5056.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;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_5056.phylip = 1 - 901;
 Partition Dummy = 1:591-tomodel-phylip/chr8_5056.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;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_5056.phylip = 1 - 901;
 Partition Dummy = 1:591-tomodel-phylip/chr8_5056.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

