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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2807.4622	5646.9244	0.2623	0.2623
GTRG	20	-2803.6758	5647.3517	0.2118	0.4741
HKYI	16	-2807.8645	5647.7290	0.1754	0.6496
GTRI	20	-2804.1475	5648.2950	0.1322	0.7817
HKYIG	17	-2807.2418	5648.4837	0.1203	0.9020
GTRIG	21	-2803.4469	5648.8939	0.0980	1.0000
F81G	15	-2831.3698	5692.7396	0.0000	1.0000
F81I	15	-2831.6929	5693.3859	0.0000	1.0000
F81IG	16	-2831.2327	5694.4655	0.0000	1.0000
GTR	19	-2832.3743	5702.7485	0.0000	1.0000
HKY	15	-2836.7991	5703.5982	0.0000	1.0000
SYMG	17	-2839.8331	5713.6663	0.0000	1.0000
SYMI	17	-2840.3735	5714.7471	0.0000	1.0000
SYMIG	18	-2839.5536	5715.1072	0.0000	1.0000
K2PG	13	-2844.8204	5715.6408	0.0000	1.0000
K2PI	13	-2845.2312	5716.4624	0.0000	1.0000
K2PIG	14	-2844.6003	5717.2006	0.0000	1.0000
F81	14	-2859.8026	5747.6053	0.0000	1.0000
JC69G	12	-2867.3603	5758.7206	0.0000	1.0000
JC69I	12	-2867.6914	5759.3828	0.0000	1.0000
JC69IG	13	-2867.2316	5760.4632	0.0000	1.0000
SYM	16	-2868.2243	5768.4486	0.0000	1.0000
K2P	12	-2873.3016	5770.6031	0.0000	1.0000
JC69	11	-2895.2722	5812.5444	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2807.4622	5727.9742	0.5864	0.5864
HKYI	16	-2807.8645	5728.7788	0.3922	0.9786
HKYIG	17	-2807.2418	5734.5991	0.0214	1.0000
GTRG	20	-2803.6758	5748.6639	0.0000	1.0000
GTRI	20	-2804.1475	5749.6073	0.0000	1.0000
GTRIG	21	-2803.4469	5755.2718	0.0000	1.0000
F81G	15	-2831.3698	5768.7238	0.0000	1.0000
F81I	15	-2831.6929	5769.3701	0.0000	1.0000
F81IG	16	-2831.2327	5775.5153	0.0000	1.0000
HKY	15	-2836.7991	5779.5824	0.0000	1.0000
K2PG	13	-2844.8204	5781.4937	0.0000	1.0000
K2PI	13	-2845.2312	5782.3154	0.0000	1.0000
K2PIG	14	-2844.6003	5788.1192	0.0000	1.0000
GTR	19	-2832.3743	5798.9952	0.0000	1.0000
SYMG	17	-2839.8331	5799.7817	0.0000	1.0000
SYMI	17	-2840.3735	5800.8625	0.0000	1.0000
SYMIG	18	-2839.5536	5806.2883	0.0000	1.0000
F81	14	-2859.8026	5818.5238	0.0000	1.0000
JC69G	12	-2867.3603	5819.5080	0.0000	1.0000
JC69I	12	-2867.6914	5820.1701	0.0000	1.0000
JC69IG	13	-2867.2316	5826.3162	0.0000	1.0000
K2P	12	-2873.3016	5831.3905	0.0000	1.0000
SYM	16	-2868.2243	5849.4984	0.0000	1.0000
JC69	11	-2895.2722	5868.2661	0.0000	1.0000



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

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

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


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

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


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

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

