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

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

Minimum AIC  model: GTRG
Minimum AICc model: SYMG
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
SYMG	61	-1789.3135	3718.4247	0.4704	0.4704
GTRG	64	-1786.1815	3720.0786	0.2057	0.6761
SYMIG	62	-1789.1470	3720.7185	0.1494	0.8255
GTRIG	65	-1785.2495	3720.8791	0.1379	0.9633
K2PG	57	-1797.4411	3724.2947	0.0250	0.9883
K2PIG	58	-1797.6884	3727.3675	0.0054	0.9937
HKYG	60	-1795.9783	3729.1397	0.0022	0.9959
SYMI	61	-1794.8708	3729.5392	0.0018	0.9977
HKYIG	61	-1795.1717	3730.1410	0.0013	0.9991
GTRI	64	-1791.6996	3731.1148	0.0008	0.9999
K2PI	57	-1803.2439	3735.9003	0.0001	1.0000
HKYI	60	-1801.2069	3739.5969	0.0000	1.0000
JC69G	56	-1826.9897	3780.8259	0.0000	1.0000
SYM	60	-1822.4985	3782.1801	0.0000	1.0000
JC69IG	57	-1826.5877	3782.5881	0.0000	1.0000
GTR	63	-1820.0088	3785.0815	0.0000	1.0000
F81G	59	-1826.0355	3786.6518	0.0000	1.0000
F81IG	60	-1825.3714	3787.9260	0.0000	1.0000
K2P	56	-1830.7732	3788.3930	0.0000	1.0000
JC69I	56	-1831.9823	3790.8111	0.0000	1.0000
HKY	59	-1829.4043	3793.3895	0.0000	1.0000
F81I	59	-1830.7848	3796.1504	0.0000	1.0000
JC69	55	-1859.3520	3842.9963	0.0000	1.0000
F81	58	-1858.4068	3848.8043	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1786.1815	3700.3629	0.3103	0.3103
GTRIG	65	-1785.2495	3700.4991	0.2899	0.6002
SYMG	61	-1789.3135	3700.6271	0.2719	0.8722
SYMIG	62	-1789.1470	3702.2940	0.1182	0.9903
K2PG	57	-1797.4411	3708.8821	0.0044	0.9947
K2PIG	58	-1797.6884	3711.3769	0.0013	0.9960
GTRI	64	-1791.6996	3711.3992	0.0012	0.9972
SYMI	61	-1794.8708	3711.7416	0.0010	0.9983
HKYG	60	-1795.9783	3711.9566	0.0009	0.9992
HKYIG	61	-1795.1717	3712.3433	0.0008	1.0000
K2PI	57	-1803.2439	3720.4877	0.0000	1.0000
HKYI	60	-1801.2069	3722.4138	0.0000	1.0000
SYM	60	-1822.4985	3764.9970	0.0000	1.0000
JC69G	56	-1826.9897	3765.9794	0.0000	1.0000
GTR	63	-1820.0088	3766.0177	0.0000	1.0000
JC69IG	57	-1826.5877	3767.1755	0.0000	1.0000
F81G	59	-1826.0355	3770.0710	0.0000	1.0000
F81IG	60	-1825.3714	3770.7429	0.0000	1.0000
K2P	56	-1830.7732	3773.5464	0.0000	1.0000
JC69I	56	-1831.9823	3775.9646	0.0000	1.0000
HKY	59	-1829.4043	3776.8087	0.0000	1.0000
F81I	59	-1830.7848	3779.5696	0.0000	1.0000
JC69	55	-1859.3520	3828.7039	0.0000	1.0000
F81	58	-1858.4068	3832.8136	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1797.4411	3947.6132	0.9488	0.9488
K2PIG	58	-1797.6884	3954.2962	0.0336	0.9824
SYMG	61	-1789.3135	3956.1112	0.0135	0.9959
K2PI	57	-1803.2439	3959.2188	0.0029	0.9988
SYMIG	62	-1789.1470	3961.9664	0.0007	0.9995
HKYG	60	-1795.9783	3963.2525	0.0004	0.9999
SYMI	61	-1794.8708	3967.2257	0.0001	0.9999
HKYIG	61	-1795.1717	3967.8274	0.0000	1.0000
GTRG	64	-1786.1815	3968.4118	0.0000	1.0000
GTRIG	65	-1785.2495	3972.7362	0.0000	1.0000
HKYI	60	-1801.2069	3973.7097	0.0000	1.0000
GTRI	64	-1791.6996	3979.4481	0.0000	1.0000
JC69G	56	-1826.9897	4000.5222	0.0000	1.0000
JC69IG	57	-1826.5877	4005.9066	0.0000	1.0000
K2P	56	-1830.7732	4008.0892	0.0000	1.0000
JC69I	56	-1831.9823	4010.5074	0.0000	1.0000
SYM	60	-1822.4985	4016.2928	0.0000	1.0000
F81G	59	-1826.0355	4017.1785	0.0000	1.0000
F81IG	60	-1825.3714	4022.0387	0.0000	1.0000
HKY	59	-1829.4043	4023.9163	0.0000	1.0000
F81I	59	-1830.7848	4026.6772	0.0000	1.0000
GTR	63	-1820.0088	4029.8783	0.0000	1.0000
JC69	55	-1859.3520	4059.0584	0.0000	1.0000
F81	58	-1858.4068	4075.7329	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/11900_hg18.phylip = 1 - 487;
 Partition Dummy = 1:261-aggressive-phylip/11900_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 (SYMG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11900_hg18.phylip = 1 - 487;
 Partition Dummy = 1:261-aggressive-phylip/11900_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=Fixed(Equal);
END;


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

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

