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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-4266.1785	8572.3570	0.6306	0.6306
GTRIG	21	-4266.1794	8574.3589	0.2317	0.8623
GTRG	20	-4267.7000	8575.4001	0.1377	1.0000
HKYI	16	-4283.4037	8598.8073	0.0000	1.0000
HKYIG	17	-4283.4068	8600.8135	0.0000	1.0000
HKYG	16	-4284.7338	8601.4676	0.0000	1.0000
GTR	19	-4284.1419	8606.2837	0.0000	1.0000
SYMI	17	-4289.0303	8612.0606	0.0000	1.0000
SYMIG	18	-4289.0335	8614.0671	0.0000	1.0000
SYMG	17	-4290.2960	8614.5920	0.0000	1.0000
HKY	15	-4300.5428	8631.0856	0.0000	1.0000
K2PI	13	-4304.5284	8635.0569	0.0000	1.0000
K2PIG	14	-4304.5307	8637.0614	0.0000	1.0000
K2PG	13	-4305.7353	8637.4707	0.0000	1.0000
SYM	16	-4305.2660	8642.5320	0.0000	1.0000
K2P	12	-4320.8983	8665.7966	0.0000	1.0000
F81I	15	-4403.5250	8837.0500	0.0000	1.0000
F81IG	16	-4403.5505	8839.1011	0.0000	1.0000
F81G	15	-4405.3303	8840.6606	0.0000	1.0000
JC69I	12	-4417.5535	8859.1071	0.0000	1.0000
JC69IG	13	-4417.5761	8861.1522	0.0000	1.0000
JC69G	12	-4418.9237	8861.8474	0.0000	1.0000
F81	14	-4417.6868	8863.3736	0.0000	1.0000
JC69	11	-4431.6836	8885.3672	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	20	-4266.1785	8671.6775	0.7709	0.7709
GTRG	20	-4267.7000	8674.7206	0.1683	0.9393
HKYI	16	-4283.4037	8678.2637	0.0286	0.9679
GTRIG	21	-4266.1794	8678.6454	0.0237	0.9915
HKYG	16	-4284.7338	8680.9240	0.0076	0.9991
HKYIG	17	-4283.4068	8685.2360	0.0009	1.0000
SYMI	17	-4289.0303	8696.4831	0.0000	1.0000
SYMG	17	-4290.2960	8699.0144	0.0000	1.0000
K2PI	13	-4304.5284	8699.6152	0.0000	1.0000
GTR	19	-4284.1419	8700.6382	0.0000	1.0000
K2PG	13	-4305.7353	8702.0290	0.0000	1.0000
SYMIG	18	-4289.0335	8703.4555	0.0000	1.0000
HKY	15	-4300.5428	8705.5760	0.0000	1.0000
K2PIG	14	-4304.5307	8706.5857	0.0000	1.0000
SYM	16	-4305.2660	8721.9884	0.0000	1.0000
K2P	12	-4320.8983	8725.3889	0.0000	1.0000
F81I	15	-4403.5250	8911.5404	0.0000	1.0000
F81G	15	-4405.3303	8915.1509	0.0000	1.0000
F81IG	16	-4403.5505	8918.5575	0.0000	1.0000
JC69I	12	-4417.5535	8918.6994	0.0000	1.0000
JC69G	12	-4418.9237	8921.4397	0.0000	1.0000
JC69IG	13	-4417.5761	8925.7105	0.0000	1.0000
F81	14	-4417.6868	8932.8979	0.0000	1.0000
JC69	11	-4431.6836	8939.9934	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr4_8796.phylip = 1 - 1060;
 Partition Dummy = 1:683-tomodel-phylip/chr4_8796.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 AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr4_8796.phylip = 1 - 1060;
 Partition Dummy = 1:683-tomodel-phylip/chr4_8796.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 (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr4_8796.phylip = 1 - 1060;
 Partition Dummy = 1:683-tomodel-phylip/chr4_8796.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;


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

