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

Input data from file "183-phylip/chr2_26677.phylip" (29 taxa, 746 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	60	-2836.0365	5802.7592	0.5211	0.5211
HKYIG	61	-2835.2386	5803.5357	0.3534	0.8745
GTRG	64	-2833.3155	5806.8484	0.0674	0.9419
GTRIG	65	-2832.5181	5807.6539	0.0451	0.9870
HKYI	60	-2839.8188	5810.3238	0.0119	0.9988
GTRI	64	-2837.3846	5814.9865	0.0012	1.0000
K2PG	57	-2889.0639	5901.7382	0.0000	1.0000
K2PIG	58	-2888.3735	5902.7091	0.0000	1.0000
HKY	59	-2888.1074	5904.5354	0.0000	1.0000
GTR	63	-2884.6434	5907.1108	0.0000	1.0000
K2PI	57	-2892.5740	5908.7584	0.0000	1.0000
SYMG	61	-2890.5387	5914.1359	0.0000	1.0000
SYMIG	62	-2889.8461	5915.1300	0.0000	1.0000
SYMI	61	-2892.4773	5918.0131	0.0000	1.0000
F81G	59	-2913.0938	5954.5083	0.0000	1.0000
F81IG	60	-2913.1514	5956.9889	0.0000	1.0000
F81I	59	-2917.3664	5963.0535	0.0000	1.0000
K2P	56	-2944.8092	6010.8841	0.0000	1.0000
SYM	60	-2945.6325	6021.9511	0.0000	1.0000
JC69G	56	-2955.5838	6032.4332	0.0000	1.0000
JC69IG	57	-2954.9052	6033.4208	0.0000	1.0000
JC69I	56	-2959.0348	6039.3353	0.0000	1.0000
F81	58	-2965.6691	6057.3003	0.0000	1.0000
JC69	55	-3010.7232	6140.3739	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-2836.0365	5792.0731	0.4255	0.4255
HKYIG	61	-2835.2386	5792.4772	0.3477	0.7732
GTRG	64	-2833.3155	5794.6310	0.1184	0.8916
GTRIG	65	-2832.5181	5795.0363	0.0967	0.9883
HKYI	60	-2839.8188	5799.6377	0.0097	0.9980
GTRI	64	-2837.3846	5802.7692	0.0020	1.0000
K2PG	57	-2889.0639	5892.1277	0.0000	1.0000
K2PIG	58	-2888.3735	5892.7469	0.0000	1.0000
HKY	59	-2888.1074	5894.2147	0.0000	1.0000
GTR	63	-2884.6434	5895.2867	0.0000	1.0000
K2PI	57	-2892.5740	5899.1479	0.0000	1.0000
SYMG	61	-2890.5387	5903.0775	0.0000	1.0000
SYMIG	62	-2889.8461	5903.6922	0.0000	1.0000
SYMI	61	-2892.4773	5906.9546	0.0000	1.0000
F81G	59	-2913.0938	5944.1876	0.0000	1.0000
F81IG	60	-2913.1514	5946.3028	0.0000	1.0000
F81I	59	-2917.3664	5952.7328	0.0000	1.0000
K2P	56	-2944.8092	6001.6185	0.0000	1.0000
SYM	60	-2945.6325	6011.2649	0.0000	1.0000
JC69G	56	-2955.5838	6023.1676	0.0000	1.0000
JC69IG	57	-2954.9052	6023.8104	0.0000	1.0000
JC69I	56	-2959.0348	6030.0697	0.0000	1.0000
F81	58	-2965.6691	6047.3381	0.0000	1.0000
JC69	55	-3010.7232	6131.4463	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2836.0365	6068.9566	0.9057	0.9057
HKYIG	61	-2835.2386	6073.9754	0.0736	0.9793
HKYI	60	-2839.8188	6076.5212	0.0206	1.0000
GTRG	64	-2833.3155	6089.9735	0.0000	1.0000
GTRIG	65	-2832.5181	6094.9934	0.0000	1.0000
GTRI	64	-2837.3846	6098.1117	0.0000	1.0000
K2PG	57	-2889.0639	6155.1671	0.0000	1.0000
K2PIG	58	-2888.3735	6160.4010	0.0000	1.0000
K2PI	57	-2892.5740	6162.1873	0.0000	1.0000
HKY	59	-2888.1074	6166.4835	0.0000	1.0000
SYMG	61	-2890.5387	6184.5757	0.0000	1.0000
GTR	63	-2884.6434	6186.0145	0.0000	1.0000
SYMI	61	-2892.4773	6188.4529	0.0000	1.0000
SYMIG	62	-2889.8461	6189.8052	0.0000	1.0000
F81G	59	-2913.0938	6216.4564	0.0000	1.0000
F81IG	60	-2913.1514	6223.1863	0.0000	1.0000
F81I	59	-2917.3664	6225.0017	0.0000	1.0000
K2P	56	-2944.8092	6260.0431	0.0000	1.0000
JC69G	56	-2955.5838	6281.5923	0.0000	1.0000
JC69IG	57	-2954.9052	6286.8497	0.0000	1.0000
SYM	60	-2945.6325	6288.1485	0.0000	1.0000
JC69I	56	-2959.0348	6288.4943	0.0000	1.0000
F81	58	-2965.6691	6314.9922	0.0000	1.0000
JC69	55	-3010.7232	6385.2562	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/chr2_26677.phylip = 1 - 746;
 Partition Dummy = 1:183-phylip/chr2_26677.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/chr2_26677.phylip = 1 - 746;
 Partition Dummy = 1:183-phylip/chr2_26677.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 183-phylip/chr2_26677.phylip = 1 - 746;
 Partition Dummy = 1:183-phylip/chr2_26677.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

