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

Input data from file "591-tomodel-phylip/chr8_4233.phylip" (7 taxa, 895 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	-3384.2157	6800.4314	0.3603	0.3603
HKYG	16	-3384.3391	6800.6781	0.3185	0.6788
HKYIG	17	-3383.7564	6801.5129	0.2098	0.8886
GTRI	20	-3382.3318	6804.6636	0.0434	0.9320
GTRG	20	-3382.3764	6804.7527	0.0415	0.9735
GTRIG	21	-3381.8269	6805.6539	0.0265	1.0000
HKY	15	-3412.2968	6854.5937	0.0000	1.0000
GTR	19	-3410.4858	6858.9716	0.0000	1.0000
SYMG	17	-3419.8154	6873.6308	0.0000	1.0000
SYMI	17	-3420.1759	6874.3519	0.0000	1.0000
SYMIG	18	-3419.4867	6874.9734	0.0000	1.0000
K2PI	13	-3426.1085	6878.2171	0.0000	1.0000
K2PG	13	-3426.3008	6878.6016	0.0000	1.0000
K2PIG	14	-3425.7349	6879.4697	0.0000	1.0000
SYM	16	-3453.4528	6938.9056	0.0000	1.0000
K2P	12	-3457.9694	6939.9388	0.0000	1.0000
F81I	15	-3518.9191	7067.8382	0.0000	1.0000
F81IG	16	-3518.9213	7069.8426	0.0000	1.0000
F81G	15	-3520.5854	7071.1709	0.0000	1.0000
F81	14	-3544.2518	7116.5037	0.0000	1.0000
JC69I	12	-3553.7086	7131.4171	0.0000	1.0000
JC69IG	13	-3553.7078	7133.4155	0.0000	1.0000
JC69G	12	-3555.1537	7134.3074	0.0000	1.0000
JC69	11	-3580.8515	7183.7031	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3384.2157	6877.1805	0.5163	0.5163
HKYG	16	-3384.3391	6877.4273	0.4564	0.9727
HKYIG	17	-3383.7564	6883.0589	0.0273	1.0000
GTRI	20	-3382.3318	6900.6001	0.0000	1.0000
GTRG	20	-3382.3764	6900.6892	0.0000	1.0000
GTRIG	21	-3381.8269	6906.3872	0.0000	1.0000
HKY	15	-3412.2968	6926.5460	0.0000	1.0000
K2PI	13	-3426.1085	6940.5758	0.0000	1.0000
K2PG	13	-3426.3008	6940.9603	0.0000	1.0000
K2PIG	14	-3425.7349	6946.6253	0.0000	1.0000
GTR	19	-3410.4858	6950.1113	0.0000	1.0000
SYMG	17	-3419.8154	6955.1768	0.0000	1.0000
SYMI	17	-3420.1759	6955.8979	0.0000	1.0000
SYMIG	18	-3419.4867	6961.3162	0.0000	1.0000
K2P	12	-3457.9694	6997.5007	0.0000	1.0000
SYM	16	-3453.4528	7015.6548	0.0000	1.0000
F81I	15	-3518.9191	7139.7906	0.0000	1.0000
F81G	15	-3520.5854	7143.1232	0.0000	1.0000
F81IG	16	-3518.9213	7146.5918	0.0000	1.0000
F81	14	-3544.2518	7183.6592	0.0000	1.0000
JC69I	12	-3553.7086	7188.9790	0.0000	1.0000
JC69G	12	-3555.1537	7191.8693	0.0000	1.0000
JC69IG	13	-3553.7078	7195.7742	0.0000	1.0000
JC69	11	-3580.8515	7236.4681	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_4233.phylip = 1 - 895;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4233.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_4233.phylip = 1 - 895;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4233.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_4233.phylip = 1 - 895;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4233.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

