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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-1849.4252	3841.1459	0.2281	0.2281
HKYG	60	-1854.4261	3841.3436	0.2066	0.4347
GTRIG	65	-1848.4233	3841.6143	0.1805	0.6152
HKYIG	61	-1853.6703	3842.2705	0.1300	0.7452
HKYI	60	-1854.9144	3842.3203	0.1268	0.8720
GTRI	64	-1850.0161	3842.3277	0.1263	0.9983
K2PG	57	-1863.8560	3852.9377	0.0006	0.9990
K2PIG	58	-1863.0499	3853.7393	0.0004	0.9994
K2PI	57	-1864.2994	3853.8247	0.0004	0.9998
SYMG	61	-1860.9578	3856.8455	0.0001	0.9999
SYMIG	62	-1860.0082	3857.3932	0.0001	0.9999
SYMI	61	-1861.5054	3857.9408	0.0001	1.0000
F81G	59	-1880.5215	3891.1044	0.0000	1.0000
F81I	59	-1880.9353	3891.9320	0.0000	1.0000
F81IG	60	-1879.9051	3892.3016	0.0000	1.0000
JC69G	56	-1888.7851	3900.3905	0.0000	1.0000
HKY	59	-1885.2033	3900.4678	0.0000	1.0000
JC69I	56	-1889.2128	3901.2460	0.0000	1.0000
JC69IG	57	-1888.0904	3901.4066	0.0000	1.0000
GTR	63	-1881.0176	3901.8671	0.0000	1.0000
K2P	56	-1894.4045	3911.6293	0.0000	1.0000
SYM	60	-1891.0928	3914.6771	0.0000	1.0000
F81	58	-1910.9433	3949.5260	0.0000	1.0000
JC69	55	-1918.9527	3958.3283	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1848.4233	3826.8466	0.2914	0.2914
GTRG	64	-1849.4252	3826.8504	0.2908	0.5822
GTRI	64	-1850.0161	3828.0322	0.1611	0.7432
HKYG	60	-1854.4261	3828.8521	0.1069	0.8501
HKYIG	61	-1853.6703	3829.3406	0.0837	0.9338
HKYI	60	-1854.9144	3829.8288	0.0656	0.9994
K2PG	57	-1863.8560	3841.7119	0.0002	0.9996
K2PIG	58	-1863.0499	3842.0998	0.0001	0.9997
K2PI	57	-1864.2994	3842.5989	0.0001	0.9999
SYMG	61	-1860.9578	3843.9156	0.0001	0.9999
SYMIG	62	-1860.0082	3844.0165	0.0001	1.0000
SYMI	61	-1861.5054	3845.0108	0.0000	1.0000
F81G	59	-1880.5215	3879.0431	0.0000	1.0000
F81IG	60	-1879.9051	3879.8101	0.0000	1.0000
F81I	59	-1880.9353	3879.8706	0.0000	1.0000
GTR	63	-1881.0176	3888.0352	0.0000	1.0000
HKY	59	-1885.2033	3888.4065	0.0000	1.0000
JC69G	56	-1888.7851	3889.5702	0.0000	1.0000
JC69IG	57	-1888.0904	3890.1808	0.0000	1.0000
JC69I	56	-1889.2128	3890.4256	0.0000	1.0000
K2P	56	-1894.4045	3900.8089	0.0000	1.0000
SYM	60	-1891.0928	3902.1856	0.0000	1.0000
F81	58	-1910.9433	3937.8866	0.0000	1.0000
JC69	55	-1918.9527	3947.9053	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1863.8560	4096.6357	0.3317	0.3317
HKYG	60	-1854.4261	4097.1929	0.2510	0.5827
K2PI	57	-1864.2994	4097.5226	0.2129	0.7956
HKYI	60	-1854.9144	4098.1696	0.1540	0.9496
K2PIG	58	-1863.0499	4101.4959	0.0292	0.9788
HKYIG	61	-1853.6703	4102.1537	0.0210	0.9998
GTRG	64	-1849.4252	4113.0806	0.0001	0.9999
GTRI	64	-1850.0161	4114.2624	0.0000	1.0000
SYMG	61	-1860.9578	4116.7287	0.0000	1.0000
GTRIG	65	-1848.4233	4117.5491	0.0000	1.0000
SYMI	61	-1861.5054	4117.8240	0.0000	1.0000
SYMIG	62	-1860.0082	4121.3020	0.0000	1.0000
JC69G	56	-1888.7851	4140.0216	0.0000	1.0000
JC69I	56	-1889.2128	4140.8770	0.0000	1.0000
F81G	59	-1880.5215	4142.9115	0.0000	1.0000
F81I	59	-1880.9353	4143.7391	0.0000	1.0000
JC69IG	57	-1888.0904	4145.1046	0.0000	1.0000
F81IG	60	-1879.9051	4148.1509	0.0000	1.0000
K2P	56	-1894.4045	4151.2603	0.0000	1.0000
HKY	59	-1885.2033	4152.2749	0.0000	1.0000
GTR	63	-1881.0176	4169.7930	0.0000	1.0000
SYM	60	-1891.0928	4170.5264	0.0000	1.0000
JC69	55	-1918.9527	4193.8844	0.0000	1.0000
F81	58	-1910.9433	4197.2826	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/chr7_4698.phylip = 1 - 647;
 Partition Dummy = 1:183-phylip/chr7_4698.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 (GTRG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr7_4698.phylip = 1 - 647;
 Partition Dummy = 1:183-phylip/chr7_4698.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr7_4698.phylip = 1 - 647;
 Partition Dummy = 1:183-phylip/chr7_4698.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

