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

Input data from file "183-phylip/chr1_28829.phylip" (29 taxa, 637 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-2139.2547	4423.5356	0.5151	0.5151
GTRG	64	-2140.5963	4423.7381	0.4655	0.9806
GTRI	64	-2143.7735	4430.0924	0.0194	1.0000
HKYG	60	-2163.0511	4458.8106	0.0000	1.0000
HKYIG	61	-2162.1993	4459.5533	0.0000	1.0000
HKYI	60	-2164.9263	4462.5609	0.0000	1.0000
SYMG	61	-2166.1828	4467.5203	0.0000	1.0000
SYMIG	62	-2165.5487	4468.7072	0.0000	1.0000
SYMI	61	-2167.8816	4470.9180	0.0000	1.0000
K2PG	57	-2182.4180	4490.2556	0.0000	1.0000
K2PIG	58	-2181.7245	4491.2898	0.0000	1.0000
K2PI	57	-2184.1609	4493.7416	0.0000	1.0000
F81G	59	-2188.9479	4508.1661	0.0000	1.0000
F81IG	60	-2188.2254	4509.1592	0.0000	1.0000
F81I	59	-2190.7978	4511.8659	0.0000	1.0000
GTR	63	-2190.3881	4520.8495	0.0000	1.0000
JC69G	56	-2206.1696	4535.3461	0.0000	1.0000
JC69IG	57	-2205.5451	4536.5099	0.0000	1.0000
JC69I	56	-2207.8421	4538.6912	0.0000	1.0000
HKY	59	-2209.0754	4548.4211	0.0000	1.0000
SYM	60	-2212.6843	4558.0770	0.0000	1.0000
K2P	56	-2228.6902	4580.3873	0.0000	1.0000
F81	58	-2235.0198	4597.8804	0.0000	1.0000
JC69	55	-2252.0755	4624.7535	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2139.2547	4408.5094	0.5746	0.5746
GTRG	64	-2140.5963	4409.1926	0.4083	0.9830
GTRI	64	-2143.7735	4415.5470	0.0170	1.0000
HKYG	60	-2163.0511	4446.1023	0.0000	1.0000
HKYIG	61	-2162.1993	4446.3985	0.0000	1.0000
HKYI	60	-2164.9263	4449.8525	0.0000	1.0000
SYMG	61	-2166.1828	4454.3656	0.0000	1.0000
SYMIG	62	-2165.5487	4455.0974	0.0000	1.0000
SYMI	61	-2167.8816	4457.7632	0.0000	1.0000
K2PG	57	-2182.4180	4478.8360	0.0000	1.0000
K2PIG	58	-2181.7245	4479.4489	0.0000	1.0000
K2PI	57	-2184.1609	4482.3219	0.0000	1.0000
F81G	59	-2188.9479	4495.8957	0.0000	1.0000
F81IG	60	-2188.2254	4496.4509	0.0000	1.0000
F81I	59	-2190.7978	4499.5956	0.0000	1.0000
GTR	63	-2190.3881	4506.7762	0.0000	1.0000
JC69G	56	-2206.1696	4524.3392	0.0000	1.0000
JC69IG	57	-2205.5451	4525.0902	0.0000	1.0000
JC69I	56	-2207.8421	4527.6843	0.0000	1.0000
HKY	59	-2209.0754	4536.1507	0.0000	1.0000
SYM	60	-2212.6843	4545.3686	0.0000	1.0000
K2P	56	-2228.6902	4569.3804	0.0000	1.0000
F81	58	-2235.0198	4586.0395	0.0000	1.0000
JC69	55	-2252.0755	4614.1511	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2140.5963	4694.4259	0.8380	0.8380
GTRIG	65	-2139.2547	4698.1994	0.1270	0.9650
GTRI	64	-2143.7735	4700.7802	0.0349	0.9999
HKYG	60	-2163.0511	4713.5085	0.0001	1.0000
HKYI	60	-2164.9263	4717.2587	0.0000	1.0000
HKYIG	61	-2162.1993	4718.2615	0.0000	1.0000
SYMG	61	-2166.1828	4726.2285	0.0000	1.0000
SYMI	61	-2167.8816	4729.6262	0.0000	1.0000
SYMIG	62	-2165.5487	4731.4171	0.0000	1.0000
K2PG	57	-2182.4180	4732.8718	0.0000	1.0000
K2PI	57	-2184.1609	4736.3578	0.0000	1.0000
K2PIG	58	-2181.7245	4737.9416	0.0000	1.0000
F81G	59	-2188.9479	4758.8451	0.0000	1.0000
F81I	59	-2190.7978	4762.5450	0.0000	1.0000
F81IG	60	-2188.2254	4763.8570	0.0000	1.0000
JC69G	56	-2206.1696	4773.9183	0.0000	1.0000
JC69I	56	-2207.8421	4777.2634	0.0000	1.0000
JC69IG	57	-2205.5451	4779.1261	0.0000	1.0000
GTR	63	-2190.3881	4787.5527	0.0000	1.0000
HKY	59	-2209.0754	4799.1001	0.0000	1.0000
SYM	60	-2212.6843	4812.7748	0.0000	1.0000
K2P	56	-2228.6902	4818.9595	0.0000	1.0000
F81	58	-2235.0198	4844.5322	0.0000	1.0000
JC69	55	-2252.0755	4859.2734	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_28829.phylip = 1 - 637;
 Partition Dummy = 1:183-phylip/chr1_28829.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_28829.phylip = 1 - 637;
 Partition Dummy = 1:183-phylip/chr1_28829.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_28829.phylip = 1 - 637;
 Partition Dummy = 1:183-phylip/chr1_28829.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

