
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:32 04/04/2011

Input data from file "591-tomodel-phylip/chr1_30414.phylip" (7 taxa, 643 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2881.9167	5796.7025	0.3932	0.3932
HKYIG	17	-2881.3680	5797.7152	0.2369	0.6301
GTRIG	21	-2877.5237	5798.5353	0.1572	0.7874
GTRG	20	-2878.6022	5798.5548	0.1557	0.9431
HKYI	16	-2884.4437	5801.7565	0.0314	0.9745
GTRI	20	-2880.4109	5802.1724	0.0255	1.0000
HKY	15	-2930.2376	5891.2408	0.0000	1.0000
GTR	19	-2927.7416	5894.7030	0.0000	1.0000
SYMG	17	-2939.4774	5913.9341	0.0000	1.0000
K2PG	13	-2944.5917	5915.7621	0.0000	1.0000
SYMIG	18	-2939.4718	5916.0397	0.0000	1.0000
K2PIG	14	-2944.5379	5917.7446	0.0000	1.0000
K2PI	13	-2948.6991	5923.9768	0.0000	1.0000
SYMI	17	-2944.5715	5924.1222	0.0000	1.0000
F81G	15	-2954.2372	5939.2399	0.0000	1.0000
F81IG	16	-2954.1859	5941.2407	0.0000	1.0000
F81I	15	-2957.2462	5945.2580	0.0000	1.0000
SYM	16	-2982.9696	5998.8082	0.0000	1.0000
K2P	12	-2987.6681	5999.8314	0.0000	1.0000
F81	14	-2992.6480	6013.9648	0.0000	1.0000
JC69G	12	-2999.2410	6022.9772	0.0000	1.0000
JC69IG	13	-2999.1758	6024.9302	0.0000	1.0000
JC69I	12	-3002.4005	6029.2962	0.0000	1.0000
JC69	11	-3036.2098	6094.8379	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2881.9167	5795.8335	0.3511	0.3511
HKYIG	17	-2881.3680	5796.7360	0.2236	0.5747
GTRIG	21	-2877.5237	5797.0473	0.1914	0.7660
GTRG	20	-2878.6022	5797.2043	0.1769	0.9430
GTRI	20	-2880.4109	5800.8219	0.0290	0.9719
HKYI	16	-2884.4437	5800.8875	0.0281	1.0000
HKY	15	-2930.2376	5890.4752	0.0000	1.0000
GTR	19	-2927.7416	5893.4831	0.0000	1.0000
SYMG	17	-2939.4774	5912.9549	0.0000	1.0000
SYMIG	18	-2939.4718	5914.9436	0.0000	1.0000
K2PG	13	-2944.5917	5915.1834	0.0000	1.0000
K2PIG	14	-2944.5379	5917.0758	0.0000	1.0000
SYMI	17	-2944.5715	5923.1430	0.0000	1.0000
K2PI	13	-2948.6991	5923.3981	0.0000	1.0000
F81G	15	-2954.2372	5938.4744	0.0000	1.0000
F81IG	16	-2954.1859	5940.3717	0.0000	1.0000
F81I	15	-2957.2462	5944.4924	0.0000	1.0000
SYM	16	-2982.9696	5997.9392	0.0000	1.0000
K2P	12	-2987.6681	5999.3362	0.0000	1.0000
F81	14	-2992.6480	6013.2960	0.0000	1.0000
JC69G	12	-2999.2410	6022.4820	0.0000	1.0000
JC69IG	13	-2999.1758	6024.3515	0.0000	1.0000
JC69I	12	-3002.4005	6028.8009	0.0000	1.0000
JC69	11	-3036.2098	6094.4196	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2881.9167	5867.2918	0.8709	0.8709
HKYI	16	-2884.4437	5872.3458	0.0696	0.9405
HKYIG	17	-2881.3680	5872.6605	0.0595	0.9999
GTRG	20	-2878.6022	5886.5272	0.0001	1.0000
GTRI	20	-2880.4109	5890.1448	0.0000	1.0000
GTRIG	21	-2877.5237	5890.8364	0.0000	1.0000
HKY	15	-2930.2376	5957.4674	0.0000	1.0000
K2PG	13	-2944.5917	5973.2433	0.0000	1.0000
GTR	19	-2927.7416	5978.3399	0.0000	1.0000
K2PIG	14	-2944.5379	5979.6018	0.0000	1.0000
K2PI	13	-2948.6991	5981.4580	0.0000	1.0000
SYMG	17	-2939.4774	5988.8793	0.0000	1.0000
SYMIG	18	-2939.4718	5995.3342	0.0000	1.0000
SYMI	17	-2944.5715	5999.0675	0.0000	1.0000
F81G	15	-2954.2372	6005.4666	0.0000	1.0000
F81I	15	-2957.2462	6011.4846	0.0000	1.0000
F81IG	16	-2954.1859	6011.8300	0.0000	1.0000
K2P	12	-2987.6681	6052.9299	0.0000	1.0000
SYM	16	-2982.9696	6069.3975	0.0000	1.0000
F81	14	-2992.6480	6075.8220	0.0000	1.0000
JC69G	12	-2999.2410	6076.0757	0.0000	1.0000
JC69I	12	-3002.4005	6082.3947	0.0000	1.0000
JC69IG	13	-2999.1758	6082.4114	0.0000	1.0000
JC69	11	-3036.2098	6143.5472	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 591-tomodel-phylip/chr1_30414.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr1_30414.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 591-tomodel-phylip/chr1_30414.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr1_30414.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_30414.phylip = 1 - 643;
 Partition Dummy = 1:591-tomodel-phylip/chr1_30414.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

