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

Input data from file "683-tomodel-phylip/chr2_9158.phylip" (7 taxa, 1094 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3320.9695	6681.9390	0.6633	0.6633
GTRIG	21	-3320.9766	6683.9531	0.2423	0.9055
GTRI	20	-3322.9193	6685.8385	0.0944	0.9999
HKYG	16	-3334.6442	6701.2885	0.0000	1.0000
HKYIG	17	-3334.6451	6703.2902	0.0000	1.0000
HKYI	16	-3335.9516	6703.9033	0.0000	1.0000
GTR	19	-3337.2738	6712.5476	0.0000	1.0000
HKY	15	-3352.1137	6734.2275	0.0000	1.0000
SYMG	17	-3356.1215	6746.2429	0.0000	1.0000
SYMIG	18	-3356.1290	6748.2580	0.0000	1.0000
SYMI	17	-3357.6411	6749.2822	0.0000	1.0000
K2PG	13	-3364.7112	6755.4224	0.0000	1.0000
K2PIG	14	-3364.7123	6757.4246	0.0000	1.0000
K2PI	13	-3366.0278	6758.0556	0.0000	1.0000
SYM	16	-3372.8733	6777.7465	0.0000	1.0000
K2P	12	-3382.3142	6788.6284	0.0000	1.0000
F81G	15	-3399.6942	6829.3883	0.0000	1.0000
F81I	15	-3400.5093	6831.0185	0.0000	1.0000
F81IG	16	-3399.6942	6831.3885	0.0000	1.0000
F81	14	-3416.1325	6860.2649	0.0000	1.0000
JC69G	12	-3425.6863	6875.3727	0.0000	1.0000
JC69I	12	-3426.4525	6876.9049	0.0000	1.0000
JC69IG	13	-3425.6933	6877.3865	0.0000	1.0000
JC69	11	-3442.5861	6907.1721	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3334.6442	6781.2500	0.4648	0.4648
GTRG	20	-3320.9695	6781.8909	0.3373	0.8021
HKYI	16	-3335.9516	6783.8648	0.1257	0.9278
GTRI	20	-3322.9193	6785.7904	0.0480	0.9758
HKYIG	17	-3334.6451	6788.2494	0.0140	0.9899
GTRIG	21	-3320.9766	6788.9027	0.0101	1.0000
GTR	19	-3337.2738	6807.5019	0.0000	1.0000
HKY	15	-3352.1137	6809.1914	0.0000	1.0000
K2PG	13	-3364.7112	6820.3911	0.0000	1.0000
K2PI	13	-3366.0278	6823.0243	0.0000	1.0000
K2PIG	14	-3364.7123	6827.3910	0.0000	1.0000
SYMG	17	-3356.1215	6831.2021	0.0000	1.0000
SYMI	17	-3357.6411	6834.2413	0.0000	1.0000
SYMIG	18	-3356.1290	6838.2148	0.0000	1.0000
K2P	12	-3382.3142	6848.5996	0.0000	1.0000
SYM	16	-3372.8733	6857.7081	0.0000	1.0000
F81G	15	-3399.6942	6904.3523	0.0000	1.0000
F81I	15	-3400.5093	6905.9824	0.0000	1.0000
F81IG	16	-3399.6942	6911.3500	0.0000	1.0000
F81	14	-3416.1325	6930.2312	0.0000	1.0000
JC69G	12	-3425.6863	6935.3439	0.0000	1.0000
JC69I	12	-3426.4525	6936.8761	0.0000	1.0000
JC69IG	13	-3425.6933	6942.3553	0.0000	1.0000
JC69	11	-3442.5861	6962.1457	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_9158.phylip = 1 - 1094;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9158.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_9158.phylip = 1 - 1094;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9158.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_9158.phylip = 1 - 1094;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9158.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

