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

Input data from file "591-tomodel-phylip/chr1_28805.phylip" (7 taxa, 1217 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	-4088.7629	8217.5259	0.6777	0.6777
GTRIG	21	-4088.7641	8219.5282	0.2490	0.9268
GTRI	20	-4090.9893	8221.9785	0.0731	0.9999
HKYG	16	-4102.0148	8236.0295	0.0001	1.0000
HKYIG	17	-4102.0265	8238.0531	0.0000	1.0000
HKYI	16	-4103.6209	8239.2418	0.0000	1.0000
GTR	19	-4117.8097	8273.6194	0.0000	1.0000
HKY	15	-4128.8611	8287.7222	0.0000	1.0000
SYMG	17	-4131.8384	8297.6768	0.0000	1.0000
SYMIG	18	-4131.8495	8299.6991	0.0000	1.0000
SYMI	17	-4133.0758	8300.1515	0.0000	1.0000
K2PG	13	-4137.2574	8300.5148	0.0000	1.0000
K2PIG	14	-4137.2733	8302.5465	0.0000	1.0000
K2PI	13	-4138.5124	8303.0247	0.0000	1.0000
SYM	16	-4157.8797	8347.7594	0.0000	1.0000
K2P	12	-4163.0544	8350.1087	0.0000	1.0000
F81G	15	-4200.5691	8431.1382	0.0000	1.0000
F81I	15	-4201.5576	8433.1153	0.0000	1.0000
F81IG	16	-4200.5818	8433.1635	0.0000	1.0000
F81	14	-4222.8978	8473.7956	0.0000	1.0000
JC69G	12	-4227.1447	8478.2894	0.0000	1.0000
JC69I	12	-4227.9833	8479.9665	0.0000	1.0000
JC69IG	13	-4227.1386	8480.2771	0.0000	1.0000
JC69	11	-4249.1831	8520.3662	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4102.0148	8317.6958	0.6003	0.6003
GTRG	20	-4088.7629	8319.6087	0.2307	0.8310
HKYI	16	-4103.6209	8320.9081	0.1205	0.9515
GTRI	20	-4090.9893	8324.0614	0.0249	0.9764
HKYIG	17	-4102.0265	8324.8235	0.0170	0.9934
GTRIG	21	-4088.7641	8326.7152	0.0066	1.0000
HKY	15	-4128.8611	8364.2844	0.0000	1.0000
K2PG	13	-4137.2574	8366.8687	0.0000	1.0000
K2PI	13	-4138.5124	8369.3786	0.0000	1.0000
GTR	19	-4117.8097	8370.5982	0.0000	1.0000
K2PIG	14	-4137.2733	8374.0045	0.0000	1.0000
SYMG	17	-4131.8384	8384.4472	0.0000	1.0000
SYMI	17	-4133.0758	8386.9220	0.0000	1.0000
SYMIG	18	-4131.8495	8391.5737	0.0000	1.0000
K2P	12	-4163.0544	8411.3584	0.0000	1.0000
SYM	16	-4157.8797	8429.4257	0.0000	1.0000
F81G	15	-4200.5691	8507.7004	0.0000	1.0000
F81I	15	-4201.5576	8509.6775	0.0000	1.0000
F81IG	16	-4200.5818	8514.8298	0.0000	1.0000
JC69G	12	-4227.1447	8539.5391	0.0000	1.0000
JC69I	12	-4227.9833	8541.2163	0.0000	1.0000
F81	14	-4222.8978	8545.2537	0.0000	1.0000
JC69IG	13	-4227.1386	8546.6310	0.0000	1.0000
JC69	11	-4249.1831	8576.5118	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/chr1_28805.phylip = 1 - 1217;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28805.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/chr1_28805.phylip = 1 - 1217;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28805.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 591-tomodel-phylip/chr1_28805.phylip = 1 - 1217;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28805.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

