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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: SYMG


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	-2356.6665	4860.3719	0.9590	0.9590
GTRIG	65	-2358.8633	4867.4055	0.0285	0.9875
SYMG	61	-2365.0347	4869.2602	0.0113	0.9987
HKYG	60	-2368.9672	4874.5330	0.0008	0.9995
SYMIG	62	-2367.3744	4876.5438	0.0003	0.9998
HKYIG	61	-2369.2182	4877.6274	0.0002	1.0000
K2PIG	58	-2376.4479	4884.3450	0.0000	1.0000
K2PG	57	-2379.2571	4887.4061	0.0000	1.0000
GTRI	64	-2371.2485	4889.5360	0.0000	1.0000
SYMI	61	-2378.9920	4897.1749	0.0000	1.0000
HKYI	60	-2382.6060	4901.8107	0.0000	1.0000
K2PI	57	-2392.8994	4914.6908	0.0000	1.0000
F81G	59	-2423.8408	4981.6997	0.0000	1.0000
F81IG	60	-2422.7287	4982.0561	0.0000	1.0000
JC69G	56	-2432.1290	4990.6041	0.0000	1.0000
JC69IG	57	-2431.1129	4991.1177	0.0000	1.0000
GTR	63	-2426.4827	4997.3764	0.0000	1.0000
SYM	60	-2432.4215	5001.4416	0.0000	1.0000
F81I	59	-2435.2034	5004.4250	0.0000	1.0000
JC69I	56	-2443.2452	5012.8365	0.0000	1.0000
HKY	59	-2442.2103	5018.4388	0.0000	1.0000
K2P	56	-2451.1465	5028.6392	0.0000	1.0000
F81	58	-2495.5505	5122.5502	0.0000	1.0000
JC69	55	-2502.1813	5128.1742	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2356.6665	4841.3330	0.9560	0.9560
GTRIG	65	-2358.8633	4847.7266	0.0391	0.9951
SYMG	61	-2365.0347	4852.0693	0.0045	0.9995
HKYG	60	-2368.9672	4857.9344	0.0002	0.9998
SYMIG	62	-2367.3744	4858.7488	0.0002	0.9999
HKYIG	61	-2369.2182	4860.4365	0.0001	1.0000
K2PIG	58	-2376.4479	4868.8958	0.0000	1.0000
GTRI	64	-2371.2485	4870.4971	0.0000	1.0000
K2PG	57	-2379.2571	4872.5142	0.0000	1.0000
SYMI	61	-2378.9920	4879.9840	0.0000	1.0000
HKYI	60	-2382.6060	4885.2120	0.0000	1.0000
K2PI	57	-2392.8994	4899.7989	0.0000	1.0000
F81IG	60	-2422.7287	4965.4575	0.0000	1.0000
F81G	59	-2423.8408	4965.6816	0.0000	1.0000
JC69IG	57	-2431.1129	4976.2258	0.0000	1.0000
JC69G	56	-2432.1290	4976.2580	0.0000	1.0000
GTR	63	-2426.4827	4978.9655	0.0000	1.0000
SYM	60	-2432.4215	4984.8430	0.0000	1.0000
F81I	59	-2435.2034	4988.4069	0.0000	1.0000
JC69I	56	-2443.2452	4998.4904	0.0000	1.0000
HKY	59	-2442.2103	5002.4207	0.0000	1.0000
K2P	56	-2451.1465	5014.2931	0.0000	1.0000
F81	58	-2495.5505	5107.1010	0.0000	1.0000
JC69	55	-2502.1813	5114.3626	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	61	-2365.0347	5109.4039	0.4683	0.4683
HKYG	60	-2368.9672	5111.0504	0.2056	0.6739
GTRG	64	-2356.6665	5111.3234	0.1793	0.8532
K2PG	57	-2379.2571	5112.9744	0.0786	0.9318
K2PIG	58	-2376.4479	5113.5746	0.0582	0.9900
HKYIG	61	-2369.2182	5117.7711	0.0071	0.9971
SYMIG	62	-2367.3744	5120.3020	0.0020	0.9991
GTRIG	65	-2358.8633	5121.9356	0.0009	1.0000
SYMI	61	-2378.9920	5137.3186	0.0000	1.0000
HKYI	60	-2382.6060	5138.3280	0.0000	1.0000
K2PI	57	-2392.8994	5140.2591	0.0000	1.0000
GTRI	64	-2371.2485	5140.4875	0.0000	1.0000
JC69G	56	-2432.1290	5212.4996	0.0000	1.0000
F81G	59	-2423.8408	5214.5790	0.0000	1.0000
JC69IG	57	-2431.1129	5216.6860	0.0000	1.0000
F81IG	60	-2422.7287	5218.5735	0.0000	1.0000
JC69I	56	-2443.2452	5234.7320	0.0000	1.0000
F81I	59	-2435.2034	5237.3043	0.0000	1.0000
SYM	60	-2432.4215	5237.9590	0.0000	1.0000
GTR	63	-2426.4827	5244.7373	0.0000	1.0000
K2P	56	-2451.1465	5250.5347	0.0000	1.0000
HKY	59	-2442.2103	5251.3181	0.0000	1.0000
JC69	55	-2502.1813	5346.3856	0.0000	1.0000
F81	58	-2495.5505	5351.7798	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/11595_hg18.phylip = 1 - 502;
 Partition Dummy = 1:261-aggressive-phylip/11595_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/11595_hg18.phylip = 1 - 502;
 Partition Dummy = 1:261-aggressive-phylip/11595_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 (SYMG)]

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


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

