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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1033.9016	2223.8623	0.6975	0.6975
HKYIG	61	-1033.6748	2226.7951	0.1610	0.8585
K2PG	57	-1041.2859	2228.6689	0.0631	0.9216
HKYI	60	-1036.4351	2228.9292	0.0554	0.9769
K2PIG	58	-1041.0791	2231.5435	0.0150	0.9919
K2PI	57	-1043.7251	2233.5472	0.0055	0.9974
GTRG	64	-1033.1172	2236.0435	0.0016	0.9990
SYMG	61	-1039.6247	2238.6949	0.0004	0.9994
GTRIG	65	-1032.8956	2239.1245	0.0003	0.9997
GTRI	64	-1035.7401	2241.2893	0.0001	0.9999
SYMIG	62	-1039.4126	2241.6908	0.0001	1.0000
SYMI	61	-1041.9353	2243.3161	0.0000	1.0000
F81G	59	-1048.7477	2250.2013	0.0000	1.0000
HKY	59	-1049.5636	2251.8331	0.0000	1.0000
F81IG	60	-1048.5340	2253.1272	0.0000	1.0000
F81I	59	-1051.3530	2255.4120	0.0000	1.0000
JC69G	56	-1056.8072	2256.4550	0.0000	1.0000
K2P	56	-1057.0233	2256.8872	0.0000	1.0000
JC69IG	57	-1056.6029	2259.3029	0.0000	1.0000
JC69I	56	-1059.2541	2261.3488	0.0000	1.0000
GTR	63	-1048.5267	2263.3734	0.0000	1.0000
SYM	60	-1054.9187	2265.8965	0.0000	1.0000
F81	58	-1063.8966	2277.1786	0.0000	1.0000
JC69	55	-1072.0507	2283.7168	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1033.9016	2187.8031	0.6160	0.6160
HKYIG	61	-1033.6748	2189.3495	0.2843	0.9003
HKYI	60	-1036.4351	2192.8701	0.0489	0.9492
GTRG	64	-1033.1172	2194.2345	0.0247	0.9739
GTRIG	65	-1032.8956	2195.7911	0.0113	0.9852
K2PG	57	-1041.2859	2196.5719	0.0077	0.9929
K2PIG	58	-1041.0791	2198.1581	0.0035	0.9964
GTRI	64	-1035.7401	2199.4802	0.0018	0.9982
SYMG	61	-1039.6247	2201.2493	0.0007	0.9989
K2PI	57	-1043.7251	2201.4502	0.0007	0.9996
SYMIG	62	-1039.4126	2202.8251	0.0003	0.9999
SYMI	61	-1041.9353	2205.8705	0.0001	1.0000
F81G	59	-1048.7477	2215.4954	0.0000	1.0000
F81IG	60	-1048.5340	2217.0680	0.0000	1.0000
HKY	59	-1049.5636	2217.1272	0.0000	1.0000
F81I	59	-1051.3530	2220.7061	0.0000	1.0000
GTR	63	-1048.5267	2223.0534	0.0000	1.0000
JC69G	56	-1056.8072	2225.6144	0.0000	1.0000
K2P	56	-1057.0233	2226.0467	0.0000	1.0000
JC69IG	57	-1056.6029	2227.2058	0.0000	1.0000
SYM	60	-1054.9187	2229.8374	0.0000	1.0000
JC69I	56	-1059.2541	2230.5082	0.0000	1.0000
F81	58	-1063.8966	2243.7932	0.0000	1.0000
JC69	55	-1072.0507	2254.1014	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1041.2859	2400.4010	0.6261	0.6261
HKYG	60	-1033.9016	2402.3601	0.2351	0.8611
K2PI	57	-1043.7251	2405.2793	0.0546	0.9157
K2PIG	58	-1041.0791	2405.5632	0.0474	0.9631
HKYI	60	-1036.4351	2407.4271	0.0187	0.9818
HKYIG	61	-1033.6748	2407.4824	0.0182	0.9999
SYMG	61	-1039.6247	2419.3822	0.0000	1.0000
GTRG	64	-1033.1172	2423.0952	0.0000	1.0000
SYMI	61	-1041.9353	2424.0034	0.0000	1.0000
SYMIG	62	-1039.4126	2424.5339	0.0000	1.0000
JC69G	56	-1056.8072	2425.8675	0.0000	1.0000
K2P	56	-1057.0233	2426.2998	0.0000	1.0000
F81G	59	-1048.7477	2426.4764	0.0000	1.0000
HKY	59	-1049.5636	2428.1082	0.0000	1.0000
GTRIG	65	-1032.8956	2428.2278	0.0000	1.0000
GTRI	64	-1035.7401	2428.3410	0.0000	1.0000
JC69I	56	-1059.2541	2430.7614	0.0000	1.0000
JC69IG	57	-1056.6029	2431.0349	0.0000	1.0000
F81IG	60	-1048.5340	2431.6250	0.0000	1.0000
F81I	59	-1051.3530	2431.6871	0.0000	1.0000
SYM	60	-1054.9187	2444.3944	0.0000	1.0000
GTR	63	-1048.5267	2448.3382	0.0000	1.0000
JC69	55	-1072.0507	2450.7786	0.0000	1.0000
F81	58	-1063.8966	2451.1983	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/chr3_24775.phylip = 1 - 264;
 Partition Dummy = 1:183-phylip/chr3_24775.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/chr3_24775.phylip = 1 - 264;
 Partition Dummy = 1:183-phylip/chr3_24775.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 (K2PG)]

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

