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

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

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


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	-1631.6243	3422.5268	0.7454	0.7454
GTRIG	65	-1631.1942	3424.7657	0.2434	0.9888
HKYG	60	-1642.8003	3432.7118	0.0046	0.9934
HKYIG	61	-1641.3253	3432.7696	0.0044	0.9978
SYMG	61	-1642.8494	3435.8177	0.0010	0.9988
SYMIG	62	-1641.4718	3436.0929	0.0008	0.9996
K2PG	57	-1649.8401	3437.8999	0.0003	1.0000
K2PIG	58	-1651.1797	3443.5212	0.0000	1.0000
GTRI	64	-1642.6206	3444.5194	0.0000	1.0000
HKYI	60	-1650.6211	3448.3533	0.0000	1.0000
SYMI	61	-1651.9796	3454.0781	0.0000	1.0000
K2PI	57	-1658.9184	3456.0567	0.0000	1.0000
JC69G	56	-1711.5011	3558.3015	0.0000	1.0000
JC69IG	57	-1710.3869	3558.9936	0.0000	1.0000
GTR	63	-1701.6440	3559.4902	0.0000	1.0000
F81IG	60	-1706.4827	3560.0766	0.0000	1.0000
F81G	59	-1708.6850	3561.4955	0.0000	1.0000
HKY	59	-1711.8593	3567.8441	0.0000	1.0000
SYM	60	-1712.0768	3571.2648	0.0000	1.0000
K2P	56	-1718.8808	3573.0608	0.0000	1.0000
F81I	59	-1714.7883	3573.7020	0.0000	1.0000
JC69I	56	-1719.4372	3574.1736	0.0000	1.0000
F81	58	-1772.5691	3686.3000	0.0000	1.0000
JC69	55	-1777.3883	3687.1765	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1631.6243	3391.2486	0.6376	0.6376
GTRIG	65	-1631.1942	3392.3883	0.3607	0.9983
HKYIG	61	-1641.3253	3404.6506	0.0008	0.9991
HKYG	60	-1642.8003	3405.6007	0.0005	0.9996
SYMIG	62	-1641.4718	3406.9436	0.0002	0.9998
SYMG	61	-1642.8494	3407.6987	0.0002	1.0000
GTRI	64	-1642.6206	3413.2412	0.0000	1.0000
K2PG	57	-1649.8401	3413.6801	0.0000	1.0000
K2PIG	58	-1651.1797	3418.3595	0.0000	1.0000
HKYI	60	-1650.6211	3421.2422	0.0000	1.0000
SYMI	61	-1651.9796	3425.9592	0.0000	1.0000
K2PI	57	-1658.9184	3431.8369	0.0000	1.0000
GTR	63	-1701.6440	3529.2880	0.0000	1.0000
F81IG	60	-1706.4827	3532.9655	0.0000	1.0000
JC69IG	57	-1710.3869	3534.7738	0.0000	1.0000
JC69G	56	-1711.5011	3535.0022	0.0000	1.0000
F81G	59	-1708.6850	3535.3701	0.0000	1.0000
HKY	59	-1711.8593	3541.7186	0.0000	1.0000
SYM	60	-1712.0768	3544.1537	0.0000	1.0000
F81I	59	-1714.7883	3547.5766	0.0000	1.0000
K2P	56	-1718.8808	3549.7615	0.0000	1.0000
JC69I	56	-1719.4372	3550.8743	0.0000	1.0000
F81	58	-1772.5691	3661.1382	0.0000	1.0000
JC69	55	-1777.3883	3664.7765	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1649.8401	3630.4008	0.7165	0.7165
HKYG	60	-1642.8003	3633.7278	0.1358	0.8523
GTRG	64	-1631.6243	3634.5842	0.0885	0.9408
HKYIG	61	-1641.3253	3636.5798	0.0326	0.9734
K2PIG	58	-1651.1797	3638.8823	0.0103	0.9837
GTRIG	65	-1631.1942	3639.5260	0.0075	0.9912
SYMG	61	-1642.8494	3639.6279	0.0071	0.9983
SYMIG	62	-1641.4718	3642.6750	0.0015	0.9999
K2PI	57	-1658.9184	3648.5576	0.0001	0.9999
HKYI	60	-1650.6211	3649.3693	0.0001	1.0000
GTRI	64	-1642.6206	3656.5768	0.0000	1.0000
SYMI	61	-1651.9796	3657.8884	0.0000	1.0000
JC69G	56	-1711.5011	3747.9208	0.0000	1.0000
JC69IG	57	-1710.3869	3751.4946	0.0000	1.0000
F81G	59	-1708.6850	3759.6950	0.0000	1.0000
F81IG	60	-1706.4827	3761.0926	0.0000	1.0000
K2P	56	-1718.8808	3762.6802	0.0000	1.0000
JC69I	56	-1719.4372	3763.7929	0.0000	1.0000
HKY	59	-1711.8593	3766.0436	0.0000	1.0000
GTR	63	-1701.6440	3768.8215	0.0000	1.0000
F81I	59	-1714.7883	3771.9016	0.0000	1.0000
SYM	60	-1712.0768	3772.2808	0.0000	1.0000
JC69	55	-1777.3883	3873.8931	0.0000	1.0000
F81	58	-1772.5691	3881.6611	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/6099_hg18.phylip = 1 - 331;
 Partition Dummy = 1:261-aggressive-phylip/6099_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/6099_hg18.phylip = 1 - 331;
 Partition Dummy = 1:261-aggressive-phylip/6099_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6099_hg18.phylip = 1 - 331;
 Partition Dummy = 1:261-aggressive-phylip/6099_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

