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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-4413.0609	8968.2920	0.7073	0.7073
GTRG	64	-4415.1360	8970.0568	0.2927	1.0000
HKYIG	61	-4431.2111	8995.0908	0.0000	1.0000
HKYG	60	-4433.0514	8996.4127	0.0000	1.0000
SYMG	61	-4445.3835	9023.4356	0.0000	1.0000
SYMIG	62	-4446.2285	9027.4909	0.0000	1.0000
K2PG	57	-4452.8075	9028.8884	0.0000	1.0000
K2PIG	58	-4452.1746	9029.9615	0.0000	1.0000
GTRI	64	-4451.5866	9042.9578	0.0000	1.0000
HKYI	60	-4460.9965	9052.3029	0.0000	1.0000
SYMI	61	-4473.3100	9079.2886	0.0000	1.0000
K2PI	57	-4482.5695	9088.4125	0.0000	1.0000
F81G	59	-4553.1784	9234.3146	0.0000	1.0000
F81IG	60	-4558.8120	9247.9338	0.0000	1.0000
JC69G	56	-4567.5105	9255.9621	0.0000	1.0000
JC69IG	57	-4574.2067	9271.6868	0.0000	1.0000
GTR	63	-4578.5784	9294.5627	0.0000	1.0000
F81I	59	-4590.2699	9308.4976	0.0000	1.0000
HKY	59	-4594.1414	9316.2406	0.0000	1.0000
JC69I	56	-4601.6184	9324.1780	0.0000	1.0000
SYM	60	-4603.7912	9337.8922	0.0000	1.0000
K2P	56	-4612.3590	9345.6591	0.0000	1.0000
F81	58	-4717.3934	9560.3991	0.0000	1.0000
JC69	55	-4728.0169	9574.6491	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4413.0609	8956.1218	0.7456	0.7456
GTRG	64	-4415.1360	8958.2721	0.2544	1.0000
HKYIG	61	-4431.2111	8984.4222	0.0000	1.0000
HKYG	60	-4433.0514	8986.1029	0.0000	1.0000
SYMG	61	-4445.3835	9012.7670	0.0000	1.0000
SYMIG	62	-4446.2285	9016.4570	0.0000	1.0000
K2PG	57	-4452.8075	9019.6149	0.0000	1.0000
K2PIG	58	-4452.1746	9020.3491	0.0000	1.0000
GTRI	64	-4451.5866	9031.1731	0.0000	1.0000
HKYI	60	-4460.9965	9041.9931	0.0000	1.0000
SYMI	61	-4473.3100	9068.6200	0.0000	1.0000
K2PI	57	-4482.5695	9079.1390	0.0000	1.0000
F81G	59	-4553.1784	9224.3568	0.0000	1.0000
F81IG	60	-4558.8120	9237.6239	0.0000	1.0000
JC69G	56	-4567.5105	9247.0209	0.0000	1.0000
JC69IG	57	-4574.2067	9262.4133	0.0000	1.0000
GTR	63	-4578.5784	9283.1567	0.0000	1.0000
F81I	59	-4590.2699	9298.5398	0.0000	1.0000
HKY	59	-4594.1414	9306.2828	0.0000	1.0000
JC69I	56	-4601.6184	9315.2368	0.0000	1.0000
SYM	60	-4603.7912	9327.5823	0.0000	1.0000
K2P	56	-4612.3590	9336.7179	0.0000	1.0000
F81	58	-4717.3934	9550.7868	0.0000	1.0000
JC69	55	-4728.0169	9566.0338	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-4415.1360	9255.7241	0.7698	0.7698
GTRIG	65	-4413.0609	9258.2215	0.2209	0.9907
HKYG	60	-4433.0514	9264.9642	0.0076	0.9983
HKYIG	61	-4431.2111	9267.9312	0.0017	1.0000
K2PG	57	-4452.8075	9284.5332	0.0000	1.0000
K2PIG	58	-4452.1746	9289.9151	0.0000	1.0000
SYMG	61	-4445.3835	9296.2760	0.0000	1.0000
SYMIG	62	-4446.2285	9304.6137	0.0000	1.0000
HKYI	60	-4460.9965	9320.8544	0.0000	1.0000
GTRI	64	-4451.5866	9328.6252	0.0000	1.0000
K2PI	57	-4482.5695	9344.0573	0.0000	1.0000
SYMI	61	-4473.3100	9352.1290	0.0000	1.0000
F81G	59	-4553.1784	9498.5704	0.0000	1.0000
JC69G	56	-4567.5105	9507.2915	0.0000	1.0000
F81IG	60	-4558.8120	9516.4852	0.0000	1.0000
JC69IG	57	-4574.2067	9527.3316	0.0000	1.0000
F81I	59	-4590.2699	9572.7535	0.0000	1.0000
JC69I	56	-4601.6184	9575.5073	0.0000	1.0000
GTR	63	-4578.5784	9575.9611	0.0000	1.0000
HKY	59	-4594.1414	9580.4964	0.0000	1.0000
K2P	56	-4612.3590	9596.9885	0.0000	1.0000
SYM	60	-4603.7912	9606.4436	0.0000	1.0000
F81	58	-4717.3934	9820.3527	0.0000	1.0000
JC69	55	-4728.0169	9821.6566	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/9789_hg18.phylip = 1 - 771;
 Partition Dummy = 1:261-aggressive-phylip/9789_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 (GTRIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9789_hg18.phylip = 1 - 771;
 Partition Dummy = 1:261-aggressive-phylip/9789_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 BIC model (GTRG)]

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

