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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2730.8178	5501.6357	0.4338	0.4338
GTRI	20	-2730.8655	5501.7311	0.4136	0.8473
GTRIG	21	-2730.8621	5503.7242	0.1527	1.0000
SYMI	17	-2751.2377	5536.4754	0.0000	1.0000
SYMG	17	-2751.4662	5536.9323	0.0000	1.0000
SYMIG	18	-2751.5008	5539.0017	0.0000	1.0000
HKYG	16	-2753.9020	5539.8041	0.0000	1.0000
HKYI	16	-2753.9363	5539.8727	0.0000	1.0000
HKYIG	17	-2753.9291	5541.8582	0.0000	1.0000
GTR	19	-2752.6545	5543.3089	0.0000	1.0000
K2PG	13	-2768.1969	5562.3939	0.0000	1.0000
K2PI	13	-2768.3270	5562.6540	0.0000	1.0000
K2PIG	14	-2768.2211	5564.4422	0.0000	1.0000
SYM	16	-2771.0634	5574.1269	0.0000	1.0000
HKY	15	-2773.1301	5576.2601	0.0000	1.0000
K2P	12	-2787.1228	5598.2456	0.0000	1.0000
F81I	15	-2802.4707	5634.9413	0.0000	1.0000
F81G	15	-2802.4900	5634.9801	0.0000	1.0000
F81IG	16	-2802.5046	5637.0092	0.0000	1.0000
JC69I	12	-2814.9785	5653.9571	0.0000	1.0000
JC69G	12	-2815.0858	5654.1717	0.0000	1.0000
JC69IG	13	-2815.0995	5656.1991	0.0000	1.0000
F81	14	-2820.2539	5668.5077	0.0000	1.0000
JC69	11	-2832.5944	5687.1888	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2730.8178	5599.3459	0.5040	0.5040
GTRI	20	-2730.8655	5599.4413	0.4805	0.9845
GTRIG	21	-2730.8621	5606.3199	0.0154	0.9999
HKYG	16	-2753.9020	5617.9722	0.0000	0.9999
HKYI	16	-2753.9363	5618.0408	0.0000	1.0000
SYMI	17	-2751.2377	5619.5291	0.0000	1.0000
SYMG	17	-2751.4662	5619.9860	0.0000	1.0000
HKYIG	17	-2753.9291	5624.9119	0.0000	1.0000
K2PG	13	-2768.1969	5625.9055	0.0000	1.0000
K2PI	13	-2768.3270	5626.1657	0.0000	1.0000
SYMIG	18	-2751.5008	5626.9409	0.0000	1.0000
K2PIG	14	-2768.2211	5632.8394	0.0000	1.0000
GTR	19	-2752.6545	5636.1336	0.0000	1.0000
HKY	15	-2773.1301	5649.5427	0.0000	1.0000
SYM	16	-2771.0634	5652.2950	0.0000	1.0000
K2P	12	-2787.1228	5656.8717	0.0000	1.0000
F81I	15	-2802.4707	5708.2240	0.0000	1.0000
F81G	15	-2802.4900	5708.2627	0.0000	1.0000
JC69I	12	-2814.9785	5712.5832	0.0000	1.0000
JC69G	12	-2815.0858	5712.7978	0.0000	1.0000
F81IG	16	-2802.5046	5715.1774	0.0000	1.0000
JC69IG	13	-2815.0995	5719.7107	0.0000	1.0000
F81	14	-2820.2539	5736.9049	0.0000	1.0000
JC69	11	-2832.5944	5740.9294	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 591-tomodel-phylip/chr11_5062.phylip = 1 - 978;
 Partition Dummy = 1:591-tomodel-phylip/chr11_5062.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 591-tomodel-phylip/chr11_5062.phylip = 1 - 978;
 Partition Dummy = 1:591-tomodel-phylip/chr11_5062.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 (GTRG)]

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


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

