
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:30 08/01/2011

Input data from file "261-aggressive-phylip/4076_hg18.phylip" (29 taxa, 327 characters)

Minimum AIC  model: GTRI
Minimum AICc model: HKYI
Minimum BIC  model: K2PI


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	-1033.6026	2214.7240	0.3834	0.3834
HKYG	60	-1033.7270	2214.9728	0.3385	0.7219
GTRI	64	-1028.9450	2217.6458	0.0890	0.8108
HKYIG	61	-1033.7268	2217.9970	0.0746	0.8855
GTRG	64	-1029.2898	2218.3354	0.0630	0.9485
GTRIG	65	-1028.8190	2220.5116	0.0212	0.9697
K2PI	57	-1041.4088	2221.3975	0.0136	0.9833
K2PG	57	-1041.4340	2221.4480	0.0133	0.9966
K2PIG	58	-1041.5325	2224.6023	0.0027	0.9994
SYMI	61	-1039.1030	2228.7493	0.0003	0.9997
SYMG	61	-1039.6438	2229.8310	0.0002	0.9999
SYMIG	62	-1038.8637	2231.3182	0.0001	1.0000
HKY	59	-1061.0963	2266.7094	0.0000	1.0000
GTR	63	-1056.1656	2268.9928	0.0000	1.0000
K2P	56	-1069.6604	2274.9652	0.0000	1.0000
F81I	59	-1067.7826	2280.0821	0.0000	1.0000
SYM	60	-1066.5359	2280.5907	0.0000	1.0000
F81G	59	-1068.1426	2280.8021	0.0000	1.0000
F81IG	60	-1068.1423	2283.8035	0.0000	1.0000
JC69I	56	-1074.6582	2284.9609	0.0000	1.0000
JC69G	56	-1075.2714	2286.1872	0.0000	1.0000
JC69IG	57	-1074.7607	2288.1014	0.0000	1.0000
F81	58	-1095.3187	2332.1746	0.0000	1.0000
JC69	55	-1102.9992	2338.7290	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-1028.9450	2185.8900	0.3048	0.3048
GTRG	64	-1029.2898	2186.5797	0.2159	0.5207
HKYI	60	-1033.6026	2187.2052	0.1579	0.6786
HKYG	60	-1033.7270	2187.4540	0.1394	0.8181
GTRIG	65	-1028.8190	2187.6380	0.1272	0.9452
HKYIG	61	-1033.7268	2189.4536	0.0513	0.9965
K2PI	57	-1041.4088	2196.8176	0.0013	0.9978
K2PG	57	-1041.4340	2196.8681	0.0013	0.9991
K2PIG	58	-1041.5325	2199.0650	0.0004	0.9995
SYMI	61	-1039.1030	2200.2059	0.0002	0.9998
SYMG	61	-1039.6438	2201.2876	0.0001	0.9999
SYMIG	62	-1038.8637	2201.7273	0.0001	1.0000
GTR	63	-1056.1656	2238.3312	0.0000	1.0000
HKY	59	-1061.0963	2240.1925	0.0000	1.0000
K2P	56	-1069.6604	2251.3208	0.0000	1.0000
SYM	60	-1066.5359	2253.0719	0.0000	1.0000
F81I	59	-1067.7826	2253.5653	0.0000	1.0000
F81G	59	-1068.1426	2254.2852	0.0000	1.0000
F81IG	60	-1068.1423	2256.2847	0.0000	1.0000
JC69I	56	-1074.6582	2261.3165	0.0000	1.0000
JC69G	56	-1075.2714	2262.5428	0.0000	1.0000
JC69IG	57	-1074.7607	2263.5215	0.0000	1.0000
F81	58	-1095.3187	2306.6373	0.0000	1.0000
JC69	55	-1102.9992	2315.9984	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-1041.4088	2412.8453	0.3537	0.3537
K2PG	57	-1041.4340	2412.8958	0.3449	0.6986
HKYI	60	-1033.6026	2414.6028	0.1469	0.8455
HKYG	60	-1033.7270	2414.8516	0.1297	0.9752
K2PIG	58	-1041.5325	2418.8827	0.0173	0.9925
HKYIG	61	-1033.7268	2420.6412	0.0072	0.9997
GTRI	64	-1028.9450	2428.4475	0.0001	0.9998
GTRG	64	-1029.2898	2429.1371	0.0001	0.9999
SYMI	61	-1039.1030	2431.3935	0.0000	1.0000
SYMG	61	-1039.6438	2432.4752	0.0000	1.0000
GTRIG	65	-1028.8190	2433.9855	0.0000	1.0000
SYMIG	62	-1038.8637	2436.7048	0.0000	1.0000
K2P	56	-1069.6604	2463.5586	0.0000	1.0000
HKY	59	-1061.0963	2463.8002	0.0000	1.0000
JC69I	56	-1074.6582	2473.5543	0.0000	1.0000
JC69G	56	-1075.2714	2474.7805	0.0000	1.0000
GTR	63	-1056.1656	2477.0987	0.0000	1.0000
F81I	59	-1067.7826	2477.1729	0.0000	1.0000
F81G	59	-1068.1426	2477.8929	0.0000	1.0000
JC69IG	57	-1074.7607	2479.5492	0.0000	1.0000
SYM	60	-1066.5359	2480.4695	0.0000	1.0000
F81IG	60	-1068.1423	2483.6823	0.0000	1.0000
JC69	55	-1102.9992	2524.4462	0.0000	1.0000
F81	58	-1095.3187	2526.4550	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4076_hg18.phylip = 1 - 327;
 Partition Dummy = 1:261-aggressive-phylip/4076_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4076_hg18.phylip = 1 - 327;
 Partition Dummy = 1:261-aggressive-phylip/4076_hg18.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 (K2PI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4076_hg18.phylip = 1 - 327;
 Partition Dummy = 1:261-aggressive-phylip/4076_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

