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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: K2PI


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
HKYIG	61	-818.8278	1797.1011	0.4749	0.4749
K2PIG	58	-824.7263	1798.8380	0.1992	0.6741
HKYI	60	-822.1984	1800.4560	0.0887	0.7628
K2PI	57	-827.2377	1800.5725	0.0837	0.8465
HKYG	60	-822.3611	1800.7814	0.0754	0.9220
K2PG	57	-827.5844	1801.2659	0.0592	0.9811
GTRIG	65	-815.7132	1804.7598	0.0103	0.9914
SYMIG	62	-821.9041	1806.6739	0.0040	0.9954
SYMI	61	-824.6715	1808.7885	0.0014	0.9968
GTRI	64	-819.6386	1809.0863	0.0012	0.9980
GTRG	64	-819.7297	1809.2685	0.0011	0.9991
SYMG	61	-825.0485	1809.5425	0.0009	1.0000
JC69IG	57	-846.1818	1838.4607	0.0000	1.0000
F81IG	60	-841.4179	1838.8949	0.0000	1.0000
JC69I	56	-848.7318	1840.3042	0.0000	1.0000
JC69G	56	-849.0309	1840.9023	0.0000	1.0000
F81I	59	-844.1454	1840.9968	0.0000	1.0000
F81G	59	-844.4216	1841.5492	0.0000	1.0000
HKY	59	-860.2757	1873.2573	0.0000	1.0000
K2P	56	-865.7217	1874.2841	0.0000	1.0000
GTR	63	-858.3814	1883.0829	0.0000	1.0000
SYM	60	-863.8845	1883.8281	0.0000	1.0000
JC69	55	-887.6549	1914.9253	0.0000	1.0000
F81	58	-883.0766	1915.5385	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-818.8278	1759.6556	0.5796	0.5796
GTRIG	65	-815.7132	1761.4265	0.2391	0.8187
HKYI	60	-822.1984	1764.3969	0.0541	0.8728
HKYG	60	-822.3611	1764.7223	0.0460	0.9188
K2PIG	58	-824.7263	1765.4527	0.0319	0.9508
GTRI	64	-819.6386	1767.2773	0.0128	0.9636
GTRG	64	-819.7297	1767.4595	0.0117	0.9753
SYMIG	62	-821.9041	1767.8082	0.0098	0.9851
K2PI	57	-827.2377	1768.4755	0.0070	0.9922
K2PG	57	-827.5844	1769.1688	0.0050	0.9972
SYMI	61	-824.6715	1771.3430	0.0017	0.9988
SYMG	61	-825.0485	1772.0969	0.0012	1.0000
F81IG	60	-841.4179	1802.8358	0.0000	1.0000
F81I	59	-844.1454	1806.2909	0.0000	1.0000
JC69IG	57	-846.1818	1806.3637	0.0000	1.0000
F81G	59	-844.4216	1806.8433	0.0000	1.0000
JC69I	56	-848.7318	1809.4636	0.0000	1.0000
JC69G	56	-849.0309	1810.0617	0.0000	1.0000
HKY	59	-860.2757	1838.5514	0.0000	1.0000
GTR	63	-858.3814	1842.7629	0.0000	1.0000
K2P	56	-865.7217	1843.4435	0.0000	1.0000
SYM	60	-863.8845	1847.7690	0.0000	1.0000
F81	58	-883.0766	1882.1532	0.0000	1.0000
JC69	55	-887.6549	1885.3099	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-827.2377	1972.3046	0.3851	0.3851
K2PIG	58	-824.7263	1972.8577	0.2920	0.6771
K2PG	57	-827.5844	1972.9979	0.2723	0.9493
HKYIG	61	-818.8278	1977.7885	0.0248	0.9742
HKYI	60	-822.1984	1978.9538	0.0139	0.9880
HKYG	60	-822.3611	1979.2792	0.0118	0.9998
SYMI	61	-824.6715	1989.4759	0.0001	0.9999
SYMIG	62	-821.9041	1989.5170	0.0001	0.9999
SYMG	61	-825.0485	1990.2298	0.0000	1.0000
GTRIG	65	-815.7132	1993.8632	0.0000	1.0000
GTRI	64	-819.6386	1996.1380	0.0000	1.0000
GTRG	64	-819.7297	1996.3202	0.0000	1.0000
JC69I	56	-848.7318	2009.7167	0.0000	1.0000
JC69IG	57	-846.1818	2010.1928	0.0000	1.0000
JC69G	56	-849.0309	2010.3149	0.0000	1.0000
F81I	59	-844.1454	2017.2719	0.0000	1.0000
F81IG	60	-841.4179	2017.3927	0.0000	1.0000
F81G	59	-844.4216	2017.8243	0.0000	1.0000
K2P	56	-865.7217	2043.6966	0.0000	1.0000
HKY	59	-860.2757	2049.5324	0.0000	1.0000
SYM	60	-863.8845	2062.3259	0.0000	1.0000
GTR	63	-858.3814	2068.0477	0.0000	1.0000
JC69	55	-887.6549	2081.9871	0.0000	1.0000
F81	58	-883.0766	2089.5582	0.0000	1.0000



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

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

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


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

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


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_3180.phylip = 1 - 264;
 Partition Dummy = 1:183-phylip/chr3_3180.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

