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

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

Minimum AIC  model: HKYG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	57	-1065.4768	2289.9332	0.7353	0.7353
HKYG	60	-1061.1650	2293.1633	0.1462	0.8815
HKYI	60	-1062.4661	2295.7656	0.0398	0.9213
K2PI	57	-1068.6373	2296.2542	0.0312	0.9525
K2PIG	58	-1066.9246	2296.7258	0.0246	0.9772
HKYIG	61	-1061.1006	2297.0963	0.0205	0.9976
GTRG	64	-1057.8425	2303.1137	0.0010	0.9986
SYMG	61	-1064.2197	2303.3345	0.0009	0.9995
GTRI	64	-1059.3838	2306.1961	0.0002	0.9998
GTRIG	65	-1057.8004	2307.3275	0.0001	0.9999
K2P	56	-1077.1980	2309.5311	0.0000	0.9999
SYMI	61	-1067.3866	2309.6683	0.0000	1.0000
SYMIG	62	-1065.5632	2310.1406	0.0000	1.0000
HKY	59	-1072.4128	2311.6532	0.0000	1.0000
GTR	63	-1069.2360	2321.6634	0.0000	1.0000
SYM	60	-1075.8797	2322.5928	0.0000	1.0000
JC69G	56	-1104.2571	2363.6493	0.0000	1.0000
F81G	59	-1098.4146	2363.6568	0.0000	1.0000
F81IG	60	-1098.3898	2367.6130	0.0000	1.0000
F81I	59	-1102.2631	2371.3537	0.0000	1.0000
JC69I	56	-1108.1898	2371.5147	0.0000	1.0000
JC69IG	57	-1106.4453	2371.8702	0.0000	1.0000
F81	58	-1109.5480	2381.9727	0.0000	1.0000
JC69	55	-1115.8313	2383.0049	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1061.1650	2242.3299	0.3572	0.3572
GTRG	64	-1057.8425	2243.6851	0.1814	0.5386
HKYIG	61	-1061.1006	2244.2012	0.1401	0.6787
HKYI	60	-1062.4661	2244.9323	0.0972	0.7759
K2PG	57	-1065.4768	2244.9536	0.0962	0.8721
GTRIG	65	-1057.8004	2245.6009	0.0696	0.9417
GTRI	64	-1059.3838	2246.7675	0.0388	0.9805
K2PIG	58	-1066.9246	2249.8491	0.0083	0.9889
SYMG	61	-1064.2197	2250.4394	0.0062	0.9950
K2PI	57	-1068.6373	2251.2746	0.0041	0.9991
SYMIG	62	-1065.5632	2255.1265	0.0006	0.9997
SYMI	61	-1067.3866	2256.7732	0.0003	1.0000
HKY	59	-1072.4128	2262.8257	0.0000	1.0000
GTR	63	-1069.2360	2264.4719	0.0000	1.0000
K2P	56	-1077.1980	2266.3959	0.0000	1.0000
SYM	60	-1075.8797	2271.7594	0.0000	1.0000
F81G	59	-1098.4146	2314.8292	0.0000	1.0000
F81IG	60	-1098.3898	2316.7797	0.0000	1.0000
JC69G	56	-1104.2571	2320.5141	0.0000	1.0000
F81I	59	-1102.2631	2322.5262	0.0000	1.0000
JC69IG	57	-1106.4453	2326.8906	0.0000	1.0000
JC69I	56	-1108.1898	2328.3795	0.0000	1.0000
F81	58	-1109.5480	2335.0960	0.0000	1.0000
JC69	55	-1115.8313	2341.6627	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1065.4768	2434.3652	0.9147	0.9147
K2PI	57	-1068.6373	2440.6862	0.0388	0.9535
HKYG	60	-1061.1650	2441.7105	0.0232	0.9767
K2PIG	58	-1066.9246	2442.5837	0.0150	0.9917
HKYI	60	-1062.4661	2444.3129	0.0063	0.9981
HKYIG	61	-1061.1006	2446.9048	0.0017	0.9998
K2P	56	-1077.1980	2452.4845	0.0001	0.9999
SYMG	61	-1064.2197	2453.1430	0.0001	1.0000
GTRG	64	-1057.8425	2456.3577	0.0000	1.0000
HKY	59	-1072.4128	2458.8832	0.0000	1.0000
GTRI	64	-1059.3838	2459.4402	0.0000	1.0000
SYMI	61	-1067.3866	2459.4768	0.0000	1.0000
SYMIG	62	-1065.5632	2461.1531	0.0000	1.0000
GTRIG	65	-1057.8004	2461.5965	0.0000	1.0000
SYM	60	-1075.8797	2471.1400	0.0000	1.0000
GTR	63	-1069.2360	2473.8215	0.0000	1.0000
JC69G	56	-1104.2571	2506.6027	0.0000	1.0000
F81G	59	-1098.4146	2510.8868	0.0000	1.0000
JC69I	56	-1108.1898	2514.4681	0.0000	1.0000
F81IG	60	-1098.3898	2516.1603	0.0000	1.0000
JC69IG	57	-1106.4453	2516.3022	0.0000	1.0000
F81I	59	-1102.2631	2518.5837	0.0000	1.0000
JC69	55	-1115.8313	2524.4282	0.0000	1.0000
F81	58	-1109.5480	2527.8306	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_7915.phylip = 1 - 205;
 Partition Dummy = 1:183-phylip/chr2_7915.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_7915.phylip = 1 - 205;
 Partition Dummy = 1:183-phylip/chr2_7915.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_7915.phylip = 1 - 205;
 Partition Dummy = 1:183-phylip/chr2_7915.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

