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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: HKYIG


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
HKYIG	61	-5065.5119	10263.2042	0.9750	0.9750
GTRIG	65	-5064.6680	10270.9462	0.0203	0.9954
SYMIG	62	-5069.7414	10274.0111	0.0044	0.9997
HKYG	60	-5075.0653	10279.9694	0.0002	1.0000
K2PG	57	-5080.9782	10284.8078	0.0000	1.0000
GTRG	64	-5073.6691	10286.5813	0.0000	1.0000
SYMG	61	-5079.6195	10291.4194	0.0000	1.0000
K2PIG	58	-5088.0874	10301.3490	0.0000	1.0000
HKYI	60	-5107.6549	10345.1485	0.0000	1.0000
K2PI	57	-5112.7782	10348.4077	0.0000	1.0000
SYMI	61	-5111.8840	10355.9484	0.0000	1.0000
GTRI	64	-5111.9637	10363.1707	0.0000	1.0000
F81IG	60	-5282.7784	10695.3955	0.0000	1.0000
JC69IG	57	-5286.6165	10696.0843	0.0000	1.0000
F81G	59	-5292.5694	10712.6421	0.0000	1.0000
JC69G	56	-5297.0900	10714.7148	0.0000	1.0000
F81I	59	-5318.0421	10763.5876	0.0000	1.0000
JC69I	56	-5328.0926	10776.7200	0.0000	1.0000
K2P	56	-5358.4446	10837.4239	0.0000	1.0000
HKY	59	-5356.1369	10839.7771	0.0000	1.0000
SYM	60	-5355.8538	10841.5462	0.0000	1.0000
GTR	63	-5354.5686	10846.0198	0.0000	1.0000
JC69	55	-5557.7812	11233.7868	0.0000	1.0000
F81	58	-5554.6236	11234.4215	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-5065.5119	10253.0238	0.9540	0.9540
GTRIG	65	-5064.6680	10259.3360	0.0406	0.9947
SYMIG	62	-5069.7414	10263.4828	0.0051	0.9998
HKYG	60	-5075.0653	10270.1307	0.0002	1.0000
GTRG	64	-5073.6691	10275.3381	0.0000	1.0000
K2PG	57	-5080.9782	10275.9564	0.0000	1.0000
SYMG	61	-5079.6195	10281.2390	0.0000	1.0000
K2PIG	58	-5088.0874	10292.1748	0.0000	1.0000
HKYI	60	-5107.6549	10335.3098	0.0000	1.0000
K2PI	57	-5112.7782	10339.5563	0.0000	1.0000
SYMI	61	-5111.8840	10345.7680	0.0000	1.0000
GTRI	64	-5111.9637	10351.9274	0.0000	1.0000
F81IG	60	-5282.7784	10685.5568	0.0000	1.0000
JC69IG	57	-5286.6165	10687.2329	0.0000	1.0000
F81G	59	-5292.5694	10703.1387	0.0000	1.0000
JC69G	56	-5297.0900	10706.1800	0.0000	1.0000
F81I	59	-5318.0421	10754.0843	0.0000	1.0000
JC69I	56	-5328.0926	10768.1852	0.0000	1.0000
K2P	56	-5358.4446	10828.8892	0.0000	1.0000
HKY	59	-5356.1369	10830.2737	0.0000	1.0000
SYM	60	-5355.8538	10831.7075	0.0000	1.0000
GTR	63	-5354.5686	10835.1372	0.0000	1.0000
F81	58	-5554.6236	11225.2473	0.0000	1.0000
JC69	55	-5557.7812	11225.5625	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-5065.5119	10539.1652	0.8874	0.8874
K2PG	57	-5080.9782	10543.3344	0.1104	0.9978
HKYG	60	-5075.0653	10551.5812	0.0018	0.9995
SYMIG	62	-5069.7414	10554.3151	0.0005	1.0000
GTRIG	65	-5064.6680	10564.2407	0.0000	1.0000
K2PIG	58	-5088.0874	10564.2436	0.0000	1.0000
SYMG	61	-5079.6195	10567.3804	0.0000	1.0000
GTRG	64	-5073.6691	10575.5520	0.0000	1.0000
K2PI	57	-5112.7782	10606.9343	0.0000	1.0000
HKYI	60	-5107.6549	10616.7603	0.0000	1.0000
SYMI	61	-5111.8840	10631.9094	0.0000	1.0000
GTRI	64	-5111.9637	10652.1413	0.0000	1.0000
JC69IG	57	-5286.6165	10954.6109	0.0000	1.0000
F81IG	60	-5282.7784	10967.0074	0.0000	1.0000
JC69G	56	-5297.0900	10968.8672	0.0000	1.0000
F81G	59	-5292.5694	10979.8984	0.0000	1.0000
F81I	59	-5318.0421	11030.8440	0.0000	1.0000
JC69I	56	-5328.0926	11030.8724	0.0000	1.0000
K2P	56	-5358.4446	11091.5763	0.0000	1.0000
HKY	59	-5356.1369	11107.0334	0.0000	1.0000
SYM	60	-5355.8538	11113.1580	0.0000	1.0000
GTR	63	-5354.5686	11130.6602	0.0000	1.0000
JC69	55	-5557.7812	11483.5588	0.0000	1.0000
F81	58	-5554.6236	11497.3161	0.0000	1.0000



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

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

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


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

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


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

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


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

