
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:57 02/17/2011

Input data from file "183-phylip/chr7_5167.phylip" (29 taxa, 674 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	64	-2576.2824	5294.2266	0.8883	0.8883
GTRIG	65	-2577.1468	5298.4055	0.1099	0.9983
SYMG	61	-2586.6706	5307.7007	0.0011	0.9993
SYMIG	62	-2586.3473	5309.4801	0.0004	0.9997
GTRI	64	-2584.4173	5310.4964	0.0003	1.0000
HKYG	60	-2594.3385	5320.6184	0.0000	1.0000
SYMI	61	-2594.5549	5323.4694	0.0000	1.0000
K2PG	57	-2599.6048	5323.9434	0.0000	1.0000
HKYIG	61	-2595.1237	5324.6069	0.0000	1.0000
K2PIG	58	-2599.8659	5326.8603	0.0000	1.0000
HKYI	60	-2602.2539	5336.4492	0.0000	1.0000
K2PI	57	-2607.5612	5339.8562	0.0000	1.0000
GTR	63	-2624.6213	5388.4623	0.0000	1.0000
SYM	60	-2636.4629	5404.8670	0.0000	1.0000
HKY	59	-2647.7740	5425.0789	0.0000	1.0000
K2P	56	-2655.1004	5432.5476	0.0000	1.0000
F81G	59	-2662.4443	5454.4195	0.0000	1.0000
JC69G	56	-2667.2386	5456.8240	0.0000	1.0000
JC69IG	57	-2667.7297	5460.1932	0.0000	1.0000
F81IG	60	-2664.4129	5460.7672	0.0000	1.0000
F81I	59	-2670.4372	5470.4053	0.0000	1.0000
JC69I	56	-2675.3883	5473.1234	0.0000	1.0000
F81	58	-2715.0222	5557.1728	0.0000	1.0000
JC69	55	-2721.6375	5563.2426	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2576.2824	5280.5649	0.8649	0.8649
GTRIG	65	-2577.1468	5284.2937	0.1340	0.9989
SYMG	61	-2586.6706	5295.3412	0.0005	0.9995
SYMIG	62	-2586.3473	5296.6945	0.0003	0.9997
GTRI	64	-2584.4173	5296.8347	0.0003	1.0000
HKYG	60	-2594.3385	5308.6771	0.0000	1.0000
SYMI	61	-2594.5549	5311.1099	0.0000	1.0000
HKYIG	61	-2595.1237	5312.2474	0.0000	1.0000
K2PG	57	-2599.6048	5313.2097	0.0000	1.0000
K2PIG	58	-2599.8659	5315.7318	0.0000	1.0000
HKYI	60	-2602.2539	5324.5079	0.0000	1.0000
K2PI	57	-2607.5612	5329.1225	0.0000	1.0000
GTR	63	-2624.6213	5375.2426	0.0000	1.0000
SYM	60	-2636.4629	5392.9257	0.0000	1.0000
HKY	59	-2647.7740	5413.5480	0.0000	1.0000
K2P	56	-2655.1004	5422.2007	0.0000	1.0000
F81G	59	-2662.4443	5442.8886	0.0000	1.0000
JC69G	56	-2667.2386	5446.4772	0.0000	1.0000
F81IG	60	-2664.4129	5448.8259	0.0000	1.0000
JC69IG	57	-2667.7297	5449.4595	0.0000	1.0000
F81I	59	-2670.4372	5458.8743	0.0000	1.0000
JC69I	56	-2675.3883	5462.7766	0.0000	1.0000
F81	58	-2715.0222	5546.0443	0.0000	1.0000
JC69	55	-2721.6375	5553.2749	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2576.2824	5569.4116	0.4546	0.4546
K2PG	57	-2599.6048	5570.4638	0.2686	0.7233
SYMG	61	-2586.6706	5570.6483	0.2450	0.9683
SYMIG	62	-2586.3473	5576.5148	0.0130	0.9813
K2PIG	58	-2599.8659	5577.4992	0.0080	0.9893
GTRIG	65	-2577.1468	5577.6536	0.0074	0.9967
HKYG	60	-2594.3385	5579.4709	0.0030	0.9996
GTRI	64	-2584.4173	5585.6814	0.0001	0.9998
K2PI	57	-2607.5612	5586.3766	0.0001	0.9999
SYMI	61	-2594.5549	5586.4169	0.0001	0.9999
HKYIG	61	-2595.1237	5587.5545	0.0001	1.0000
HKYI	60	-2602.2539	5595.3017	0.0000	1.0000
GTR	63	-2624.6213	5659.5761	0.0000	1.0000
SYM	60	-2636.4629	5663.7195	0.0000	1.0000
K2P	56	-2655.1004	5674.9416	0.0000	1.0000
HKY	59	-2647.7740	5679.8285	0.0000	1.0000
JC69G	56	-2667.2386	5699.2180	0.0000	1.0000
JC69IG	57	-2667.7297	5706.7136	0.0000	1.0000
F81G	59	-2662.4443	5709.1692	0.0000	1.0000
JC69I	56	-2675.3883	5715.5175	0.0000	1.0000
F81IG	60	-2664.4129	5719.6197	0.0000	1.0000
F81I	59	-2670.4372	5725.1549	0.0000	1.0000
JC69	55	-2721.6375	5801.5026	0.0000	1.0000
F81	58	-2715.0222	5807.8117	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 183-phylip/chr7_5167.phylip = 1 - 674;
 Partition Dummy = 1:183-phylip/chr7_5167.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 183-phylip/chr7_5167.phylip = 1 - 674;
 Partition Dummy = 1:183-phylip/chr7_5167.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 183-phylip/chr7_5167.phylip = 1 - 674;
 Partition Dummy = 1:183-phylip/chr7_5167.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

