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

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

Minimum AIC  model: GTRIG
Minimum AICc model: HKYIG
Minimum BIC  model: HKYG


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
HKYIG	61	-1087.5180	2334.2972	0.4984	0.4984
HKYG	60	-1089.6401	2335.1626	0.3233	0.8217
GTRG	64	-1084.0107	2337.6214	0.0946	0.9162
GTRIG	65	-1082.3774	2337.8703	0.0835	0.9997
HKYI	60	-1096.8242	2349.5308	0.0002	1.0000
GTRI	64	-1092.2165	2354.0330	0.0000	1.0000
K2PIG	58	-1104.5857	2358.3948	0.0000	1.0000
SYMG	61	-1099.6441	2358.5492	0.0000	1.0000
SYMIG	62	-1097.9462	2358.5657	0.0000	1.0000
K2PG	57	-1106.8107	2359.5634	0.0000	1.0000
F81IG	60	-1105.8117	2367.5058	0.0000	1.0000
F81G	59	-1107.7784	2368.0933	0.0000	1.0000
K2PI	57	-1113.0359	2372.0139	0.0000	1.0000
SYMI	61	-1107.1241	2373.5093	0.0000	1.0000
F81I	59	-1113.8981	2380.3328	0.0000	1.0000
JC69IG	57	-1120.0862	2386.1144	0.0000	1.0000
JC69G	56	-1122.1182	2386.9286	0.0000	1.0000
JC69I	56	-1128.0188	2398.7299	0.0000	1.0000
GTR	63	-1125.0112	2416.1417	0.0000	1.0000
HKY	59	-1132.0200	2416.5765	0.0000	1.0000
SYM	60	-1143.6768	2443.2360	0.0000	1.0000
K2P	56	-1151.5723	2445.8369	0.0000	1.0000
F81	58	-1149.7425	2448.7084	0.0000	1.0000
JC69	55	-1166.0738	2471.6213	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1082.3774	2294.7547	0.5116	0.5116
GTRG	64	-1084.0107	2296.0214	0.2716	0.7831
HKYIG	61	-1087.5180	2297.0361	0.1635	0.9466
HKYG	60	-1089.6401	2299.2803	0.0532	0.9999
GTRI	64	-1092.2165	2312.4330	0.0001	1.0000
HKYI	60	-1096.8242	2313.6484	0.0000	1.0000
SYMIG	62	-1097.9462	2319.8925	0.0000	1.0000
SYMG	61	-1099.6441	2321.2881	0.0000	1.0000
K2PIG	58	-1104.5857	2325.1715	0.0000	1.0000
K2PG	57	-1106.8107	2327.6214	0.0000	1.0000
F81IG	60	-1105.8117	2331.6234	0.0000	1.0000
F81G	59	-1107.7784	2333.5567	0.0000	1.0000
SYMI	61	-1107.1241	2336.2482	0.0000	1.0000
K2PI	57	-1113.0359	2340.0718	0.0000	1.0000
F81I	59	-1113.8981	2345.7962	0.0000	1.0000
JC69IG	57	-1120.0862	2354.1724	0.0000	1.0000
JC69G	56	-1122.1182	2356.2363	0.0000	1.0000
JC69I	56	-1128.0188	2368.0376	0.0000	1.0000
GTR	63	-1125.0112	2376.0223	0.0000	1.0000
HKY	59	-1132.0200	2382.0399	0.0000	1.0000
SYM	60	-1143.6768	2407.3537	0.0000	1.0000
K2P	56	-1151.5723	2415.1446	0.0000	1.0000
F81	58	-1149.7425	2415.4851	0.0000	1.0000
JC69	55	-1166.0738	2442.1476	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1089.6401	2514.0640	0.6583	0.6583
HKYIG	61	-1087.5180	2515.3996	0.3376	0.9959
GTRG	64	-1084.0107	2525.1241	0.0026	0.9985
GTRIG	65	-1082.3774	2527.4372	0.0008	0.9993
HKYI	60	-1096.8242	2528.4322	0.0005	0.9998
K2PG	57	-1106.8107	2531.6660	0.0001	0.9999
K2PIG	58	-1104.5857	2532.7958	0.0001	1.0000
SYMG	61	-1099.6441	2539.6517	0.0000	1.0000
GTRI	64	-1092.2165	2541.5357	0.0000	1.0000
SYMIG	62	-1097.9462	2541.8357	0.0000	1.0000
K2PI	57	-1113.0359	2544.1164	0.0000	1.0000
F81G	59	-1107.7784	2544.7608	0.0000	1.0000
F81IG	60	-1105.8117	2546.4072	0.0000	1.0000
SYMI	61	-1107.1241	2554.6117	0.0000	1.0000
JC69G	56	-1122.1182	2556.7012	0.0000	1.0000
F81I	59	-1113.8981	2557.0003	0.0000	1.0000
JC69IG	57	-1120.0862	2558.2170	0.0000	1.0000
JC69I	56	-1128.0188	2568.5025	0.0000	1.0000
HKY	59	-1132.0200	2593.2440	0.0000	1.0000
GTR	63	-1125.0112	2601.5453	0.0000	1.0000
K2P	56	-1151.5723	2615.6095	0.0000	1.0000
SYM	60	-1143.6768	2622.1374	0.0000	1.0000
F81	58	-1149.7425	2623.1094	0.0000	1.0000
JC69	55	-1166.0738	2639.0327	0.0000	1.0000



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

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

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


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

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


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12408_hg18.phylip = 1 - 265;
 Partition Dummy = 1:261-aggressive-phylip/12408_hg18.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;


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

