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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-2584.3525	5209.8930	0.5557	0.5557
GTRI	20	-2585.2338	5211.6557	0.2302	0.7858
GTRIG	21	-2584.2455	5211.7998	0.2142	1.0000
HKYG	16	-2604.8445	5242.4541	0.0000	1.0000
HKYIG	17	-2604.8026	5244.4671	0.0000	1.0000
HKYI	16	-2606.3515	5245.4682	0.0000	1.0000
GTR	19	-2623.3444	5285.7622	0.0000	1.0000
SYMG	17	-2633.7540	5302.3700	0.0000	1.0000
SYMIG	18	-2633.7569	5304.4785	0.0000	1.0000
HKY	15	-2637.5897	5305.8535	0.0000	1.0000
SYMI	17	-2635.6808	5306.2236	0.0000	1.0000
K2PG	13	-2640.7950	5308.0998	0.0000	1.0000
K2PIG	14	-2640.7804	5310.1499	0.0000	1.0000
K2PI	13	-2643.4480	5313.4058	0.0000	1.0000
SYM	16	-2666.1229	5365.0110	0.0000	1.0000
K2P	12	-2671.6204	5367.6773	0.0000	1.0000
F81G	15	-2680.0413	5390.7568	0.0000	1.0000
F81IG	16	-2680.0304	5392.8258	0.0000	1.0000
F81I	15	-2681.3545	5393.3832	0.0000	1.0000
JC69G	12	-2705.0805	5434.5974	0.0000	1.0000
JC69IG	13	-2705.0830	5436.6758	0.0000	1.0000
JC69I	12	-2706.6629	5437.7622	0.0000	1.0000
F81	14	-2707.0826	5442.7542	0.0000	1.0000
JC69	11	-2731.5758	5485.5203	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2584.3525	5208.7049	0.5484	0.5484
GTRI	20	-2585.2338	5210.4676	0.2271	0.7755
GTRIG	21	-2584.2455	5210.4910	0.2245	1.0000
HKYG	16	-2604.8445	5241.6890	0.0000	1.0000
HKYIG	17	-2604.8026	5243.6052	0.0000	1.0000
HKYI	16	-2606.3515	5244.7031	0.0000	1.0000
GTR	19	-2623.3444	5284.6887	0.0000	1.0000
SYMG	17	-2633.7540	5301.5080	0.0000	1.0000
SYMIG	18	-2633.7569	5303.5138	0.0000	1.0000
HKY	15	-2637.5897	5305.1794	0.0000	1.0000
SYMI	17	-2635.6808	5305.3616	0.0000	1.0000
K2PG	13	-2640.7950	5307.5900	0.0000	1.0000
K2PIG	14	-2640.7804	5309.5608	0.0000	1.0000
K2PI	13	-2643.4480	5312.8960	0.0000	1.0000
SYM	16	-2666.1229	5364.2458	0.0000	1.0000
K2P	12	-2671.6204	5367.2409	0.0000	1.0000
F81G	15	-2680.0413	5390.0827	0.0000	1.0000
F81IG	16	-2680.0304	5392.0607	0.0000	1.0000
F81I	15	-2681.3545	5392.7090	0.0000	1.0000
JC69G	12	-2705.0805	5434.1611	0.0000	1.0000
JC69IG	13	-2705.0830	5436.1660	0.0000	1.0000
JC69I	12	-2706.6629	5437.3258	0.0000	1.0000
F81	14	-2707.0826	5442.1651	0.0000	1.0000
JC69	11	-2731.5758	5485.1515	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2584.3525	5300.5109	0.6867	0.6867
GTRI	20	-2585.2338	5302.2736	0.2844	0.9711
GTRIG	21	-2584.2455	5306.8873	0.0283	0.9994
HKYG	16	-2604.8445	5315.1338	0.0005	0.9999
HKYI	16	-2606.3515	5318.1479	0.0001	1.0000
HKYIG	17	-2604.8026	5321.6403	0.0000	1.0000
K2PG	13	-2640.7950	5367.2639	0.0000	1.0000
GTR	19	-2623.3444	5371.9045	0.0000	1.0000
K2PI	13	-2643.4480	5372.5699	0.0000	1.0000
K2PIG	14	-2640.7804	5373.8250	0.0000	1.0000
HKY	15	-2637.5897	5374.0339	0.0000	1.0000
SYMG	17	-2633.7540	5379.5431	0.0000	1.0000
SYMI	17	-2635.6808	5383.3968	0.0000	1.0000
SYMIG	18	-2633.7569	5386.1392	0.0000	1.0000
K2P	12	-2671.6204	5422.3245	0.0000	1.0000
SYM	16	-2666.1229	5437.6907	0.0000	1.0000
F81G	15	-2680.0413	5458.9372	0.0000	1.0000
F81I	15	-2681.3545	5461.5636	0.0000	1.0000
F81IG	16	-2680.0304	5465.5055	0.0000	1.0000
JC69G	12	-2705.0805	5489.2447	0.0000	1.0000
JC69I	12	-2706.6629	5492.4095	0.0000	1.0000
JC69IG	13	-2705.0830	5495.8399	0.0000	1.0000
F81	14	-2707.0826	5506.4294	0.0000	1.0000
JC69	11	-2731.5758	5535.6449	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_5845.phylip = 1 - 728;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5845.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_5845.phylip = 1 - 728;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5845.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 (GTRG)]

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


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

