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

Input data from file "591-tomodel-phylip/chr2_7957.phylip" (7 taxa, 1009 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	-2577.7215	5195.4430	0.5029	0.5029
GTRIG	21	-2577.2956	5196.5912	0.2833	0.7862
GTRI	20	-2578.6035	5197.2069	0.2082	0.9944
HKYG	16	-2586.9050	5205.8100	0.0028	0.9972
HKYIG	17	-2586.4822	5206.9644	0.0016	0.9988
HKYI	16	-2587.7279	5207.4557	0.0012	1.0000
SYMG	17	-2598.4060	5230.8120	0.0000	1.0000
SYMIG	18	-2598.0776	5232.1552	0.0000	1.0000
SYMI	17	-2599.3066	5232.6133	0.0000	1.0000
K2PG	13	-2604.8195	5235.6390	0.0000	1.0000
K2PIG	14	-2604.4690	5236.9380	0.0000	1.0000
K2PI	13	-2605.6409	5237.2818	0.0000	1.0000
GTR	19	-2606.3567	5250.7134	0.0000	1.0000
HKY	15	-2614.2349	5258.4698	0.0000	1.0000
SYM	16	-2623.7182	5279.4364	0.0000	1.0000
K2P	12	-2630.3683	5284.7366	0.0000	1.0000
F81G	15	-2635.2339	5300.4678	0.0000	1.0000
F81I	15	-2635.7452	5301.4903	0.0000	1.0000
F81IG	16	-2635.0639	5302.1278	0.0000	1.0000
JC69G	12	-2649.5284	5323.0567	0.0000	1.0000
JC69I	12	-2650.0479	5324.0958	0.0000	1.0000
JC69IG	13	-2649.3761	5324.7521	0.0000	1.0000
F81	14	-2659.6032	5347.2064	0.0000	1.0000
JC69	11	-2673.1930	5368.3861	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2586.9050	5284.4774	0.6658	0.6658
HKYI	16	-2587.7279	5286.1232	0.2924	0.9582
HKYIG	17	-2586.4822	5290.5486	0.0320	0.9902
GTRG	20	-2577.7215	5293.7773	0.0064	0.9965
GTRI	20	-2578.6035	5295.5412	0.0026	0.9992
K2PG	13	-2604.8195	5299.5563	0.0004	0.9995
GTRIG	21	-2577.2956	5299.8422	0.0003	0.9998
K2PI	13	-2605.6409	5301.1991	0.0002	1.0000
K2PIG	14	-2604.4690	5305.7720	0.0000	1.0000
SYMG	17	-2598.4060	5314.3961	0.0000	1.0000
SYMI	17	-2599.3066	5316.1975	0.0000	1.0000
SYMIG	18	-2598.0776	5320.6561	0.0000	1.0000
HKY	15	-2614.2349	5332.2205	0.0000	1.0000
K2P	12	-2630.3683	5343.7372	0.0000	1.0000
GTR	19	-2606.3567	5344.1310	0.0000	1.0000
SYM	16	-2623.7182	5358.1039	0.0000	1.0000
F81G	15	-2635.2339	5374.2185	0.0000	1.0000
F81I	15	-2635.7452	5375.2410	0.0000	1.0000
F81IG	16	-2635.0639	5380.7953	0.0000	1.0000
JC69G	12	-2649.5284	5382.0573	0.0000	1.0000
JC69I	12	-2650.0479	5383.0964	0.0000	1.0000
JC69IG	13	-2649.3761	5388.6694	0.0000	1.0000
F81	14	-2659.6032	5416.0404	0.0000	1.0000
JC69	11	-2673.1930	5422.4700	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/chr2_7957.phylip = 1 - 1009;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7957.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/chr2_7957.phylip = 1 - 1009;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7957.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/chr2_7957.phylip = 1 - 1009;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7957.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

