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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: K2PG


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
SYMG	61	-2394.7739	4921.7281	0.6449	0.6449
GTRG	64	-2392.1588	4923.5609	0.2579	0.9028
SYMIG	62	-2395.8116	4926.1514	0.0706	0.9734
GTRIG	65	-2393.3585	4928.3272	0.0238	0.9972
K2PG	57	-2405.0280	4932.9075	0.0024	0.9996
K2PIG	58	-2405.9901	4937.1544	0.0003	0.9999
HKYG	60	-2404.7138	4939.2663	0.0001	1.0000
HKYIG	61	-2405.6773	4943.5348	0.0000	1.0000
JC69G	56	-2422.5901	4965.7150	0.0000	1.0000
F81G	59	-2422.3985	4972.3005	0.0000	1.0000
JC69IG	57	-2426.8471	4976.5455	0.0000	1.0000
F81IG	60	-2426.6182	4983.0752	0.0000	1.0000
SYM	60	-2443.3495	5016.5376	0.0000	1.0000
GTR	63	-2441.6479	5020.1784	0.0000	1.0000
SYMI	61	-2446.4308	5025.0420	0.0000	1.0000
K2P	56	-2453.0989	5026.7325	0.0000	1.0000
GTRI	64	-2444.8501	5028.9434	0.0000	1.0000
HKY	59	-2452.7929	5033.0892	0.0000	1.0000
K2PI	57	-2456.0425	5034.9364	0.0000	1.0000
HKYI	60	-2455.7517	5041.3422	0.0000	1.0000
JC69	55	-2469.3942	5057.0127	0.0000	1.0000
F81	58	-2469.2357	5063.6456	0.0000	1.0000
JC69I	56	-2476.6562	5073.8472	0.0000	1.0000
F81I	59	-2476.4686	5080.4407	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	61	-2394.7739	4911.5478	0.5295	0.5295
GTRG	64	-2392.1588	4912.3176	0.3603	0.8898
SYMIG	62	-2395.8116	4915.6231	0.0690	0.9588
GTRIG	65	-2393.3585	4916.7169	0.0399	0.9988
K2PG	57	-2405.0280	4924.0561	0.0010	0.9998
K2PIG	58	-2405.9901	4927.9802	0.0001	0.9999
HKYG	60	-2404.7138	4929.4276	0.0001	1.0000
HKYIG	61	-2405.6773	4933.3545	0.0000	1.0000
JC69G	56	-2422.5901	4957.1802	0.0000	1.0000
F81G	59	-2422.3985	4962.7971	0.0000	1.0000
JC69IG	57	-2426.8471	4967.6941	0.0000	1.0000
F81IG	60	-2426.6182	4973.2365	0.0000	1.0000
SYM	60	-2443.3495	5006.6989	0.0000	1.0000
GTR	63	-2441.6479	5009.2958	0.0000	1.0000
SYMI	61	-2446.4308	5014.8616	0.0000	1.0000
GTRI	64	-2444.8501	5017.7001	0.0000	1.0000
K2P	56	-2453.0989	5018.1978	0.0000	1.0000
HKY	59	-2452.7929	5023.5858	0.0000	1.0000
K2PI	57	-2456.0425	5026.0850	0.0000	1.0000
HKYI	60	-2455.7517	5031.5035	0.0000	1.0000
JC69	55	-2469.3942	5048.7884	0.0000	1.0000
F81	58	-2469.2357	5054.4713	0.0000	1.0000
JC69I	56	-2476.6562	5065.3125	0.0000	1.0000
F81I	59	-2476.4686	5070.9373	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2405.0280	5191.4341	0.9452	0.9452
SYMG	61	-2394.7739	5197.6891	0.0414	0.9867
K2PIG	58	-2405.9901	5200.0490	0.0127	0.9994
SYMIG	62	-2395.8116	5206.4554	0.0005	0.9999
HKYG	60	-2404.7138	5210.8782	0.0001	1.0000
GTRG	64	-2392.1588	5212.5315	0.0000	1.0000
HKYIG	61	-2405.6773	5219.4959	0.0000	1.0000
JC69G	56	-2422.5901	5219.8674	0.0000	1.0000
GTRIG	65	-2393.3585	5221.6216	0.0000	1.0000
JC69IG	57	-2426.8471	5235.0721	0.0000	1.0000
F81G	59	-2422.3985	5239.5568	0.0000	1.0000
F81IG	60	-2426.6182	5254.6870	0.0000	1.0000
K2P	56	-2453.0989	5280.8849	0.0000	1.0000
SYM	60	-2443.3495	5288.1495	0.0000	1.0000
K2PI	57	-2456.0425	5293.4630	0.0000	1.0000
HKY	59	-2452.7929	5300.3455	0.0000	1.0000
SYMI	61	-2446.4308	5301.0030	0.0000	1.0000
GTR	63	-2441.6479	5304.8188	0.0000	1.0000
JC69	55	-2469.3942	5306.7847	0.0000	1.0000
HKYI	60	-2455.7517	5312.9540	0.0000	1.0000
GTRI	64	-2444.8501	5317.9140	0.0000	1.0000
F81	58	-2469.2357	5326.5402	0.0000	1.0000
JC69I	56	-2476.6562	5327.9996	0.0000	1.0000
F81I	59	-2476.4686	5347.6970	0.0000	1.0000



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

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

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


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

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


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

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


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

