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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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	-1020.0089	2189.2979	0.4677	0.4677
HKYG	60	-1020.1358	2189.5517	0.4120	0.8797
HKYIG	61	-1020.1503	2192.6781	0.0863	0.9660
GTRI	64	-1017.1741	2196.1694	0.0151	0.9810
GTRG	64	-1017.2418	2196.3047	0.0141	0.9951
GTRIG	65	-1017.2437	2199.5078	0.0028	0.9979
HKY	59	-1027.0192	2200.2456	0.0020	0.9999
GTR	63	-1024.1613	2206.9704	0.0001	1.0000
K2PI	57	-1034.6648	2209.4639	0.0000	1.0000
K2PG	57	-1034.9341	2210.0025	0.0000	1.0000
K2PIG	58	-1034.6718	2212.5024	0.0000	1.0000
SYMI	61	-1031.7120	2215.8016	0.0000	1.0000
SYMG	61	-1032.0464	2216.4704	0.0000	1.0000
SYMIG	62	-1031.7219	2218.9437	0.0000	1.0000
K2P	56	-1041.8029	2220.7396	0.0000	1.0000
SYM	60	-1038.9073	2227.0946	0.0000	1.0000
F81I	59	-1045.4980	2237.2031	0.0000	1.0000
F81G	59	-1045.6889	2237.5851	0.0000	1.0000
F81IG	60	-1045.7018	2240.6836	0.0000	1.0000
F81	58	-1052.2508	2247.6603	0.0000	1.0000
JC69I	56	-1058.2682	2253.6702	0.0000	1.0000
JC69G	56	-1058.5728	2254.2794	0.0000	1.0000
JC69IG	57	-1058.2764	2256.6871	0.0000	1.0000
JC69	55	-1065.2209	2264.5986	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-1020.0089	2160.0179	0.3432	0.3432
HKYG	60	-1020.1358	2160.2717	0.3023	0.6455
HKYIG	61	-1020.1503	2162.3005	0.1096	0.7551
GTRI	64	-1017.1741	2162.3483	0.1070	0.8621
GTRG	64	-1017.2418	2162.4836	0.1000	0.9622
GTRIG	65	-1017.2437	2164.4874	0.0367	0.9989
HKY	59	-1027.0192	2172.0384	0.0008	0.9997
GTR	63	-1024.1613	2174.3227	0.0003	1.0000
K2PI	57	-1034.6648	2183.3295	0.0000	1.0000
K2PG	57	-1034.9341	2183.8682	0.0000	1.0000
K2PIG	58	-1034.6718	2185.3437	0.0000	1.0000
SYMI	61	-1031.7120	2185.4241	0.0000	1.0000
SYMG	61	-1032.0464	2186.0928	0.0000	1.0000
SYMIG	62	-1031.7219	2187.4437	0.0000	1.0000
K2P	56	-1041.8029	2195.6058	0.0000	1.0000
SYM	60	-1038.9073	2197.8146	0.0000	1.0000
F81I	59	-1045.4980	2208.9960	0.0000	1.0000
F81G	59	-1045.6889	2209.3779	0.0000	1.0000
F81IG	60	-1045.7018	2211.4036	0.0000	1.0000
F81	58	-1052.2508	2220.5016	0.0000	1.0000
JC69I	56	-1058.2682	2228.5364	0.0000	1.0000
JC69G	56	-1058.5728	2229.1456	0.0000	1.0000
JC69IG	57	-1058.2764	2230.5527	0.0000	1.0000
JC69	55	-1065.2209	2240.4418	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1020.0089	2384.4055	0.5126	0.5126
HKYG	60	-1020.1358	2384.6593	0.4515	0.9642
HKYIG	61	-1020.1503	2390.4279	0.0252	0.9894
HKY	59	-1027.0192	2392.6862	0.0082	0.9976
K2PI	57	-1034.6648	2396.4977	0.0012	0.9988
K2PG	57	-1034.9341	2397.0364	0.0009	0.9997
GTRI	64	-1017.1741	2401.6950	0.0001	0.9998
GTRG	64	-1017.2418	2401.8303	0.0001	0.9999
K2PIG	58	-1034.6718	2402.2517	0.0001	1.0000
K2P	56	-1041.8029	2405.0342	0.0000	1.0000
GTRIG	65	-1017.2437	2407.5739	0.0000	1.0000
GTR	63	-1024.1613	2409.9296	0.0000	1.0000
SYMI	61	-1031.7120	2413.5514	0.0000	1.0000
SYMG	61	-1032.0464	2414.2202	0.0000	1.0000
SYMIG	62	-1031.7219	2419.3109	0.0000	1.0000
SYM	60	-1038.9073	2422.2022	0.0000	1.0000
F81I	59	-1045.4980	2429.6437	0.0000	1.0000
F81G	59	-1045.6889	2430.0257	0.0000	1.0000
F81IG	60	-1045.7018	2435.7912	0.0000	1.0000
F81	58	-1052.2508	2437.4096	0.0000	1.0000
JC69I	56	-1058.2682	2437.9648	0.0000	1.0000
JC69G	56	-1058.5728	2438.5740	0.0000	1.0000
JC69IG	57	-1058.2764	2443.7209	0.0000	1.0000
JC69	55	-1065.2209	2446.1304	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/chrZ_6778.phylip = 1 - 311;
 Partition Dummy = 1:183-phylip/chrZ_6778.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/chrZ_6778.phylip = 1 - 311;
 Partition Dummy = 1:183-phylip/chrZ_6778.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 (HKYI)]

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


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

