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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1874.7919	3884.5226	0.5815	0.5815
HKYIG	61	-1874.0460	3885.5604	0.3461	0.9276
HKYI	60	-1877.9165	3890.7719	0.0256	0.9531
GTRG	64	-1872.9167	3890.9528	0.0233	0.9765
GTRIG	65	-1872.1694	3892.0296	0.0136	0.9901
K2PG	57	-1883.2393	3893.8904	0.0054	0.9955
K2PIG	58	-1882.4794	3894.8693	0.0033	0.9988
GTRI	64	-1876.3184	3897.7561	0.0008	0.9996
K2PI	57	-1886.0646	3899.5410	0.0003	0.9999
SYMG	61	-1882.5706	3902.6095	0.0001	1.0000
SYMIG	62	-1881.8213	3903.6509	0.0000	1.0000
SYMI	61	-1885.4801	3908.4285	0.0000	1.0000
HKY	59	-1901.7219	3935.8634	0.0000	1.0000
GTR	63	-1899.7774	3942.1133	0.0000	1.0000
K2P	56	-1910.6739	3946.2709	0.0000	1.0000
F81G	59	-1907.9709	3948.3614	0.0000	1.0000
F81IG	60	-1907.0076	3948.9539	0.0000	1.0000
F81I	59	-1910.8042	3954.0279	0.0000	1.0000
SYM	60	-1909.6841	3954.3071	0.0000	1.0000
JC69G	56	-1917.0464	3959.0159	0.0000	1.0000
JC69IG	57	-1916.1287	3959.6691	0.0000	1.0000
JC69I	56	-1919.8028	3964.5287	0.0000	1.0000
F81	58	-1935.2082	4000.3269	0.0000	1.0000
JC69	55	-1944.2768	4010.9980	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1874.7919	3869.5839	0.4892	0.4892
HKYIG	61	-1874.0460	3870.0921	0.3795	0.8687
GTRG	64	-1872.9167	3873.8335	0.0584	0.9271
GTRIG	65	-1872.1694	3874.3389	0.0454	0.9725
HKYI	60	-1877.9165	3875.8331	0.0215	0.9940
K2PG	57	-1883.2393	3880.4786	0.0021	0.9961
GTRI	64	-1876.3184	3880.6368	0.0019	0.9981
K2PIG	58	-1882.4794	3880.9587	0.0017	0.9997
K2PI	57	-1886.0646	3886.1292	0.0001	0.9999
SYMG	61	-1882.5706	3887.1412	0.0001	0.9999
SYMIG	62	-1881.8213	3887.6427	0.0001	1.0000
SYMI	61	-1885.4801	3892.9602	0.0000	1.0000
HKY	59	-1901.7219	3921.4439	0.0000	1.0000
GTR	63	-1899.7774	3925.5548	0.0000	1.0000
K2P	56	-1910.6739	3933.3478	0.0000	1.0000
F81G	59	-1907.9709	3933.9418	0.0000	1.0000
F81IG	60	-1907.0076	3934.0151	0.0000	1.0000
SYM	60	-1909.6841	3939.3683	0.0000	1.0000
F81I	59	-1910.8042	3939.6083	0.0000	1.0000
JC69G	56	-1917.0464	3946.0928	0.0000	1.0000
JC69IG	57	-1916.1287	3946.2573	0.0000	1.0000
JC69I	56	-1919.8028	3951.6056	0.0000	1.0000
F81	58	-1935.2082	3986.4164	0.0000	1.0000
JC69	55	-1944.2768	3998.5536	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1883.2393	4126.2475	0.6414	0.6414
HKYG	60	-1874.7919	4128.2879	0.2312	0.8726
K2PIG	58	-1882.4794	4131.0393	0.0584	0.9310
K2PI	57	-1886.0646	4131.8981	0.0380	0.9691
HKYIG	61	-1874.0460	4133.1079	0.0208	0.9898
HKYI	60	-1877.9165	4134.5372	0.0102	1.0000
GTRG	64	-1872.9167	4149.7845	0.0000	1.0000
SYMG	61	-1882.5706	4150.1570	0.0000	1.0000
GTRIG	65	-1872.1694	4154.6016	0.0000	1.0000
SYMIG	62	-1881.8213	4154.9702	0.0000	1.0000
SYMI	61	-1885.4801	4155.9760	0.0000	1.0000
GTRI	64	-1876.3184	4156.5878	0.0000	1.0000
K2P	56	-1910.6739	4174.8050	0.0000	1.0000
HKY	59	-1901.7219	4175.8362	0.0000	1.0000
JC69G	56	-1917.0464	4187.5499	0.0000	1.0000
F81G	59	-1907.9709	4188.3342	0.0000	1.0000
JC69IG	57	-1916.1287	4192.0262	0.0000	1.0000
F81IG	60	-1907.0076	4192.7192	0.0000	1.0000
JC69I	56	-1919.8028	4193.0628	0.0000	1.0000
F81I	59	-1910.8042	4194.0007	0.0000	1.0000
GTR	63	-1899.7774	4197.1941	0.0000	1.0000
SYM	60	-1909.6841	4198.0724	0.0000	1.0000
JC69	55	-1944.2768	4235.6990	0.0000	1.0000
F81	58	-1935.2082	4236.4970	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 183-phylip/chrZ_7924.phylip = 1 - 551;
 Partition Dummy = 1:183-phylip/chrZ_7924.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 183-phylip/chrZ_7924.phylip = 1 - 551;
 Partition Dummy = 1:183-phylip/chrZ_7924.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 (K2PG)]

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

