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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: K2PIG


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
GTRIG	65	-1702.6419	3555.5674	0.9278	0.9278
GTRG	64	-1706.9840	3561.5907	0.0457	0.9735
HKYIG	61	-1711.7330	3563.1802	0.0206	0.9941
SYMIG	62	-1712.6964	3567.7308	0.0021	0.9962
K2PIG	58	-1718.0460	3568.0083	0.0018	0.9981
HKYG	60	-1715.7286	3568.5601	0.0014	0.9995
GTRI	64	-1712.1096	3571.8419	0.0003	0.9997
SYMG	61	-1716.7519	3573.2181	0.0001	0.9999
K2PG	57	-1722.0932	3573.5274	0.0001	1.0000
HKYI	60	-1721.0373	3579.1773	0.0000	1.0000
SYMI	61	-1721.6853	3583.0850	0.0000	1.0000
K2PI	57	-1727.2560	3583.8532	0.0000	1.0000
F81IG	60	-1754.2333	3645.5694	0.0000	1.0000
JC69IG	57	-1760.2706	3649.8822	0.0000	1.0000
F81G	59	-1758.0318	3650.5671	0.0000	1.0000
JC69G	56	-1764.1463	3655.0704	0.0000	1.0000
F81I	59	-1762.5828	3659.6691	0.0000	1.0000
JC69I	56	-1768.6840	3664.1458	0.0000	1.0000
GTR	63	-1764.4523	3673.8787	0.0000	1.0000
SYM	60	-1773.7752	3684.6532	0.0000	1.0000
HKY	59	-1778.3576	3691.2187	0.0000	1.0000
K2P	56	-1783.9838	3694.7454	0.0000	1.0000
F81	58	-1815.7804	3763.4770	0.0000	1.0000
JC69	55	-1821.7810	3767.7884	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1702.6419	3535.2837	0.9586	0.9586
GTRG	64	-1706.9840	3541.9680	0.0339	0.9925
HKYIG	61	-1711.7330	3545.4659	0.0059	0.9984
SYMIG	62	-1712.6964	3549.3928	0.0008	0.9992
HKYG	60	-1715.7286	3551.4573	0.0003	0.9995
K2PIG	58	-1718.0460	3552.0920	0.0002	0.9997
GTRI	64	-1712.1096	3552.2193	0.0002	0.9999
SYMG	61	-1716.7519	3555.5039	0.0000	1.0000
K2PG	57	-1722.0932	3558.1863	0.0000	1.0000
HKYI	60	-1721.0373	3562.0745	0.0000	1.0000
SYMI	61	-1721.6853	3565.3707	0.0000	1.0000
K2PI	57	-1727.2560	3568.5121	0.0000	1.0000
F81IG	60	-1754.2333	3628.4666	0.0000	1.0000
F81G	59	-1758.0318	3634.0636	0.0000	1.0000
JC69IG	57	-1760.2706	3634.5412	0.0000	1.0000
JC69G	56	-1764.1463	3640.2926	0.0000	1.0000
F81I	59	-1762.5828	3643.1656	0.0000	1.0000
JC69I	56	-1768.6840	3649.3680	0.0000	1.0000
GTR	63	-1764.4523	3654.9046	0.0000	1.0000
SYM	60	-1773.7752	3667.5504	0.0000	1.0000
HKY	59	-1778.3576	3674.7152	0.0000	1.0000
K2P	56	-1783.9838	3679.9677	0.0000	1.0000
F81	58	-1815.7804	3747.5607	0.0000	1.0000
JC69	55	-1821.7810	3753.5620	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	58	-1718.0460	3795.2490	0.6830	0.6830
K2PG	57	-1722.0932	3797.1510	0.2639	0.9469
HKYIG	61	-1711.7330	3801.2001	0.0348	0.9817
HKYG	60	-1715.7286	3802.9990	0.0142	0.9959
K2PI	57	-1727.2560	3807.4768	0.0015	0.9974
GTRIG	65	-1702.6419	3807.7873	0.0013	0.9987
SYMIG	62	-1712.6964	3809.3193	0.0006	0.9993
GTRG	64	-1706.9840	3810.2792	0.0004	0.9997
SYMG	61	-1716.7519	3811.2380	0.0002	0.9999
HKYI	60	-1721.0373	3813.6162	0.0001	1.0000
GTRI	64	-1712.1096	3820.5305	0.0000	1.0000
SYMI	61	-1721.6853	3821.1048	0.0000	1.0000
JC69IG	57	-1760.2706	3873.5058	0.0000	1.0000
JC69G	56	-1764.1463	3875.0649	0.0000	1.0000
F81IG	60	-1754.2333	3880.0083	0.0000	1.0000
F81G	59	-1758.0318	3881.4130	0.0000	1.0000
JC69I	56	-1768.6840	3884.1403	0.0000	1.0000
F81I	59	-1762.5828	3890.5150	0.0000	1.0000
K2P	56	-1783.9838	3914.7400	0.0000	1.0000
GTR	63	-1764.4523	3919.0235	0.0000	1.0000
SYM	60	-1773.7752	3919.0921	0.0000	1.0000
HKY	59	-1778.3576	3922.0646	0.0000	1.0000
JC69	55	-1821.7810	3984.1420	0.0000	1.0000
F81	58	-1815.7804	3990.7177	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9114_hg18.phylip = 1 - 489;
 Partition Dummy = 1:261-aggressive-phylip/9114_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9114_hg18.phylip = 1 - 489;
 Partition Dummy = 1:261-aggressive-phylip/9114_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

