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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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	-1784.7054	3717.3628	0.2916	0.2916
GTRIG	65	-1783.3859	3717.3968	0.2867	0.5783
HKYG	60	-1790.1551	3717.6974	0.2467	0.8249
HKYIG	61	-1789.1894	3718.3883	0.1746	0.9996
HKYI	60	-1796.7604	3730.9081	0.0003	0.9999
GTRI	64	-1792.6192	3733.1904	0.0001	1.0000
SYMIG	62	-1812.0048	3766.6539	0.0000	1.0000
SYMI	61	-1815.9219	3771.8534	0.0000	1.0000
K2PIG	58	-1820.0314	3772.2425	0.0000	1.0000
SYMG	61	-1821.8184	3783.6464	0.0000	1.0000
K2PI	57	-1829.5923	3788.7790	0.0000	1.0000
K2PG	57	-1830.5676	3790.7296	0.0000	1.0000
F81IG	60	-1830.8948	3799.1768	0.0000	1.0000
F81G	59	-1832.5681	3799.9134	0.0000	1.0000
HKY	59	-1836.6729	3808.1231	0.0000	1.0000
F81I	59	-1837.5658	3809.9088	0.0000	1.0000
GTR	63	-1834.1122	3813.5162	0.0000	1.0000
JC69IG	57	-1862.5811	3854.7565	0.0000	1.0000
JC69G	56	-1869.1662	3865.3535	0.0000	1.0000
SYM	60	-1867.0546	3871.4964	0.0000	1.0000
K2P	56	-1872.6205	3872.2622	0.0000	1.0000
JC69I	56	-1873.9303	3874.8818	0.0000	1.0000
F81	58	-1876.9797	3886.1392	0.0000	1.0000
JC69	55	-1908.6259	3941.7119	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1783.3859	3696.7718	0.4849	0.4849
GTRG	64	-1784.7054	3697.4107	0.3523	0.8372
HKYG	60	-1790.1551	3700.3102	0.0827	0.9199
HKYIG	61	-1789.1894	3700.3788	0.0799	0.9998
GTRI	64	-1792.6192	3713.2384	0.0001	0.9999
HKYI	60	-1796.7604	3713.5209	0.0001	1.0000
SYMIG	62	-1812.0048	3748.0095	0.0000	1.0000
SYMI	61	-1815.9219	3753.8439	0.0000	1.0000
K2PIG	58	-1820.0314	3756.0628	0.0000	1.0000
SYMG	61	-1821.8184	3765.6369	0.0000	1.0000
K2PI	57	-1829.5923	3773.1846	0.0000	1.0000
K2PG	57	-1830.5676	3775.1353	0.0000	1.0000
F81IG	60	-1830.8948	3781.7896	0.0000	1.0000
F81G	59	-1832.5681	3783.1361	0.0000	1.0000
HKY	59	-1836.6729	3791.3459	0.0000	1.0000
F81I	59	-1837.5658	3793.1315	0.0000	1.0000
GTR	63	-1834.1122	3794.2244	0.0000	1.0000
JC69IG	57	-1862.5811	3839.1621	0.0000	1.0000
JC69G	56	-1869.1662	3850.3323	0.0000	1.0000
SYM	60	-1867.0546	3854.1092	0.0000	1.0000
K2P	56	-1872.6205	3857.2410	0.0000	1.0000
JC69I	56	-1873.9303	3859.8606	0.0000	1.0000
F81	58	-1876.9797	3869.9595	0.0000	1.0000
JC69	55	-1908.6259	3927.2518	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1790.1551	3950.9869	0.8911	0.8911
HKYIG	61	-1789.1894	3955.2334	0.1066	0.9977
HKYI	60	-1796.7604	3964.1975	0.0012	0.9990
GTRG	64	-1784.7054	3964.7991	0.0009	0.9998
GTRIG	65	-1783.3859	3968.3381	0.0002	1.0000
GTRI	64	-1792.6192	3980.6268	0.0000	1.0000
K2PIG	58	-1820.0314	3998.3836	0.0000	1.0000
SYMIG	62	-1812.0048	4007.0421	0.0000	1.0000
SYMI	61	-1815.9219	4008.6985	0.0000	1.0000
K2PI	57	-1829.5923	4011.3275	0.0000	1.0000
K2PG	57	-1830.5676	4013.2781	0.0000	1.0000
SYMG	61	-1821.8184	4020.4915	0.0000	1.0000
F81G	59	-1832.5681	4029.6348	0.0000	1.0000
F81IG	60	-1830.8948	4032.4663	0.0000	1.0000
HKY	59	-1836.6729	4037.8446	0.0000	1.0000
F81I	59	-1837.5658	4039.6302	0.0000	1.0000
GTR	63	-1834.1122	4057.4349	0.0000	1.0000
JC69IG	57	-1862.5811	4077.3050	0.0000	1.0000
JC69G	56	-1869.1662	4084.2972	0.0000	1.0000
K2P	56	-1872.6205	4091.2059	0.0000	1.0000
JC69I	56	-1873.9303	4093.8255	0.0000	1.0000
SYM	60	-1867.0546	4104.7858	0.0000	1.0000
F81	58	-1876.9797	4112.2803	0.0000	1.0000
JC69	55	-1908.6259	4157.0388	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/11023_hg18.phylip = 1 - 482;
 Partition Dummy = 1:261-aggressive-phylip/11023_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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11023_hg18.phylip = 1 - 482;
 Partition Dummy = 1:261-aggressive-phylip/11023_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 (HKYG)]

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


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

