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

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

Minimum AIC  model: GTRG
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	-1080.5740	2311.5214	0.6316	0.6316
HKYIG	61	-1080.5083	2314.5333	0.1401	0.7717
HKYI	60	-1082.1316	2314.6367	0.1330	0.9047
GTRG	64	-1076.4426	2315.9907	0.0676	0.9723
GTRIG	65	-1076.3774	2319.1108	0.0142	0.9865
GTRI	64	-1078.1662	2319.4379	0.0121	0.9986
K2PG	57	-1091.7558	2324.6100	0.0009	0.9995
K2PI	57	-1093.1662	2327.4307	0.0002	0.9997
K2PIG	58	-1091.6890	2327.5427	0.0002	0.9999
HKY	59	-1091.8117	2330.8795	0.0000	1.0000
SYMG	61	-1089.4266	2332.3698	0.0000	1.0000
SYMI	61	-1090.9298	2335.3762	0.0000	1.0000
SYMIG	62	-1089.3728	2335.4318	0.0000	1.0000
GTR	63	-1088.0547	2335.9917	0.0000	1.0000
K2P	56	-1103.0564	2344.1699	0.0000	1.0000
SYM	60	-1100.4273	2351.2281	0.0000	1.0000
F81G	59	-1111.3323	2369.9208	0.0000	1.0000
F81I	59	-1112.6836	2372.6234	0.0000	1.0000
F81IG	60	-1111.2664	2372.9062	0.0000	1.0000
JC69G	56	-1119.7296	2377.5164	0.0000	1.0000
JC69I	56	-1121.1284	2380.3139	0.0000	1.0000
JC69IG	57	-1119.6632	2380.4247	0.0000	1.0000
F81	58	-1121.6302	2387.4249	0.0000	1.0000
JC69	55	-1130.5612	2396.1630	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1076.4426	2280.8852	0.3358	0.3358
HKYG	60	-1080.5740	2281.1480	0.2945	0.6303
GTRIG	65	-1076.3774	2282.7549	0.1319	0.7622
HKYIG	61	-1080.5083	2283.0167	0.1157	0.8779
HKYI	60	-1082.1316	2284.2633	0.0620	0.9399
GTRI	64	-1078.1662	2284.3324	0.0599	0.9998
K2PG	57	-1091.7558	2297.5116	0.0001	0.9999
K2PIG	58	-1091.6890	2299.3781	0.0000	0.9999
K2PI	57	-1093.1662	2300.3324	0.0000	1.0000
SYMG	61	-1089.4266	2300.8531	0.0000	1.0000
HKY	59	-1091.8117	2301.6233	0.0000	1.0000
GTR	63	-1088.0547	2302.1093	0.0000	1.0000
SYMIG	62	-1089.3728	2302.7456	0.0000	1.0000
SYMI	61	-1090.9298	2303.8595	0.0000	1.0000
K2P	56	-1103.0564	2318.1128	0.0000	1.0000
SYM	60	-1100.4273	2320.8547	0.0000	1.0000
F81G	59	-1111.3323	2340.6646	0.0000	1.0000
F81IG	60	-1111.2664	2342.5328	0.0000	1.0000
F81I	59	-1112.6836	2343.3672	0.0000	1.0000
JC69G	56	-1119.7296	2351.4593	0.0000	1.0000
JC69IG	57	-1119.6632	2353.3264	0.0000	1.0000
JC69I	56	-1121.1284	2354.2567	0.0000	1.0000
F81	58	-1121.6302	2359.2603	0.0000	1.0000
JC69	55	-1130.5612	2371.1224	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1080.5740	2503.7736	0.7307	0.7307
HKYI	60	-1082.1316	2506.8889	0.1539	0.8846
K2PG	57	-1091.7558	2509.0059	0.0534	0.9380
HKYIG	61	-1080.5083	2509.3527	0.0449	0.9829
K2PI	57	-1093.1662	2511.8267	0.0130	0.9959
K2PIG	58	-1091.6890	2514.5828	0.0033	0.9992
GTRG	64	-1076.4426	2518.3525	0.0005	0.9997
HKY	59	-1091.8117	2520.5385	0.0002	0.9999
GTRI	64	-1078.1662	2521.7997	0.0001	1.0000
GTRIG	65	-1076.3774	2523.9326	0.0000	1.0000
K2P	56	-1103.0564	2525.8967	0.0000	1.0000
SYMG	61	-1089.4266	2527.1892	0.0000	1.0000
SYMI	61	-1090.9298	2530.1956	0.0000	1.0000
SYMIG	62	-1089.3728	2532.7921	0.0000	1.0000
GTR	63	-1088.0547	2535.8662	0.0000	1.0000
SYM	60	-1100.4273	2543.4803	0.0000	1.0000
JC69G	56	-1119.7296	2559.2432	0.0000	1.0000
F81G	59	-1111.3323	2559.5798	0.0000	1.0000
JC69I	56	-1121.1284	2562.0406	0.0000	1.0000
F81I	59	-1112.6836	2562.2824	0.0000	1.0000
JC69IG	57	-1119.6632	2564.8207	0.0000	1.0000
F81IG	60	-1111.2664	2565.1584	0.0000	1.0000
F81	58	-1121.6302	2574.4651	0.0000	1.0000
JC69	55	-1130.5612	2575.1959	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/chr2_25816.phylip = 1 - 302;
 Partition Dummy = 1:183-phylip/chr2_25816.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 (HKYG)]

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

