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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-4201.2986	8444.5971	0.3701	0.3701
GTRI	20	-4202.3254	8444.6508	0.3603	0.7304
GTRG	20	-4202.8287	8445.6573	0.2178	0.9482
HKYI	16	-4208.9586	8449.9172	0.0259	0.9741
HKYIG	17	-4208.3403	8450.6805	0.0177	0.9918
HKYG	16	-4210.1094	8452.2188	0.0082	1.0000
SYMG	17	-4232.0514	8498.1027	0.0000	1.0000
SYMIG	18	-4231.3555	8498.7109	0.0000	1.0000
K2PI	13	-4236.5459	8499.0919	0.0000	1.0000
SYMI	17	-4232.6967	8499.3934	0.0000	1.0000
K2PIG	14	-4235.8118	8499.6236	0.0000	1.0000
K2PG	13	-4237.0791	8500.1582	0.0000	1.0000
GTR	19	-4249.2342	8536.4684	0.0000	1.0000
HKY	15	-4253.7634	8537.5268	0.0000	1.0000
SYM	16	-4272.7332	8577.4664	0.0000	1.0000
K2P	12	-4276.8919	8577.7838	0.0000	1.0000
F81I	15	-4320.7908	8671.5817	0.0000	1.0000
F81IG	16	-4320.4132	8672.8264	0.0000	1.0000
F81G	15	-4321.8289	8673.6578	0.0000	1.0000
JC69I	12	-4338.3441	8700.6883	0.0000	1.0000
JC69IG	13	-4337.9491	8701.8982	0.0000	1.0000
JC69G	12	-4339.2481	8702.4962	0.0000	1.0000
F81	14	-4356.2666	8740.5331	0.0000	1.0000
JC69	11	-4372.2502	8766.5005	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-4208.9586	8528.2966	0.7264	0.7264
HKYG	16	-4210.1094	8530.5983	0.2298	0.9562
HKYIG	17	-4208.3403	8533.9587	0.0428	0.9990
GTRI	20	-4202.3254	8542.6251	0.0006	0.9996
GTRG	20	-4202.8287	8543.6316	0.0003	1.0000
GTRIG	21	-4201.2986	8547.4701	0.0000	1.0000
K2PI	13	-4236.5459	8562.7752	0.0000	1.0000
K2PG	13	-4237.0791	8563.8415	0.0000	1.0000
K2PIG	14	-4235.8118	8568.2056	0.0000	1.0000
SYMG	17	-4232.0514	8581.3809	0.0000	1.0000
SYMI	17	-4232.6967	8582.6715	0.0000	1.0000
SYMIG	18	-4231.3555	8586.8878	0.0000	1.0000
HKY	15	-4253.7634	8611.0075	0.0000	1.0000
GTR	19	-4249.2342	8629.5440	0.0000	1.0000
K2P	12	-4276.8919	8636.5684	0.0000	1.0000
SYM	16	-4272.7332	8655.8459	0.0000	1.0000
F81I	15	-4320.7908	8745.0624	0.0000	1.0000
F81G	15	-4321.8289	8747.1386	0.0000	1.0000
F81IG	16	-4320.4132	8751.2058	0.0000	1.0000
JC69I	12	-4338.3441	8759.4729	0.0000	1.0000
JC69G	12	-4339.2481	8761.2808	0.0000	1.0000
JC69IG	13	-4337.9491	8765.5815	0.0000	1.0000
F81	14	-4356.2666	8809.1151	0.0000	1.0000
JC69	11	-4372.2502	8820.3863	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_14682.phylip = 1 - 991;
 Partition Dummy = 1:591-tomodel-phylip/chr5_14682.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

