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

Input data from file "683-tomodel-phylip/chr1_32290.phylip" (7 taxa, 1084 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	-3069.3156	6178.6312	0.4372	0.4372
GTRG	20	-3069.9182	6179.8364	0.2393	0.6765
GTRIG	21	-3069.3265	6180.6530	0.1591	0.8356
HKYI	16	-3074.9690	6181.9381	0.0837	0.9192
HKYG	16	-3075.4785	6182.9571	0.0503	0.9695
HKYIG	17	-3074.9783	6183.9566	0.0305	1.0000
GTR	19	-3089.6625	6217.3250	0.0000	1.0000
HKY	15	-3095.6106	6221.2213	0.0000	1.0000
K2PG	13	-3149.5673	6325.1347	0.0000	1.0000
K2PI	13	-3149.5995	6325.1989	0.0000	1.0000
SYMG	17	-3146.2780	6326.5560	0.0000	1.0000
SYMI	17	-3146.3440	6326.6879	0.0000	1.0000
K2PIG	14	-3149.6168	6327.2337	0.0000	1.0000
SYMIG	18	-3146.3328	6328.6655	0.0000	1.0000
F81I	15	-3158.6800	6347.3600	0.0000	1.0000
F81G	15	-3159.0291	6348.0582	0.0000	1.0000
F81IG	16	-3158.6850	6349.3700	0.0000	1.0000
K2P	12	-3172.3274	6368.6547	0.0000	1.0000
SYM	16	-3169.2544	6370.5088	0.0000	1.0000
F81	14	-3177.7709	6383.5419	0.0000	1.0000
JC69I	12	-3223.0565	6470.1130	0.0000	1.0000
JC69G	12	-3223.1852	6470.3705	0.0000	1.0000
JC69IG	13	-3223.2211	6472.4421	0.0000	1.0000
JC69	11	-3243.4085	6508.8171	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3074.9690	6261.7527	0.6130	0.6130
HKYG	16	-3075.4785	6262.7717	0.3683	0.9813
HKYIG	17	-3074.9783	6268.7596	0.0184	0.9998
GTRI	20	-3069.3156	6278.3994	0.0001	0.9999
GTRG	20	-3069.9182	6279.6047	0.0001	1.0000
GTRIG	21	-3069.3265	6285.4097	0.0000	1.0000
HKY	15	-3095.6106	6296.0475	0.0000	1.0000
GTR	19	-3089.6625	6312.1048	0.0000	1.0000
K2PG	13	-3149.5673	6389.9840	0.0000	1.0000
K2PI	13	-3149.5995	6390.0483	0.0000	1.0000
K2PIG	14	-3149.6168	6397.0715	0.0000	1.0000
SYMG	17	-3146.2780	6411.3590	0.0000	1.0000
SYMI	17	-3146.3440	6411.4909	0.0000	1.0000
SYMIG	18	-3146.3328	6418.4570	0.0000	1.0000
F81I	15	-3158.6800	6422.1862	0.0000	1.0000
F81G	15	-3159.0291	6422.8844	0.0000	1.0000
K2P	12	-3172.3274	6428.5157	0.0000	1.0000
F81IG	16	-3158.6850	6429.1846	0.0000	1.0000
SYM	16	-3169.2544	6450.3234	0.0000	1.0000
F81	14	-3177.7709	6453.3796	0.0000	1.0000
JC69I	12	-3223.0565	6529.9739	0.0000	1.0000
JC69G	12	-3223.1852	6530.2314	0.0000	1.0000
JC69IG	13	-3223.2211	6537.2915	0.0000	1.0000
JC69	11	-3243.4085	6563.6896	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_32290.phylip = 1 - 1084;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32290.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_32290.phylip = 1 - 1084;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32290.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_32290.phylip = 1 - 1084;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32290.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

