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

Input data from file "261-aggressive-phylip/13336_hg18.phylip" (29 taxa, 1147 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	-2059.1299	4245.0002	0.4981	0.4981
HKYIG	61	-2059.0019	4246.9753	0.1855	0.6837
GTRI	64	-2055.7880	4247.2654	0.1605	0.8442
HKYG	60	-2061.1762	4249.0928	0.0644	0.9085
GTRIG	65	-2055.6453	4249.2277	0.0602	0.9687
GTRG	64	-2057.5132	4250.7158	0.0286	0.9973
F81I	59	-2065.8386	4256.1905	0.0019	0.9991
F81IG	60	-2065.7364	4258.2131	0.0007	0.9998
F81G	59	-2068.1401	4260.7936	0.0002	1.0000
HKY	59	-2075.2483	4275.0100	0.0000	1.0000
GTR	63	-2071.4546	4276.3551	0.0000	1.0000
F81	58	-2082.4030	4287.0965	0.0000	1.0000
K2PI	57	-2101.0842	4322.2401	0.0000	1.0000
K2PIG	58	-2101.0019	4324.2942	0.0000	1.0000
K2PG	57	-2103.3994	4326.8705	0.0000	1.0000
SYMI	61	-2099.4029	4327.7773	0.0000	1.0000
SYMIG	62	-2099.3414	4329.8894	0.0000	1.0000
JC69I	56	-2106.9789	4331.8147	0.0000	1.0000
SYMG	61	-2101.4663	4331.9041	0.0000	1.0000
JC69IG	57	-2106.8855	4333.8425	0.0000	1.0000
JC69G	56	-2109.4177	4336.6922	0.0000	1.0000
K2P	56	-2117.6494	4353.1557	0.0000	1.0000
SYM	60	-2115.4657	4357.6717	0.0000	1.0000
JC69	55	-2123.7605	4363.1671	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-2059.1299	4238.2599	0.4198	0.4198
GTRI	64	-2055.7880	4239.5759	0.2174	0.6372
HKYIG	61	-2059.0019	4240.0039	0.1755	0.8127
GTRIG	65	-2055.6453	4241.2906	0.0922	0.9049
HKYG	60	-2061.1762	4242.3525	0.0542	0.9592
GTRG	64	-2057.5132	4243.0263	0.0387	0.9979
F81I	59	-2065.8386	4249.6771	0.0014	0.9993
F81IG	60	-2065.7364	4251.4727	0.0006	0.9999
F81G	59	-2068.1401	4254.2803	0.0001	1.0000
HKY	59	-2075.2483	4268.4966	0.0000	1.0000
GTR	63	-2071.4546	4268.9091	0.0000	1.0000
F81	58	-2082.4030	4280.8060	0.0000	1.0000
K2PI	57	-2101.0842	4316.1684	0.0000	1.0000
K2PIG	58	-2101.0019	4318.0038	0.0000	1.0000
K2PG	57	-2103.3994	4320.7989	0.0000	1.0000
SYMI	61	-2099.4029	4320.8059	0.0000	1.0000
SYMIG	62	-2099.3414	4322.6827	0.0000	1.0000
SYMG	61	-2101.4663	4324.9327	0.0000	1.0000
JC69I	56	-2106.9789	4325.9578	0.0000	1.0000
JC69IG	57	-2106.8855	4327.7709	0.0000	1.0000
JC69G	56	-2109.4177	4330.8353	0.0000	1.0000
K2P	56	-2117.6494	4347.2988	0.0000	1.0000
SYM	60	-2115.4657	4350.9313	0.0000	1.0000
JC69	55	-2123.7605	4357.5209	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-2059.1299	4540.9542	0.8267	0.8267
HKYG	60	-2061.1762	4545.0468	0.1068	0.9335
F81I	59	-2065.8386	4547.3265	0.0342	0.9677
HKYIG	61	-2059.0019	4547.7431	0.0277	0.9954
F81G	59	-2068.1401	4551.9297	0.0034	0.9989
F81IG	60	-2065.7364	4554.1670	0.0011	1.0000
GTRI	64	-2055.7880	4562.4498	0.0000	1.0000
GTRG	64	-2057.5132	4565.9002	0.0000	1.0000
HKY	59	-2075.2483	4566.1460	0.0000	1.0000
GTRIG	65	-2055.6453	4569.2095	0.0000	1.0000
F81	58	-2082.4030	4573.4105	0.0000	1.0000
GTR	63	-2071.4546	4586.7382	0.0000	1.0000
K2PI	57	-2101.0842	4603.7280	0.0000	1.0000
K2PG	57	-2103.3994	4608.3585	0.0000	1.0000
JC69I	56	-2106.9789	4608.4725	0.0000	1.0000
K2PIG	58	-2101.0019	4610.6083	0.0000	1.0000
JC69G	56	-2109.4177	4613.3500	0.0000	1.0000
JC69IG	57	-2106.8855	4615.3305	0.0000	1.0000
SYMI	61	-2099.4029	4628.5451	0.0000	1.0000
K2P	56	-2117.6494	4629.8135	0.0000	1.0000
SYMG	61	-2101.4663	4632.6719	0.0000	1.0000
JC69	55	-2123.7605	4634.9907	0.0000	1.0000
SYMIG	62	-2099.3414	4635.4669	0.0000	1.0000
SYM	60	-2115.4657	4653.6256	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 261-aggressive-phylip/13336_hg18.phylip = 1 - 1147;
 Partition Dummy = 1:261-aggressive-phylip/13336_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 AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13336_hg18.phylip = 1 - 1147;
 Partition Dummy = 1:261-aggressive-phylip/13336_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 (HKYI)]

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


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

