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

Input data from file "683-tomodel-phylip/chr3_13321.phylip" (7 taxa, 877 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	-2248.3954	4536.7909	0.3998	0.3998
GTRI	20	-2248.5963	4537.1925	0.3270	0.7268
GTRIG	21	-2248.3905	4538.7809	0.1478	0.8746
HKYG	16	-2254.3378	4540.6757	0.0573	0.9319
HKYI	16	-2254.5534	4541.1069	0.0462	0.9781
HKYIG	17	-2254.3358	4542.6715	0.0211	0.9992
GTR	19	-2255.7681	4549.5362	0.0007	0.9999
HKY	15	-2261.9166	4553.8333	0.0001	1.0000
K2PG	13	-2279.3646	4584.7292	0.0000	1.0000
K2PI	13	-2279.4831	4584.9661	0.0000	1.0000
SYMG	17	-2275.8600	4585.7200	0.0000	1.0000
SYMI	17	-2275.9050	4585.8099	0.0000	1.0000
K2PIG	14	-2279.3663	4586.7326	0.0000	1.0000
SYMIG	18	-2275.8598	4587.7195	0.0000	1.0000
SYM	16	-2282.2789	4596.5577	0.0000	1.0000
K2P	12	-2286.5103	4597.0207	0.0000	1.0000
F81G	15	-2305.7599	4641.5199	0.0000	1.0000
F81I	15	-2305.8180	4641.6361	0.0000	1.0000
F81IG	16	-2305.7608	4643.5217	0.0000	1.0000
F81	14	-2311.9163	4651.8325	0.0000	1.0000
JC69G	12	-2326.7731	4677.5461	0.0000	1.0000
JC69I	12	-2326.8098	4677.6197	0.0000	1.0000
JC69IG	13	-2326.7801	4679.5602	0.0000	1.0000
JC69	11	-2333.2767	4688.5533	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2254.3378	4617.0998	0.5388	0.5388
HKYI	16	-2254.5534	4617.5310	0.4343	0.9731
HKYIG	17	-2254.3358	4623.8721	0.0182	0.9913
HKY	15	-2261.9166	4625.4809	0.0082	0.9995
GTRG	20	-2248.3954	4632.3210	0.0003	0.9998
GTRI	20	-2248.5963	4632.7226	0.0002	1.0000
GTRIG	21	-2248.3905	4639.0876	0.0000	1.0000
GTR	19	-2255.7681	4640.2899	0.0000	1.0000
K2PG	13	-2279.3646	4646.8238	0.0000	1.0000
K2PI	13	-2279.4831	4647.0607	0.0000	1.0000
K2PIG	14	-2279.3663	4653.6037	0.0000	1.0000
K2P	12	-2286.5103	4654.3388	0.0000	1.0000
SYMG	17	-2275.8600	4666.9207	0.0000	1.0000
SYMI	17	-2275.9050	4667.0105	0.0000	1.0000
SYM	16	-2282.2789	4672.9818	0.0000	1.0000
SYMIG	18	-2275.8598	4673.6967	0.0000	1.0000
F81G	15	-2305.7599	4713.1675	0.0000	1.0000
F81I	15	-2305.8180	4713.2837	0.0000	1.0000
F81	14	-2311.9163	4718.7036	0.0000	1.0000
F81IG	16	-2305.7608	4719.9458	0.0000	1.0000
JC69G	12	-2326.7731	4734.8642	0.0000	1.0000
JC69I	12	-2326.8098	4734.9377	0.0000	1.0000
JC69	11	-2333.2767	4741.0949	0.0000	1.0000
JC69IG	13	-2326.7801	4741.6548	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 683-tomodel-phylip/chr3_13321.phylip = 1 - 877;
 Partition Dummy = 1:683-tomodel-phylip/chr3_13321.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 683-tomodel-phylip/chr3_13321.phylip = 1 - 877;
 Partition Dummy = 1:683-tomodel-phylip/chr3_13321.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 683-tomodel-phylip/chr3_13321.phylip = 1 - 877;
 Partition Dummy = 1:683-tomodel-phylip/chr3_13321.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

