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

Input data from file "183-phylip/chr4_4009.phylip" (29 taxa, 343 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	-1035.6424	2217.2423	0.5718	0.5718
HKYI	60	-1036.5064	2218.9702	0.2410	0.8128
HKYIG	61	-1035.5574	2220.0330	0.1417	0.9545
GTRG	64	-1032.6483	2223.2246	0.0287	0.9832
GTRI	64	-1033.7189	2225.3658	0.0098	0.9930
GTRIG	65	-1032.5430	2226.0606	0.0070	1.0000
HKY	59	-1048.8899	2240.7974	0.0000	1.0000
K2PG	57	-1053.3638	2243.9276	0.0000	1.0000
K2PI	57	-1054.4329	2246.0659	0.0000	1.0000
K2PIG	58	-1053.1306	2246.3599	0.0000	1.0000
GTR	63	-1045.9483	2246.7998	0.0000	1.0000
SYMG	61	-1052.3294	2253.5769	0.0000	1.0000
SYMI	61	-1053.3323	2255.5828	0.0000	1.0000
SYMIG	62	-1052.1314	2256.1629	0.0000	1.0000
F81G	59	-1058.3087	2259.6351	0.0000	1.0000
F81I	59	-1059.3858	2261.7892	0.0000	1.0000
F81IG	60	-1058.1176	2262.1927	0.0000	1.0000
K2P	56	-1068.2614	2270.8446	0.0000	1.0000
SYM	60	-1066.9374	2279.8322	0.0000	1.0000
F81	58	-1072.1066	2284.3118	0.0000	1.0000
JC69G	56	-1076.1842	2286.6901	0.0000	1.0000
JC69I	56	-1077.2458	2288.8133	0.0000	1.0000
JC69IG	57	-1075.9542	2289.1083	0.0000	1.0000
JC69	55	-1090.6931	2312.8497	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1035.6424	2191.2848	0.4061	0.4061
HKYI	60	-1036.5064	2193.0128	0.1712	0.5772
HKYIG	61	-1035.5574	2193.1148	0.1626	0.7399
GTRG	64	-1032.6483	2193.2965	0.1485	0.8884
GTRIG	65	-1032.5430	2195.0859	0.0607	0.9491
GTRI	64	-1033.7189	2195.4377	0.0509	1.0000
HKY	59	-1048.8899	2215.7797	0.0000	1.0000
GTR	63	-1045.9483	2217.8966	0.0000	1.0000
K2PG	57	-1053.3638	2220.7276	0.0000	1.0000
K2PIG	58	-1053.1306	2222.2613	0.0000	1.0000
K2PI	57	-1054.4329	2222.8659	0.0000	1.0000
SYMG	61	-1052.3294	2226.6587	0.0000	1.0000
SYMIG	62	-1052.1314	2228.2629	0.0000	1.0000
SYMI	61	-1053.3323	2228.6646	0.0000	1.0000
F81G	59	-1058.3087	2234.6174	0.0000	1.0000
F81IG	60	-1058.1176	2236.2353	0.0000	1.0000
F81I	59	-1059.3858	2236.7715	0.0000	1.0000
K2P	56	-1068.2614	2248.5229	0.0000	1.0000
SYM	60	-1066.9374	2253.8748	0.0000	1.0000
F81	58	-1072.1066	2260.2132	0.0000	1.0000
JC69G	56	-1076.1842	2264.3684	0.0000	1.0000
JC69IG	57	-1075.9542	2265.9083	0.0000	1.0000
JC69I	56	-1077.2458	2266.4916	0.0000	1.0000
JC69	55	-1090.6931	2291.3863	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1035.6424	2421.5486	0.6753	0.6753
HKYI	60	-1036.5064	2423.2766	0.2846	0.9600
HKYIG	61	-1035.5574	2427.2164	0.0397	0.9997
GTRG	64	-1032.6483	2438.9113	0.0001	0.9998
K2PG	57	-1053.3638	2439.4782	0.0001	0.9999
GTRI	64	-1033.7189	2441.0525	0.0000	0.9999
K2PI	57	-1054.4329	2441.6165	0.0000	1.0000
HKY	59	-1048.8899	2442.2058	0.0000	1.0000
GTRIG	65	-1032.5430	2444.5384	0.0000	1.0000
K2PIG	58	-1053.1306	2444.8496	0.0000	1.0000
GTR	63	-1045.9483	2459.6736	0.0000	1.0000
SYMG	61	-1052.3294	2460.7603	0.0000	1.0000
F81G	59	-1058.3087	2461.0435	0.0000	1.0000
SYMI	61	-1053.3323	2462.7662	0.0000	1.0000
F81I	59	-1059.3858	2463.1976	0.0000	1.0000
K2P	56	-1068.2614	2463.4358	0.0000	1.0000
SYMIG	62	-1052.1314	2466.2022	0.0000	1.0000
F81IG	60	-1058.1176	2466.4991	0.0000	1.0000
JC69G	56	-1076.1842	2479.2813	0.0000	1.0000
JC69I	56	-1077.2458	2481.4045	0.0000	1.0000
F81	58	-1072.1066	2482.8016	0.0000	1.0000
SYM	60	-1066.9374	2484.1386	0.0000	1.0000
JC69IG	57	-1075.9542	2484.6590	0.0000	1.0000
JC69	55	-1090.6931	2502.4615	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/chr4_4009.phylip = 1 - 343;
 Partition Dummy = 1:183-phylip/chr4_4009.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/chr4_4009.phylip = 1 - 343;
 Partition Dummy = 1:183-phylip/chr4_4009.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/chr4_4009.phylip = 1 - 343;
 Partition Dummy = 1:183-phylip/chr4_4009.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

