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

Input data from file "591-tomodel-phylip/chr2_7985.phylip" (7 taxa, 1157 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	-3798.6256	7637.2512	0.4766	0.4766
GTRG	20	-3799.0664	7638.1327	0.3067	0.7833
GTRIG	21	-3798.4138	7638.8275	0.2167	1.0000
HKYG	16	-3812.8236	7657.6472	0.0000	1.0000
HKYI	16	-3813.0309	7658.0617	0.0000	1.0000
HKYIG	17	-3812.6620	7659.3241	0.0000	1.0000
GTR	19	-3841.3991	7720.7982	0.0000	1.0000
HKY	15	-3853.0109	7736.0218	0.0000	1.0000
SYMG	17	-3853.1940	7740.3879	0.0000	1.0000
K2PG	13	-3857.3819	7740.7637	0.0000	1.0000
SYMI	17	-3853.9953	7741.9905	0.0000	1.0000
SYMIG	18	-3853.2037	7742.4073	0.0000	1.0000
K2PI	13	-3858.4014	7742.8028	0.0000	1.0000
K2PIG	14	-3857.4090	7742.8180	0.0000	1.0000
F81G	15	-3877.7831	7785.5663	0.0000	1.0000
F81I	15	-3878.7525	7787.5050	0.0000	1.0000
F81IG	16	-3877.8256	7787.6512	0.0000	1.0000
SYM	16	-3892.5813	7817.1626	0.0000	1.0000
K2P	12	-3897.2379	7818.4759	0.0000	1.0000
JC69G	12	-3913.4231	7850.8462	0.0000	1.0000
JC69IG	13	-3913.4568	7852.9136	0.0000	1.0000
JC69I	12	-3914.7206	7853.4412	0.0000	1.0000
F81	14	-3914.7799	7857.5599	0.0000	1.0000
JC69	11	-3950.7282	7923.4564	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	20	-3798.6256	7738.3230	0.2970	0.2970
HKYG	16	-3812.8236	7738.5046	0.2712	0.5682
HKYI	16	-3813.0309	7738.9191	0.2205	0.7887
GTRG	20	-3799.0664	7739.2044	0.1911	0.9798
GTRIG	21	-3798.4138	7744.9528	0.0108	0.9906
HKYIG	17	-3812.6620	7745.2350	0.0094	1.0000
K2PG	13	-3857.3819	7806.4603	0.0000	1.0000
K2PI	13	-3858.4014	7808.4994	0.0000	1.0000
HKY	15	-3853.0109	7811.8256	0.0000	1.0000
K2PIG	14	-3857.4090	7813.5682	0.0000	1.0000
GTR	19	-3841.3991	7816.8163	0.0000	1.0000
SYMG	17	-3853.1940	7826.2989	0.0000	1.0000
SYMI	17	-3853.9953	7827.9015	0.0000	1.0000
SYMIG	18	-3853.2037	7833.3719	0.0000	1.0000
F81G	15	-3877.7831	7861.3700	0.0000	1.0000
F81I	15	-3878.7525	7863.3088	0.0000	1.0000
F81IG	16	-3877.8256	7868.5086	0.0000	1.0000
K2P	12	-3897.2379	7879.1189	0.0000	1.0000
SYM	16	-3892.5813	7898.0200	0.0000	1.0000
JC69G	12	-3913.4231	7911.4892	0.0000	1.0000
JC69I	12	-3914.7206	7914.0842	0.0000	1.0000
JC69IG	13	-3913.4568	7918.6102	0.0000	1.0000
F81	14	-3914.7799	7928.3101	0.0000	1.0000
JC69	11	-3950.7282	7979.0459	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 591-tomodel-phylip/chr2_7985.phylip = 1 - 1157;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7985.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 591-tomodel-phylip/chr2_7985.phylip = 1 - 1157;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7985.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 591-tomodel-phylip/chr2_7985.phylip = 1 - 1157;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7985.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

