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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	64	-1052.7743	2258.2371	0.7295	0.7295
GTRIG	65	-1052.6840	2260.9036	0.1923	0.9218
GTRG	64	-1055.0158	2262.7201	0.0775	0.9994
HKYI	60	-1066.2780	2274.0222	0.0003	0.9996
SYMI	61	-1065.3770	2275.0010	0.0002	0.9998
HKYIG	61	-1066.0304	2276.3079	0.0001	0.9999
SYMIG	62	-1065.3399	2277.7241	0.0000	0.9999
SYMG	61	-1066.8036	2277.8543	0.0000	1.0000
HKYG	60	-1068.4471	2278.3606	0.0000	1.0000
K2PI	57	-1077.3497	2287.9204	0.0000	1.0000
K2PIG	58	-1076.8064	2289.5662	0.0000	1.0000
K2PG	57	-1078.9459	2291.1127	0.0000	1.0000
GTR	63	-1079.3156	2308.4892	0.0000	1.0000
HKY	59	-1091.9056	2322.5129	0.0000	1.0000
SYM	60	-1090.5390	2322.5443	0.0000	1.0000
K2P	56	-1104.4404	2339.3852	0.0000	1.0000
F81I	59	-1101.2230	2341.1477	0.0000	1.0000
F81G	59	-1101.5823	2341.8664	0.0000	1.0000
F81IG	60	-1100.8945	2343.2552	0.0000	1.0000
JC69G	56	-1111.1372	2352.7787	0.0000	1.0000
JC69I	56	-1111.1563	2352.8170	0.0000	1.0000
JC69IG	57	-1110.8258	2354.8726	0.0000	1.0000
F81	58	-1127.3957	2390.7448	0.0000	1.0000
JC69	55	-1136.6997	2401.2029	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-1052.7743	2233.5487	0.6626	0.6626
GTRIG	65	-1052.6840	2235.3679	0.2668	0.9294
GTRG	64	-1055.0158	2238.0316	0.0704	0.9998
HKYI	60	-1066.2780	2252.5559	0.0000	0.9999
SYMI	61	-1065.3770	2252.7539	0.0000	0.9999
HKYIG	61	-1066.0304	2254.0608	0.0000	1.0000
SYMIG	62	-1065.3399	2254.6798	0.0000	1.0000
SYMG	61	-1066.8036	2255.6072	0.0000	1.0000
HKYG	60	-1068.4471	2256.8943	0.0000	1.0000
K2PI	57	-1077.3497	2268.6995	0.0000	1.0000
K2PIG	58	-1076.8064	2269.6129	0.0000	1.0000
K2PG	57	-1078.9459	2271.8918	0.0000	1.0000
GTR	63	-1079.3156	2284.6312	0.0000	1.0000
SYM	60	-1090.5390	2301.0780	0.0000	1.0000
HKY	59	-1091.9056	2301.8112	0.0000	1.0000
F81I	59	-1101.2230	2320.4459	0.0000	1.0000
K2P	56	-1104.4404	2320.8809	0.0000	1.0000
F81G	59	-1101.5823	2321.1646	0.0000	1.0000
F81IG	60	-1100.8945	2321.7889	0.0000	1.0000
JC69G	56	-1111.1372	2334.2744	0.0000	1.0000
JC69I	56	-1111.1563	2334.3127	0.0000	1.0000
JC69IG	57	-1110.8258	2335.6517	0.0000	1.0000
F81	58	-1127.3957	2370.7915	0.0000	1.0000
JC69	55	-1136.6997	2383.3995	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	64	-1052.7743	2489.3216	0.6704	0.6704
HKYI	60	-1066.2780	2492.3431	0.1480	0.8184
GTRG	64	-1055.0158	2493.8046	0.0713	0.8897
GTRIG	65	-1052.6840	2495.1373	0.0366	0.9263
K2PI	57	-1077.3497	2496.4972	0.0185	0.9448
SYMI	61	-1065.3770	2496.5375	0.0182	0.9630
HKYG	60	-1068.4471	2496.6814	0.0169	0.9799
HKYIG	61	-1066.0304	2497.8444	0.0095	0.9893
SYMG	61	-1066.8036	2499.3908	0.0044	0.9937
K2PG	57	-1078.9459	2499.6895	0.0038	0.9975
K2PIG	58	-1076.8064	2501.4071	0.0016	0.9991
SYMIG	62	-1065.3399	2502.4598	0.0009	1.0000
GTR	63	-1079.3156	2536.4077	0.0000	1.0000
HKY	59	-1091.9056	2537.6018	0.0000	1.0000
SYM	60	-1090.5390	2540.8651	0.0000	1.0000
K2P	56	-1104.4404	2544.6822	0.0000	1.0000
F81I	59	-1101.2230	2556.2366	0.0000	1.0000
F81G	59	-1101.5823	2556.9553	0.0000	1.0000
JC69G	56	-1111.1372	2558.0757	0.0000	1.0000
JC69I	56	-1111.1563	2558.1140	0.0000	1.0000
F81IG	60	-1100.8945	2561.5761	0.0000	1.0000
JC69IG	57	-1110.8258	2563.4494	0.0000	1.0000
F81	58	-1127.3957	2602.5857	0.0000	1.0000
JC69	55	-1136.6997	2603.2043	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/4983_hg18.phylip = 1 - 402;
 Partition Dummy = 1:261-aggressive-phylip/4983_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 (GTRI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4983_hg18.phylip = 1 - 402;
 Partition Dummy = 1:261-aggressive-phylip/4983_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 BIC model (GTRI)]

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


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

