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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	60	-1066.4564	2279.9239	0.4131	0.4131
HKYG	60	-1066.6461	2280.3032	0.3417	0.7548
HKYIG	61	-1066.0202	2282.0551	0.1423	0.8971
GTRI	64	-1062.7750	2284.7111	0.0377	0.9349
GTRG	64	-1062.9066	2284.9742	0.0331	0.9679
GTRIG	65	-1062.3495	2286.9547	0.0123	0.9802
K2PG	57	-1074.8026	2287.7367	0.0083	0.9885
K2PI	57	-1074.8662	2287.8637	0.0078	0.9963
K2PIG	58	-1074.2146	2289.4989	0.0034	0.9998
SYMG	61	-1073.1395	2296.2938	0.0001	0.9999
SYMI	61	-1073.3708	2296.7564	0.0001	1.0000
SYMIG	62	-1072.7002	2298.4413	0.0000	1.0000
HKY	59	-1083.1659	2310.3612	0.0000	1.0000
GTR	63	-1079.0949	2314.2794	0.0000	1.0000
F81I	59	-1085.9379	2315.9052	0.0000	1.0000
F81G	59	-1086.6497	2317.3288	0.0000	1.0000
F81IG	60	-1085.4108	2317.8326	0.0000	1.0000
K2P	56	-1092.4765	2320.1676	0.0000	1.0000
JC69I	56	-1093.3853	2321.9851	0.0000	1.0000
JC69G	56	-1093.8831	2322.9808	0.0000	1.0000
JC69IG	57	-1092.7347	2323.6008	0.0000	1.0000
SYM	60	-1090.2147	2327.4405	0.0000	1.0000
F81	58	-1103.4280	2347.9255	0.0000	1.0000
JC69	55	-1111.2369	2354.7926	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-1066.4564	2252.9128	0.2390	0.2390
HKYG	60	-1066.6461	2253.2921	0.1977	0.4367
GTRI	64	-1062.7750	2253.5500	0.1738	0.6105
GTRG	64	-1062.9066	2253.8131	0.1524	0.7629
HKYIG	61	-1066.0202	2254.0403	0.1360	0.8989
GTRIG	65	-1062.3495	2254.6990	0.0978	0.9968
K2PG	57	-1074.8026	2263.6053	0.0011	0.9979
K2PI	57	-1074.8662	2263.7323	0.0011	0.9990
K2PIG	58	-1074.2146	2264.4293	0.0008	0.9997
SYMG	61	-1073.1395	2268.2790	0.0001	0.9998
SYMI	61	-1073.3708	2268.7416	0.0001	0.9999
SYMIG	62	-1072.7002	2269.4004	0.0001	1.0000
GTR	63	-1079.0949	2284.1898	0.0000	1.0000
HKY	59	-1083.1659	2284.3317	0.0000	1.0000
F81I	59	-1085.9379	2289.8757	0.0000	1.0000
F81IG	60	-1085.4108	2290.8216	0.0000	1.0000
F81G	59	-1086.6497	2291.2993	0.0000	1.0000
K2P	56	-1092.4765	2296.9531	0.0000	1.0000
JC69I	56	-1093.3853	2298.7706	0.0000	1.0000
JC69IG	57	-1092.7347	2299.4694	0.0000	1.0000
JC69G	56	-1093.8831	2299.7663	0.0000	1.0000
SYM	60	-1090.2147	2300.4294	0.0000	1.0000
F81	58	-1103.4280	2322.8559	0.0000	1.0000
JC69	55	-1111.2369	2332.4737	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1074.8026	2480.4980	0.2967	0.2967
K2PI	57	-1074.8662	2480.6250	0.2785	0.5752
HKYI	60	-1066.4564	2481.2209	0.2067	0.7819
HKYG	60	-1066.6461	2481.6002	0.1710	0.9530
K2PIG	58	-1074.2146	2485.1271	0.0293	0.9823
HKYIG	61	-1066.0202	2486.1535	0.0176	0.9998
GTRI	64	-1062.7750	2497.0787	0.0001	0.9999
GTRG	64	-1062.9066	2497.3418	0.0001	1.0000
SYMG	61	-1073.1395	2500.3922	0.0000	1.0000
SYMI	61	-1073.3708	2500.8549	0.0000	1.0000
GTRIG	65	-1062.3495	2502.0328	0.0000	1.0000
SYMIG	62	-1072.7002	2505.3187	0.0000	1.0000
HKY	59	-1083.1659	2508.8347	0.0000	1.0000
K2P	56	-1092.4765	2510.0406	0.0000	1.0000
JC69I	56	-1093.3853	2511.8582	0.0000	1.0000
JC69G	56	-1093.8831	2512.8539	0.0000	1.0000
F81I	59	-1085.9379	2514.3787	0.0000	1.0000
F81G	59	-1086.6497	2515.8023	0.0000	1.0000
JC69IG	57	-1092.7347	2516.3621	0.0000	1.0000
F81IG	60	-1085.4108	2519.1297	0.0000	1.0000
GTR	63	-1079.0949	2523.9133	0.0000	1.0000
SYM	60	-1090.2147	2528.7375	0.0000	1.0000
JC69	55	-1111.2369	2541.7562	0.0000	1.0000
F81	58	-1103.4280	2543.5537	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_9719.phylip = 1 - 332;
 Partition Dummy = 1:183-phylip/chr3_9719.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_9719.phylip = 1 - 332;
 Partition Dummy = 1:183-phylip/chr3_9719.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

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

