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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1574.5752	3313.1677	0.9753	0.9753
GTRIG	65	-1576.6705	3320.6452	0.0232	0.9985
GTRI	64	-1581.6990	3327.4153	0.0008	0.9993
HKYG	60	-1589.1600	3329.4689	0.0003	0.9996
SYMG	61	-1588.1213	3330.5674	0.0002	0.9998
HKYIG	61	-1588.1636	3330.6521	0.0002	0.9999
SYMIG	62	-1587.3363	3332.2004	0.0001	1.0000
K2PIG	58	-1599.4956	3343.8687	0.0000	1.0000
SYMI	61	-1596.9746	3348.2741	0.0000	1.0000
HKYI	60	-1602.0022	3355.1533	0.0000	1.0000
K2PG	57	-1607.0094	3355.8003	0.0000	1.0000
K2PI	57	-1617.5306	3376.8428	0.0000	1.0000
F81G	59	-1622.8994	3393.7988	0.0000	1.0000
F81IG	60	-1622.1378	3395.4246	0.0000	1.0000
JC69G	56	-1634.8876	3408.4864	0.0000	1.0000
JC69IG	57	-1635.0664	3411.9143	0.0000	1.0000
F81I	59	-1634.7645	3417.5290	0.0000	1.0000
JC69I	56	-1642.6961	3424.1036	0.0000	1.0000
GTR	63	-1672.0697	3504.8980	0.0000	1.0000
SYM	60	-1688.7229	3528.5947	0.0000	1.0000
HKY	59	-1699.3909	3546.7817	0.0000	1.0000
K2P	56	-1715.4664	3569.6441	0.0000	1.0000
F81	58	-1721.5631	3588.0039	0.0000	1.0000
JC69	55	-1734.1921	3604.0508	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1574.5752	3277.1503	0.9559	0.9559
GTRIG	65	-1576.6705	3283.3409	0.0433	0.9991
GTRI	64	-1581.6990	3291.3980	0.0008	0.9999
SYMG	61	-1588.1213	3298.2426	0.0000	0.9999
HKYG	60	-1589.1600	3298.3200	0.0000	1.0000
HKYIG	61	-1588.1636	3298.3273	0.0000	1.0000
SYMIG	62	-1587.3363	3298.6725	0.0000	1.0000
K2PIG	58	-1599.4956	3314.9911	0.0000	1.0000
SYMI	61	-1596.9746	3315.9493	0.0000	1.0000
HKYI	60	-1602.0022	3324.0044	0.0000	1.0000
K2PG	57	-1607.0094	3328.0188	0.0000	1.0000
K2PI	57	-1617.5306	3349.0613	0.0000	1.0000
F81G	59	-1622.8994	3363.7988	0.0000	1.0000
F81IG	60	-1622.1378	3364.2757	0.0000	1.0000
JC69G	56	-1634.8876	3381.7751	0.0000	1.0000
JC69IG	57	-1635.0664	3384.1328	0.0000	1.0000
F81I	59	-1634.7645	3387.5290	0.0000	1.0000
JC69I	56	-1642.6961	3397.3923	0.0000	1.0000
GTR	63	-1672.0697	3470.1394	0.0000	1.0000
SYM	60	-1688.7229	3497.4457	0.0000	1.0000
HKY	59	-1699.3909	3516.7817	0.0000	1.0000
K2P	56	-1715.4664	3542.9328	0.0000	1.0000
F81	58	-1721.5631	3559.1263	0.0000	1.0000
JC69	55	-1734.1921	3578.3842	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1574.5752	3513.3333	0.9409	0.9409
HKYG	60	-1589.1600	3519.7415	0.0382	0.9791
GTRIG	65	-1576.6705	3523.2143	0.0067	0.9858
SYMG	61	-1588.1213	3523.3546	0.0063	0.9921
HKYIG	61	-1588.1636	3523.4392	0.0060	0.9981
SYMIG	62	-1587.3363	3527.4748	0.0008	0.9989
GTRI	64	-1581.6990	3527.5810	0.0008	0.9996
K2PIG	58	-1599.4956	3529.0319	0.0004	1.0000
K2PG	57	-1607.0094	3538.3693	0.0000	1.0000
SYMI	61	-1596.9746	3541.0612	0.0000	1.0000
HKYI	60	-1602.0022	3545.4259	0.0000	1.0000
K2PI	57	-1617.5306	3559.4118	0.0000	1.0000
F81G	59	-1622.8994	3581.5300	0.0000	1.0000
F81IG	60	-1622.1378	3585.6972	0.0000	1.0000
JC69G	56	-1634.8876	3588.4353	0.0000	1.0000
JC69IG	57	-1635.0664	3594.4832	0.0000	1.0000
JC69I	56	-1642.6961	3604.0524	0.0000	1.0000
F81I	59	-1634.7645	3605.2602	0.0000	1.0000
GTR	63	-1672.0697	3702.6320	0.0000	1.0000
SYM	60	-1688.7229	3718.8673	0.0000	1.0000
HKY	59	-1699.3909	3734.5129	0.0000	1.0000
K2P	56	-1715.4664	3749.5929	0.0000	1.0000
F81	58	-1721.5631	3773.1671	0.0000	1.0000
JC69	55	-1734.1921	3781.3539	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/3573_hg18.phylip = 1 - 296;
 Partition Dummy = 1:261-aggressive-phylip/3573_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/3573_hg18.phylip = 1 - 296;
 Partition Dummy = 1:261-aggressive-phylip/3573_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 (GTRG)]

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


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

