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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1815.6379	3769.7606	0.7558	0.7558
HKYIG	61	-1815.6450	3772.4394	0.1980	0.9538
GTRG	64	-1813.3027	3775.8300	0.0363	0.9902
GTRIG	65	-1813.3100	3778.5637	0.0093	0.9995
K2PG	57	-1827.0797	3784.7309	0.0004	0.9999
K2PIG	58	-1827.0869	3787.3699	0.0001	1.0000
SYMG	61	-1826.5095	3794.1683	0.0000	1.0000
SYMIG	62	-1826.5165	3796.8604	0.0000	1.0000
HKY	59	-1851.7587	3839.3512	0.0000	1.0000
HKYI	60	-1851.3209	3841.1267	0.0000	1.0000
GTR	63	-1849.7313	3845.9818	0.0000	1.0000
GTRI	64	-1849.0356	3847.2958	0.0000	1.0000
K2P	56	-1864.0473	3856.0546	0.0000	1.0000
K2PI	57	-1863.6297	3857.8309	0.0000	1.0000
SYM	60	-1863.2935	3865.0718	0.0000	1.0000
SYMI	61	-1862.8557	3866.8608	0.0000	1.0000
F81G	59	-1869.6097	3875.0532	0.0000	1.0000
F81IG	60	-1869.6169	3877.7186	0.0000	1.0000
JC69G	56	-1878.1974	3884.3549	0.0000	1.0000
JC69IG	57	-1878.2048	3886.9809	0.0000	1.0000
F81	58	-1904.3003	3941.7967	0.0000	1.0000
F81I	59	-1904.1010	3944.0357	0.0000	1.0000
JC69	55	-1913.2517	3951.8651	0.0000	1.0000
JC69I	56	-1913.0367	3954.0335	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1815.6379	3751.2758	0.6158	0.6158
HKYIG	61	-1815.6450	3753.2901	0.2249	0.8407
GTRG	64	-1813.3027	3754.6055	0.1165	0.9573
GTRIG	65	-1813.3100	3756.6200	0.0426	0.9998
K2PG	57	-1827.0797	3768.1595	0.0001	0.9999
K2PIG	58	-1827.0869	3770.1739	0.0000	1.0000
SYMG	61	-1826.5095	3775.0190	0.0000	1.0000
SYMIG	62	-1826.5165	3777.0330	0.0000	1.0000
HKY	59	-1851.7587	3821.5174	0.0000	1.0000
HKYI	60	-1851.3209	3822.6419	0.0000	1.0000
GTR	63	-1849.7313	3825.4627	0.0000	1.0000
GTRI	64	-1849.0356	3826.0713	0.0000	1.0000
K2P	56	-1864.0473	3840.0946	0.0000	1.0000
K2PI	57	-1863.6297	3841.2595	0.0000	1.0000
SYM	60	-1863.2935	3846.5870	0.0000	1.0000
SYMI	61	-1862.8557	3847.7115	0.0000	1.0000
F81G	59	-1869.6097	3857.2194	0.0000	1.0000
F81IG	60	-1869.6169	3859.2338	0.0000	1.0000
JC69G	56	-1878.1974	3868.3949	0.0000	1.0000
JC69IG	57	-1878.2048	3870.4095	0.0000	1.0000
F81	58	-1904.3003	3924.6007	0.0000	1.0000
F81I	59	-1904.1010	3926.2020	0.0000	1.0000
JC69	55	-1913.2517	3936.5035	0.0000	1.0000
JC69I	56	-1913.0367	3938.0735	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1815.6379	3998.7568	0.8649	0.8649
K2PG	57	-1827.0797	4003.2664	0.0907	0.9556
HKYIG	61	-1815.6450	4004.8957	0.0402	0.9957
K2PIG	58	-1827.0869	4009.4055	0.0042	1.0000
GTRG	64	-1813.3027	4018.5852	0.0000	1.0000
GTRIG	65	-1813.3100	4024.7244	0.0000	1.0000
SYMG	61	-1826.5095	4026.6247	0.0000	1.0000
SYMIG	62	-1826.5165	4032.7634	0.0000	1.0000
HKY	59	-1851.7587	4064.8737	0.0000	1.0000
HKYI	60	-1851.3209	4070.1229	0.0000	1.0000
K2P	56	-1864.0473	4071.0769	0.0000	1.0000
K2PI	57	-1863.6297	4076.3664	0.0000	1.0000
GTR	63	-1849.7313	4085.3177	0.0000	1.0000
GTRI	64	-1849.0356	4090.0510	0.0000	1.0000
SYM	60	-1863.2935	4094.0680	0.0000	1.0000
SYMI	61	-1862.8557	4099.3171	0.0000	1.0000
JC69G	56	-1878.1974	4099.3771	0.0000	1.0000
F81G	59	-1869.6097	4100.5757	0.0000	1.0000
JC69IG	57	-1878.2048	4105.5165	0.0000	1.0000
F81IG	60	-1869.6169	4106.7148	0.0000	1.0000
JC69	55	-1913.2517	4163.3611	0.0000	1.0000
F81	58	-1904.3003	4163.8323	0.0000	1.0000
JC69I	56	-1913.0367	4169.0558	0.0000	1.0000
F81I	59	-1904.1010	4169.5583	0.0000	1.0000



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

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

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


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

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


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

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

