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

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

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


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	-5132.9029	10400.6199	0.6818	0.6818
GTRIG	65	-5132.5554	10402.1437	0.3182	1.0000
HKYG	60	-5155.2305	10436.4366	0.0000	1.0000
SYMG	61	-5154.2958	10436.7714	0.0000	1.0000
HKYIG	61	-5155.1993	10438.5784	0.0000	1.0000
SYMIG	62	-5154.2515	10438.8905	0.0000	1.0000
K2PG	57	-5173.6825	10466.7493	0.0000	1.0000
K2PIG	58	-5173.6835	10468.9447	0.0000	1.0000
GTRI	64	-5208.1587	10551.1315	0.0000	1.0000
HKYI	60	-5226.6860	10579.3474	0.0000	1.0000
SYMI	61	-5228.2539	10584.6876	0.0000	1.0000
K2PI	57	-5245.1807	10609.7458	0.0000	1.0000
F81G	59	-5299.1457	10722.0663	0.0000	1.0000
F81IG	60	-5299.1472	10724.2699	0.0000	1.0000
JC69G	56	-5316.1371	10749.4687	0.0000	1.0000
JC69IG	57	-5316.1386	10751.6615	0.0000	1.0000
F81I	59	-5366.8811	10857.5371	0.0000	1.0000
JC69I	56	-5383.3098	10883.8142	0.0000	1.0000
GTR	63	-5394.1059	10920.8108	0.0000	1.0000
HKY	59	-5411.2533	10946.2816	0.0000	1.0000
SYM	60	-5411.0977	10948.1710	0.0000	1.0000
K2P	56	-5428.6526	10974.4997	0.0000	1.0000
F81	58	-5544.0231	11209.6240	0.0000	1.0000
JC69	55	-5558.7731	11232.5543	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-5132.9029	10393.8059	0.6576	0.6576
GTRIG	65	-5132.5554	10395.1109	0.3424	1.0000
HKYG	60	-5155.2305	10430.4611	0.0000	1.0000
SYMG	61	-5154.2958	10430.5916	0.0000	1.0000
HKYIG	61	-5155.1993	10432.3986	0.0000	1.0000
SYMIG	62	-5154.2515	10432.5029	0.0000	1.0000
K2PG	57	-5173.6825	10461.3649	0.0000	1.0000
K2PIG	58	-5173.6835	10463.3669	0.0000	1.0000
GTRI	64	-5208.1587	10544.3174	0.0000	1.0000
HKYI	60	-5226.6860	10573.3719	0.0000	1.0000
SYMI	61	-5228.2539	10578.5079	0.0000	1.0000
K2PI	57	-5245.1807	10604.3614	0.0000	1.0000
F81G	59	-5299.1457	10716.2914	0.0000	1.0000
F81IG	60	-5299.1472	10718.2944	0.0000	1.0000
JC69G	56	-5316.1371	10744.2742	0.0000	1.0000
JC69IG	57	-5316.1386	10746.2772	0.0000	1.0000
F81I	59	-5366.8811	10851.7622	0.0000	1.0000
JC69I	56	-5383.3098	10878.6197	0.0000	1.0000
GTR	63	-5394.1059	10914.2118	0.0000	1.0000
HKY	59	-5411.2533	10940.5067	0.0000	1.0000
SYM	60	-5411.0977	10942.1955	0.0000	1.0000
K2P	56	-5428.6526	10969.3053	0.0000	1.0000
F81	58	-5544.0231	11204.0462	0.0000	1.0000
JC69	55	-5558.7731	11227.5462	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-5132.9029	10724.0005	0.9617	0.9617
GTRIG	65	-5132.5554	10730.4649	0.0380	0.9996
HKYG	60	-5155.2305	10740.0186	0.0003	1.0000
SYMG	61	-5154.2958	10745.3084	0.0000	1.0000
HKYIG	61	-5155.1993	10747.1154	0.0000	1.0000
SYMIG	62	-5154.2515	10752.3790	0.0000	1.0000
K2PG	57	-5173.6825	10755.4446	0.0000	1.0000
K2PIG	58	-5173.6835	10762.6058	0.0000	1.0000
GTRI	64	-5208.1587	10874.5121	0.0000	1.0000
HKYI	60	-5226.6860	10882.9294	0.0000	1.0000
SYMI	61	-5228.2539	10893.2247	0.0000	1.0000
K2PI	57	-5245.1807	10898.4411	0.0000	1.0000
F81G	59	-5299.1457	11020.6896	0.0000	1.0000
F81IG	60	-5299.1472	11027.8519	0.0000	1.0000
JC69G	56	-5316.1371	11033.1946	0.0000	1.0000
JC69IG	57	-5316.1386	11040.3568	0.0000	1.0000
F81I	59	-5366.8811	11156.1605	0.0000	1.0000
JC69I	56	-5383.3098	11167.5400	0.0000	1.0000
GTR	63	-5394.1059	11239.2472	0.0000	1.0000
HKY	59	-5411.2533	11244.9049	0.0000	1.0000
SYM	60	-5411.0977	11251.7530	0.0000	1.0000
K2P	56	-5428.6526	11258.2256	0.0000	1.0000
F81	58	-5544.0231	11503.2851	0.0000	1.0000
JC69	55	-5558.7731	11511.3073	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/524_hg18.phylip = 1 - 1286;
 Partition Dummy = 1:261-aggressive-phylip/524_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/524_hg18.phylip = 1 - 1286;
 Partition Dummy = 1:261-aggressive-phylip/524_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 (GTRG)]

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


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

