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

Input data from file "591-tomodel-phylip/chr5_10164.phylip" (7 taxa, 706 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	16	-2660.4174	5353.6245	0.3414	0.3414
HKYI	16	-2660.9162	5354.6219	0.2073	0.5487
GTRG	20	-2656.9129	5355.0521	0.1672	0.7159
HKYIG	17	-2660.3733	5355.6361	0.1249	0.8407
GTRI	20	-2657.4367	5356.0997	0.0990	0.9397
GTRIG	21	-2656.8718	5357.0945	0.0602	0.9999
SYMG	17	-2669.3492	5373.5879	0.0000	1.0000
SYMI	17	-2669.6024	5374.0943	0.0000	1.0000
K2PG	13	-2673.9112	5374.3485	0.0000	1.0000
K2PI	13	-2674.1168	5374.7596	0.0000	1.0000
SYMIG	18	-2669.2387	5375.4730	0.0000	1.0000
K2PIG	14	-2673.7749	5376.1575	0.0000	1.0000
HKY	15	-2678.2616	5387.2190	0.0000	1.0000
GTR	19	-2674.6900	5388.4879	0.0000	1.0000
SYM	16	-2687.1337	5407.0569	0.0000	1.0000
K2P	12	-2691.5674	5407.5851	0.0000	1.0000
F81I	15	-2721.8055	5474.3067	0.0000	1.0000
F81G	15	-2722.6146	5475.9249	0.0000	1.0000
F81IG	16	-2721.8083	5476.4062	0.0000	1.0000
JC69I	12	-2733.2770	5491.0042	0.0000	1.0000
JC69G	12	-2734.1428	5492.7358	0.0000	1.0000
JC69IG	13	-2733.2821	5493.0902	0.0000	1.0000
F81	14	-2737.7609	5504.1296	0.0000	1.0000
JC69	11	-2749.1585	5520.6974	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2660.4174	5352.8349	0.3129	0.3129
GTRG	20	-2656.9129	5353.8259	0.1907	0.5036
HKYI	16	-2660.9162	5353.8324	0.1900	0.6936
HKYIG	17	-2660.3733	5354.7466	0.1203	0.8139
GTRI	20	-2657.4367	5354.8734	0.1129	0.9269
GTRIG	21	-2656.8718	5355.7436	0.0731	0.9999
SYMG	17	-2669.3492	5372.6984	0.0000	1.0000
SYMI	17	-2669.6024	5373.2048	0.0000	1.0000
K2PG	13	-2673.9112	5373.8225	0.0000	1.0000
K2PI	13	-2674.1168	5374.2336	0.0000	1.0000
SYMIG	18	-2669.2387	5374.4774	0.0000	1.0000
K2PIG	14	-2673.7749	5375.5497	0.0000	1.0000
HKY	15	-2678.2616	5386.5233	0.0000	1.0000
GTR	19	-2674.6900	5387.3800	0.0000	1.0000
SYM	16	-2687.1337	5406.2674	0.0000	1.0000
K2P	12	-2691.5674	5407.1349	0.0000	1.0000
F81I	15	-2721.8055	5473.6110	0.0000	1.0000
F81G	15	-2722.6146	5475.2292	0.0000	1.0000
F81IG	16	-2721.8083	5475.6166	0.0000	1.0000
JC69I	12	-2733.2770	5490.5540	0.0000	1.0000
JC69G	12	-2734.1428	5492.2855	0.0000	1.0000
JC69IG	13	-2733.2821	5492.5642	0.0000	1.0000
F81	14	-2737.7609	5503.5218	0.0000	1.0000
JC69	11	-2749.1585	5520.3170	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2660.4174	5425.7887	0.5900	0.5900
HKYI	16	-2660.9162	5426.7862	0.3583	0.9484
HKYIG	17	-2660.3733	5432.2600	0.0232	0.9716
K2PG	13	-2673.9112	5433.0975	0.0153	0.9868
K2PI	13	-2674.1168	5433.5086	0.0124	0.9993
K2PIG	14	-2673.7749	5439.3843	0.0007	0.9999
GTRG	20	-2656.9129	5445.0182	0.0000	1.0000
GTRI	20	-2657.4367	5446.0657	0.0000	1.0000
SYMG	17	-2669.3492	5450.2118	0.0000	1.0000
SYMI	17	-2669.6024	5450.7183	0.0000	1.0000
GTRIG	21	-2656.8718	5451.4955	0.0000	1.0000
HKY	15	-2678.2616	5454.9175	0.0000	1.0000
SYMIG	18	-2669.2387	5456.5505	0.0000	1.0000
K2P	12	-2691.5674	5461.8503	0.0000	1.0000
GTR	19	-2674.6900	5474.0127	0.0000	1.0000
SYM	16	-2687.1337	5479.2212	0.0000	1.0000
F81I	15	-2721.8055	5542.0053	0.0000	1.0000
F81G	15	-2722.6146	5543.6234	0.0000	1.0000
JC69I	12	-2733.2770	5545.2694	0.0000	1.0000
JC69G	12	-2734.1428	5547.0009	0.0000	1.0000
F81IG	16	-2721.8083	5548.5704	0.0000	1.0000
JC69IG	13	-2733.2821	5551.8392	0.0000	1.0000
F81	14	-2737.7609	5567.3564	0.0000	1.0000
JC69	11	-2749.1585	5570.4727	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10164.phylip = 1 - 706;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10164.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;


[Mrbayes block for the best AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10164.phylip = 1 - 706;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10164.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;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10164.phylip = 1 - 706;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10164.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

