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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1048.7850	2254.5596	0.4363	0.4363
GTRI	64	-1048.8172	2254.6238	0.4225	0.8588
GTRIG	65	-1048.6499	2257.2999	0.1108	0.9696
HKYG	60	-1058.0206	2261.1958	0.0158	0.9854
HKYI	60	-1058.4786	2262.1119	0.0100	0.9954
HKYIG	61	-1057.8663	2263.8154	0.0043	0.9997
K2PG	57	-1067.4385	2271.3668	0.0001	0.9998
SYMG	61	-1061.8926	2271.8680	0.0001	0.9999
K2PI	57	-1068.1394	2272.7686	0.0000	0.9999
SYMI	61	-1062.3667	2272.8162	0.0000	0.9999
K2PIG	58	-1067.1386	2273.6356	0.0000	1.0000
SYMIG	62	-1061.7379	2274.5070	0.0000	1.0000
GTR	63	-1064.1694	2282.3389	0.0000	1.0000
HKY	59	-1073.6681	2289.5827	0.0000	1.0000
SYM	60	-1077.7245	2300.6036	0.0000	1.0000
K2P	56	-1083.9814	2301.6035	0.0000	1.0000
F81G	59	-1085.5138	2313.2743	0.0000	1.0000
F81I	59	-1086.0712	2314.3889	0.0000	1.0000
F81IG	60	-1085.3131	2315.7809	0.0000	1.0000
JC69G	56	-1093.3565	2320.3536	0.0000	1.0000
JC69I	56	-1094.0221	2321.6850	0.0000	1.0000
JC69IG	57	-1093.0856	2322.6609	0.0000	1.0000
F81	58	-1101.2461	2341.8506	0.0000	1.0000
JC69	55	-1109.5735	2349.9578	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1048.7850	2225.5701	0.4166	0.4166
GTRI	64	-1048.8172	2225.6343	0.4034	0.8200
GTRIG	65	-1048.6499	2227.2999	0.1754	0.9954
HKYG	60	-1058.0206	2236.0412	0.0022	0.9976
HKYI	60	-1058.4786	2236.9573	0.0014	0.9990
HKYIG	61	-1057.8663	2237.7327	0.0010	1.0000
SYMG	61	-1061.8926	2245.7853	0.0000	1.0000
SYMI	61	-1062.3667	2246.7334	0.0000	1.0000
SYMIG	62	-1061.7379	2247.4759	0.0000	1.0000
K2PG	57	-1067.4385	2248.8770	0.0000	1.0000
K2PIG	58	-1067.1386	2250.2772	0.0000	1.0000
K2PI	57	-1068.1394	2250.2788	0.0000	1.0000
GTR	63	-1064.1694	2254.3389	0.0000	1.0000
HKY	59	-1073.6681	2265.3362	0.0000	1.0000
SYM	60	-1077.7245	2275.4489	0.0000	1.0000
K2P	56	-1083.9814	2279.9628	0.0000	1.0000
F81G	59	-1085.5138	2289.0277	0.0000	1.0000
F81I	59	-1086.0712	2290.1423	0.0000	1.0000
F81IG	60	-1085.3131	2290.6262	0.0000	1.0000
JC69G	56	-1093.3565	2298.7129	0.0000	1.0000
JC69I	56	-1094.0221	2300.0443	0.0000	1.0000
JC69IG	57	-1093.0856	2300.1711	0.0000	1.0000
F81	58	-1101.2461	2318.4922	0.0000	1.0000
JC69	55	-1109.5735	2329.1470	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1058.0206	2467.8591	0.3695	0.3695
HKYI	60	-1058.4786	2468.7751	0.2337	0.6033
K2PG	57	-1067.4385	2469.1040	0.1983	0.8016
K2PI	57	-1068.1394	2470.5058	0.0984	0.9000
GTRG	64	-1048.7850	2472.8425	0.0306	0.9306
GTRI	64	-1048.8172	2472.9067	0.0296	0.9602
HKYIG	61	-1057.8663	2473.4142	0.0230	0.9832
K2PIG	58	-1067.1386	2474.3678	0.0143	0.9974
GTRIG	65	-1048.6499	2478.4359	0.0019	0.9993
SYMG	61	-1061.8926	2481.4668	0.0004	0.9997
SYMI	61	-1062.3667	2482.4149	0.0003	1.0000
SYMIG	62	-1061.7379	2487.0210	0.0000	1.0000
HKY	59	-1073.6681	2493.2904	0.0000	1.0000
K2P	56	-1083.9814	2496.3262	0.0000	1.0000
GTR	63	-1064.1694	2497.7476	0.0000	1.0000
SYM	60	-1077.7245	2507.2668	0.0000	1.0000
JC69G	56	-1093.3565	2515.0763	0.0000	1.0000
JC69I	56	-1094.0221	2516.4076	0.0000	1.0000
F81G	59	-1085.5138	2516.9819	0.0000	1.0000
F81I	59	-1086.0712	2518.0966	0.0000	1.0000
JC69IG	57	-1093.0856	2520.3981	0.0000	1.0000
F81IG	60	-1085.3131	2522.4441	0.0000	1.0000
JC69	55	-1109.5735	2541.6467	0.0000	1.0000
F81	58	-1101.2461	2542.5828	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/2377_hg18.phylip = 1 - 352;
 Partition Dummy = 1:261-aggressive-phylip/2377_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/2377_hg18.phylip = 1 - 352;
 Partition Dummy = 1:261-aggressive-phylip/2377_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


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

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

