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

Input data from file "261-aggressive-phylip/9589_hg18.phylip" (29 taxa, 422 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	-1620.5010	3392.3074	0.5907	0.5907
GTRIG	65	-1619.5338	3393.1688	0.3840	0.9747
GTRI	64	-1624.9065	3401.1183	0.0072	0.9820
K2PG	57	-1634.6648	3401.4944	0.0060	0.9879
K2PIG	58	-1633.9366	3402.7273	0.0032	0.9912
SYMG	61	-1629.9237	3402.8586	0.0030	0.9942
HKYG	60	-1631.5326	3403.3421	0.0024	0.9966
SYMIG	62	-1629.0004	3403.7613	0.0019	0.9985
HKYIG	61	-1630.8939	3404.7990	0.0011	0.9996
K2PI	57	-1637.9949	3408.1546	0.0002	0.9998
HKYI	60	-1634.6858	3409.6487	0.0001	0.9999
SYMI	61	-1633.8585	3410.7281	0.0001	1.0000
F81G	59	-1655.8878	3449.3336	0.0000	1.0000
F81IG	60	-1654.9463	3450.1696	0.0000	1.0000
JC69G	56	-1661.7926	3453.0756	0.0000	1.0000
JC69IG	57	-1660.8738	3453.9125	0.0000	1.0000
F81I	59	-1659.3357	3456.2294	0.0000	1.0000
JC69I	56	-1665.2294	3459.9491	0.0000	1.0000
GTR	63	-1655.9283	3460.3817	0.0000	1.0000
SYM	60	-1665.7841	3471.8452	0.0000	1.0000
HKY	59	-1667.9620	3473.4820	0.0000	1.0000
K2P	56	-1672.6363	3474.7630	0.0000	1.0000
F81	58	-1693.0955	3521.0449	0.0000	1.0000
JC69	55	-1699.2197	3525.2701	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1620.5010	3369.0020	0.5035	0.5035
GTRIG	65	-1619.5338	3369.0676	0.4873	0.9908
GTRI	64	-1624.9065	3377.8130	0.0061	0.9969
SYMG	61	-1629.9237	3381.8475	0.0008	0.9978
SYMIG	62	-1629.0004	3382.0009	0.0008	0.9985
HKYG	60	-1631.5326	3383.0651	0.0004	0.9990
K2PG	57	-1634.6648	3383.3295	0.0004	0.9993
HKYIG	61	-1630.8939	3383.7879	0.0003	0.9997
K2PIG	58	-1633.9366	3383.8733	0.0003	1.0000
HKYI	60	-1634.6858	3389.3717	0.0000	1.0000
SYMI	61	-1633.8585	3389.7170	0.0000	1.0000
K2PI	57	-1637.9949	3389.9897	0.0000	1.0000
F81G	59	-1655.8878	3429.7756	0.0000	1.0000
F81IG	60	-1654.9463	3429.8926	0.0000	1.0000
JC69G	56	-1661.7926	3435.5851	0.0000	1.0000
JC69IG	57	-1660.8738	3435.7477	0.0000	1.0000
F81I	59	-1659.3357	3436.6714	0.0000	1.0000
GTR	63	-1655.9283	3437.8566	0.0000	1.0000
JC69I	56	-1665.2294	3442.4587	0.0000	1.0000
SYM	60	-1665.7841	3451.5682	0.0000	1.0000
HKY	59	-1667.9620	3453.9240	0.0000	1.0000
K2P	56	-1672.6363	3457.2726	0.0000	1.0000
F81	58	-1693.0955	3502.1909	0.0000	1.0000
JC69	55	-1699.2197	3508.4395	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1634.6648	3613.8948	0.8761	0.8761
K2PIG	58	-1633.9366	3618.4836	0.0883	0.9645
K2PI	57	-1637.9949	3620.5550	0.0314	0.9958
HKYG	60	-1631.5326	3625.7654	0.0023	0.9981
GTRG	64	-1620.5010	3627.8824	0.0008	0.9989
SYMG	61	-1629.9237	3628.5928	0.0006	0.9995
HKYIG	61	-1630.8939	3630.5332	0.0002	0.9997
GTRIG	65	-1619.5338	3631.9930	0.0001	0.9998
HKYI	60	-1634.6858	3632.0720	0.0001	0.9999
SYMIG	62	-1629.0004	3632.7912	0.0001	1.0000
SYMI	61	-1633.8585	3636.4623	0.0000	1.0000
GTRI	64	-1624.9065	3636.6933	0.0000	1.0000
JC69G	56	-1661.7926	3662.1054	0.0000	1.0000
JC69IG	57	-1660.8738	3666.3130	0.0000	1.0000
F81G	59	-1655.8878	3668.4309	0.0000	1.0000
JC69I	56	-1665.2294	3668.9790	0.0000	1.0000
F81IG	60	-1654.9463	3672.5929	0.0000	1.0000
F81I	59	-1659.3357	3675.3267	0.0000	1.0000
K2P	56	-1672.6363	3683.7929	0.0000	1.0000
HKY	59	-1667.9620	3692.5793	0.0000	1.0000
GTR	63	-1655.9283	3692.6919	0.0000	1.0000
SYM	60	-1665.7841	3694.2685	0.0000	1.0000
JC69	55	-1699.2197	3730.9148	0.0000	1.0000
F81	58	-1693.0955	3736.8012	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/9589_hg18.phylip = 1 - 422;
 Partition Dummy = 1:261-aggressive-phylip/9589_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/9589_hg18.phylip = 1 - 422;
 Partition Dummy = 1:261-aggressive-phylip/9589_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/9589_hg18.phylip = 1 - 422;
 Partition Dummy = 1:261-aggressive-phylip/9589_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

