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

Input data from file "683-tomodel-phylip/chr1_22430.phylip" (7 taxa, 1117 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	-3114.9484	6261.8968	0.4243	0.4243
HKYG	16	-3115.1766	6262.3533	0.3377	0.7620
HKYIG	17	-3115.1886	6264.3772	0.1228	0.8848
GTRI	20	-3113.0390	6266.0780	0.0524	0.9372
GTRG	20	-3113.1662	6266.3323	0.0462	0.9834
GTRIG	21	-3113.1890	6268.3781	0.0166	1.0000
HKY	15	-3136.7598	6303.5197	0.0000	1.0000
GTR	19	-3135.3848	6308.7695	0.0000	1.0000
K2PIG	14	-3156.1756	6340.3512	0.0000	1.0000
K2PG	13	-3157.3503	6340.7006	0.0000	1.0000
K2PI	13	-3157.5804	6341.1609	0.0000	1.0000
SYMI	17	-3154.3670	6342.7340	0.0000	1.0000
SYMIG	18	-3153.8916	6343.7833	0.0000	1.0000
SYMG	17	-3155.0786	6344.1571	0.0000	1.0000
F81I	15	-3166.1475	6362.2951	0.0000	1.0000
F81G	15	-3166.2682	6362.5363	0.0000	1.0000
F81IG	16	-3166.2789	6364.5578	0.0000	1.0000
K2P	12	-3177.3845	6378.7691	0.0000	1.0000
SYM	16	-3175.2249	6382.4498	0.0000	1.0000
F81	14	-3187.5246	6403.0493	0.0000	1.0000
JC69IG	13	-3203.5342	6433.0685	0.0000	1.0000
JC69I	12	-3205.1722	6434.3444	0.0000	1.0000
JC69G	12	-3205.2520	6434.5040	0.0000	1.0000
JC69	11	-3225.5161	6473.0321	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3114.9484	6342.1912	0.5496	0.5496
HKYG	16	-3115.1766	6342.6477	0.4375	0.9871
HKYIG	17	-3115.1886	6349.6900	0.0129	1.0000
GTRI	20	-3113.0390	6366.4461	0.0000	1.0000
GTRG	20	-3113.1662	6366.7004	0.0000	1.0000
GTRIG	21	-3113.1890	6373.7645	0.0000	1.0000
HKY	15	-3136.7598	6378.7957	0.0000	1.0000
GTR	19	-3135.3848	6404.1192	0.0000	1.0000
K2PG	13	-3157.3503	6405.9398	0.0000	1.0000
K2PI	13	-3157.5804	6406.4001	0.0000	1.0000
K2PIG	14	-3156.1756	6410.6089	0.0000	1.0000
SYMI	17	-3154.3670	6428.0468	0.0000	1.0000
SYMG	17	-3155.0786	6429.4700	0.0000	1.0000
SYMIG	18	-3153.8916	6434.1145	0.0000	1.0000
F81I	15	-3166.1475	6437.5711	0.0000	1.0000
F81G	15	-3166.2682	6437.8123	0.0000	1.0000
K2P	12	-3177.3845	6438.9899	0.0000	1.0000
F81IG	16	-3166.2789	6444.8523	0.0000	1.0000
SYM	16	-3175.2249	6462.7443	0.0000	1.0000
F81	14	-3187.5246	6473.3069	0.0000	1.0000
JC69I	12	-3205.1722	6494.5652	0.0000	1.0000
JC69G	12	-3205.2520	6494.7249	0.0000	1.0000
JC69IG	13	-3203.5342	6498.3077	0.0000	1.0000
JC69	11	-3225.5161	6528.2345	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 683-tomodel-phylip/chr1_22430.phylip = 1 - 1117;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22430.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 683-tomodel-phylip/chr1_22430.phylip = 1 - 1117;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22430.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 683-tomodel-phylip/chr1_22430.phylip = 1 - 1117;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22430.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

