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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-843.8408	1849.5027	0.2857	0.2857
GTRI	64	-844.3130	1850.4471	0.1782	0.4639
GTRIG	65	-842.8859	1850.7922	0.1499	0.6138
HKYG	60	-850.7793	1850.8385	0.1465	0.7603
HKYIG	61	-849.3710	1851.1195	0.1273	0.8876
HKYI	60	-851.0446	1851.3691	0.1124	1.0000
K2PIG	58	-871.1168	1885.3923	0.0000	1.0000
K2PG	57	-872.8582	1885.8507	0.0000	1.0000
K2PI	57	-872.9953	1886.1250	0.0000	1.0000
F81IG	60	-871.4890	1892.2580	0.0000	1.0000
F81G	59	-873.1392	1892.4855	0.0000	1.0000
F81I	59	-873.2906	1892.7884	0.0000	1.0000
SYMIG	62	-868.7770	1893.0539	0.0000	1.0000
SYMG	61	-870.4893	1893.3560	0.0000	1.0000
SYMI	61	-870.6521	1893.6818	0.0000	1.0000
HKY	59	-874.7971	1895.8015	0.0000	1.0000
GTR	63	-868.8617	1896.3712	0.0000	1.0000
JC69IG	57	-892.0987	1924.3319	0.0000	1.0000
JC69G	56	-893.6530	1924.4398	0.0000	1.0000
JC69I	56	-893.8005	1924.7348	0.0000	1.0000
K2P	56	-897.5695	1932.2728	0.0000	1.0000
SYM	60	-894.4522	1938.1843	0.0000	1.0000
F81	58	-897.6847	1938.5280	0.0000	1.0000
JC69	55	-918.1290	1970.4149	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-843.8408	1815.6815	0.3633	0.3633
GTRIG	65	-842.8859	1815.7718	0.3472	0.7105
GTRI	64	-844.3130	1816.6260	0.2265	0.9371
HKYIG	61	-849.3710	1820.7420	0.0289	0.9660
HKYG	60	-850.7793	1821.5585	0.0192	0.9852
HKYI	60	-851.0446	1822.0891	0.0148	1.0000
K2PIG	58	-871.1168	1858.2336	0.0000	1.0000
K2PG	57	-872.8582	1859.7163	0.0000	1.0000
K2PI	57	-872.9953	1859.9906	0.0000	1.0000
SYMIG	62	-868.7770	1861.5539	0.0000	1.0000
F81IG	60	-871.4890	1862.9780	0.0000	1.0000
SYMG	61	-870.4893	1862.9785	0.0000	1.0000
SYMI	61	-870.6521	1863.3042	0.0000	1.0000
GTR	63	-868.8617	1863.7234	0.0000	1.0000
F81G	59	-873.1392	1864.2784	0.0000	1.0000
F81I	59	-873.2906	1864.5812	0.0000	1.0000
HKY	59	-874.7971	1867.5943	0.0000	1.0000
JC69IG	57	-892.0987	1898.1975	0.0000	1.0000
JC69G	56	-893.6530	1899.3060	0.0000	1.0000
JC69I	56	-893.8005	1899.6010	0.0000	1.0000
K2P	56	-897.5695	1907.1389	0.0000	1.0000
SYM	60	-894.4522	1908.9043	0.0000	1.0000
F81	58	-897.6847	1911.3693	0.0000	1.0000
JC69	55	-918.1290	1946.2580	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-850.7793	2045.9461	0.4956	0.4956
HKYI	60	-851.0446	2046.4767	0.3801	0.8757
HKYIG	61	-849.3710	2048.8694	0.1149	0.9906
GTRG	64	-843.8408	2055.0283	0.0053	0.9959
GTRI	64	-844.3130	2055.9727	0.0033	0.9992
GTRIG	65	-842.8859	2058.8584	0.0008	1.0000
K2PG	57	-872.8582	2072.8845	0.0000	1.0000
K2PI	57	-872.9953	2073.1588	0.0000	1.0000
K2PIG	58	-871.1168	2075.1416	0.0000	1.0000
F81G	59	-873.1392	2084.9261	0.0000	1.0000
F81I	59	-873.2906	2085.2290	0.0000	1.0000
F81IG	60	-871.4890	2087.3656	0.0000	1.0000
HKY	59	-874.7971	2088.2421	0.0000	1.0000
SYMG	61	-870.4893	2091.1059	0.0000	1.0000
SYMI	61	-870.6521	2091.4316	0.0000	1.0000
SYMIG	62	-868.7770	2093.4211	0.0000	1.0000
GTR	63	-868.8617	2099.3304	0.0000	1.0000
JC69G	56	-893.6530	2108.7344	0.0000	1.0000
JC69I	56	-893.8005	2109.0294	0.0000	1.0000
JC69IG	57	-892.0987	2111.3657	0.0000	1.0000
K2P	56	-897.5695	2116.5673	0.0000	1.0000
F81	58	-897.6847	2128.2773	0.0000	1.0000
SYM	60	-894.4522	2133.2919	0.0000	1.0000
JC69	55	-918.1290	2151.9466	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/chr3_226.phylip = 1 - 311;
 Partition Dummy = 1:183-phylip/chr3_226.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/chr3_226.phylip = 1 - 311;
 Partition Dummy = 1:183-phylip/chr3_226.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 (HKYG)]

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

