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

Input data from file "261-aggressive-phylip/11562_hg18.phylip" (29 taxa, 407 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	-1733.1202	3607.3965	0.9152	0.9152
GTRG	64	-1730.6038	3613.5350	0.0425	0.9577
HKYIG	61	-1734.9068	3613.7382	0.0384	0.9961
GTRIG	65	-1732.2806	3619.7226	0.0019	0.9980
HKYI	60	-1739.2980	3619.7521	0.0019	0.9999
GTRI	64	-1737.0207	3626.3689	0.0001	1.0000
F81G	59	-1749.3251	3637.0536	0.0000	1.0000
F81IG	60	-1750.4392	3642.0344	0.0000	1.0000
F81I	59	-1754.7514	3647.9063	0.0000	1.0000
SYMIG	62	-1759.7333	3666.1760	0.0000	1.0000
K2PG	57	-1767.3459	3667.6373	0.0000	1.0000
SYMG	61	-1762.1009	3668.1264	0.0000	1.0000
K2PIG	58	-1769.5243	3674.7153	0.0000	1.0000
K2PI	57	-1777.0637	3687.0729	0.0000	1.0000
SYMI	61	-1772.3947	3688.7140	0.0000	1.0000
JC69G	56	-1782.2690	3694.7780	0.0000	1.0000
JC69IG	57	-1784.9291	3702.8037	0.0000	1.0000
JC69I	56	-1791.4726	3713.1851	0.0000	1.0000
HKY	59	-1827.1689	3792.7412	0.0000	1.0000
GTR	63	-1822.8343	3795.1788	0.0000	1.0000
F81	58	-1842.3069	3820.2804	0.0000	1.0000
SYM	60	-1844.4353	3830.0268	0.0000	1.0000
K2P	56	-1854.9235	3840.0871	0.0000	1.0000
JC69	55	-1868.7144	3864.9786	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1733.1202	3586.2404	0.7654	0.7654
GTRG	64	-1730.6038	3589.2075	0.1736	0.9390
HKYIG	61	-1734.9068	3591.8136	0.0472	0.9862
GTRIG	65	-1732.2806	3594.5613	0.0119	0.9981
HKYI	60	-1739.2980	3598.5960	0.0016	0.9997
GTRI	64	-1737.0207	3602.0415	0.0003	1.0000
F81G	59	-1749.3251	3616.6502	0.0000	1.0000
F81IG	60	-1750.4392	3620.8784	0.0000	1.0000
F81I	59	-1754.7514	3627.5029	0.0000	1.0000
SYMIG	62	-1759.7333	3643.4667	0.0000	1.0000
SYMG	61	-1762.1009	3646.2018	0.0000	1.0000
K2PG	57	-1767.3459	3648.6918	0.0000	1.0000
K2PIG	58	-1769.5243	3655.0486	0.0000	1.0000
SYMI	61	-1772.3947	3666.7893	0.0000	1.0000
K2PI	57	-1777.0637	3668.1273	0.0000	1.0000
JC69G	56	-1782.2690	3676.5380	0.0000	1.0000
JC69IG	57	-1784.9291	3683.8581	0.0000	1.0000
JC69I	56	-1791.4726	3694.9451	0.0000	1.0000
GTR	63	-1822.8343	3771.6686	0.0000	1.0000
HKY	59	-1827.1689	3772.3377	0.0000	1.0000
F81	58	-1842.3069	3800.6137	0.0000	1.0000
SYM	60	-1844.4353	3808.8707	0.0000	1.0000
K2P	56	-1854.9235	3821.8471	0.0000	1.0000
JC69	55	-1868.7144	3847.4288	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1733.1202	3826.7692	0.9897	0.9897
HKYIG	61	-1734.9068	3836.3512	0.0082	0.9979
HKYI	60	-1739.2980	3839.1248	0.0021	0.9999
GTRG	64	-1730.6038	3845.7716	0.0001	1.0000
F81G	59	-1749.3251	3853.1701	0.0000	1.0000
GTRIG	65	-1732.2806	3855.1341	0.0000	1.0000
GTRI	64	-1737.0207	3858.6055	0.0000	1.0000
F81IG	60	-1750.4392	3861.4072	0.0000	1.0000
F81I	59	-1754.7514	3864.0229	0.0000	1.0000
K2PG	57	-1767.3459	3877.1941	0.0000	1.0000
K2PIG	58	-1769.5243	3887.5598	0.0000	1.0000
SYMG	61	-1762.1009	3890.7394	0.0000	1.0000
SYMIG	62	-1759.7333	3892.0131	0.0000	1.0000
K2PI	57	-1777.0637	3896.6297	0.0000	1.0000
JC69G	56	-1782.2690	3901.0315	0.0000	1.0000
SYMI	61	-1772.3947	3911.3269	0.0000	1.0000
JC69IG	57	-1784.9291	3912.3605	0.0000	1.0000
JC69I	56	-1791.4726	3919.4387	0.0000	1.0000
HKY	59	-1827.1689	4008.8577	0.0000	1.0000
GTR	63	-1822.8343	4024.2238	0.0000	1.0000
F81	58	-1842.3069	4033.1249	0.0000	1.0000
K2P	56	-1854.9235	4046.3406	0.0000	1.0000
SYM	60	-1844.4353	4049.3995	0.0000	1.0000
JC69	55	-1868.7144	4067.9135	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/11562_hg18.phylip = 1 - 407;
 Partition Dummy = 1:261-aggressive-phylip/11562_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/11562_hg18.phylip = 1 - 407;
 Partition Dummy = 1:261-aggressive-phylip/11562_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/11562_hg18.phylip = 1 - 407;
 Partition Dummy = 1:261-aggressive-phylip/11562_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

