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

Input data from file "261-aggressive-phylip/7055_hg18.phylip" (29 taxa, 685 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	-2361.2503	4863.9199	0.5686	0.5686
GTRIG	65	-2360.3319	4864.5250	0.4202	0.9888
GTRI	64	-2365.1758	4871.7710	0.0112	1.0000
HKYG	60	-2386.7090	4905.1487	0.0000	1.0000
HKYIG	61	-2386.2569	4906.6550	0.0000	1.0000
HKYI	60	-2388.9100	4909.5508	0.0000	1.0000
SYMG	61	-2388.8006	4911.7425	0.0000	1.0000
SYMIG	62	-2388.3792	4913.3179	0.0000	1.0000
SYMI	61	-2390.6754	4915.4921	0.0000	1.0000
K2PG	57	-2406.0046	4936.5546	0.0000	1.0000
K2PIG	58	-2405.5687	4938.0704	0.0000	1.0000
K2PI	57	-2408.0788	4940.7030	0.0000	1.0000
GTR	63	-2405.9916	4950.9688	0.0000	1.0000
F81G	59	-2418.5928	4966.5136	0.0000	1.0000
F81IG	60	-2418.1005	4967.9317	0.0000	1.0000
F81I	59	-2420.9698	4971.2675	0.0000	1.0000
HKY	59	-2426.7635	4982.8550	0.0000	1.0000
JC69G	56	-2434.6046	4991.3748	0.0000	1.0000
SYM	60	-2430.2527	4992.2362	0.0000	1.0000
JC69IG	57	-2434.1555	4992.8564	0.0000	1.0000
JC69I	56	-2436.7737	4995.7130	0.0000	1.0000
K2P	56	-2446.2906	5014.7468	0.0000	1.0000
F81	58	-2457.8810	5042.6950	0.0000	1.0000
JC69	55	-2473.9410	5067.6754	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2361.2503	4850.5005	0.5151	0.5151
GTRIG	65	-2360.3319	4850.6639	0.4747	0.9898
GTRI	64	-2365.1758	4858.3516	0.0102	1.0000
HKYG	60	-2386.7090	4893.4179	0.0000	1.0000
HKYIG	61	-2386.2569	4894.5138	0.0000	1.0000
HKYI	60	-2388.9100	4897.8200	0.0000	1.0000
SYMG	61	-2388.8006	4899.6013	0.0000	1.0000
SYMIG	62	-2388.3792	4900.7584	0.0000	1.0000
SYMI	61	-2390.6754	4903.3509	0.0000	1.0000
K2PG	57	-2406.0046	4926.0091	0.0000	1.0000
K2PIG	58	-2405.5687	4927.1375	0.0000	1.0000
K2PI	57	-2408.0788	4930.1575	0.0000	1.0000
GTR	63	-2405.9916	4937.9833	0.0000	1.0000
F81G	59	-2418.5928	4955.1856	0.0000	1.0000
F81IG	60	-2418.1005	4956.2009	0.0000	1.0000
F81I	59	-2420.9698	4959.9395	0.0000	1.0000
HKY	59	-2426.7635	4971.5270	0.0000	1.0000
SYM	60	-2430.2527	4980.5054	0.0000	1.0000
JC69G	56	-2434.6046	4981.2092	0.0000	1.0000
JC69IG	57	-2434.1555	4982.3109	0.0000	1.0000
JC69I	56	-2436.7737	4985.5474	0.0000	1.0000
K2P	56	-2446.2906	5004.5812	0.0000	1.0000
F81	58	-2457.8810	5031.7621	0.0000	1.0000
JC69	55	-2473.9410	5057.8821	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2361.2503	5140.3833	0.8965	0.8965
GTRIG	65	-2360.3319	5145.0761	0.0858	0.9823
GTRI	64	-2365.1758	5148.2344	0.0177	1.0000
HKYG	60	-2386.7090	5165.1831	0.0000	1.0000
HKYI	60	-2388.9100	5169.5852	0.0000	1.0000
HKYIG	61	-2386.2569	5170.8083	0.0000	1.0000
SYMG	61	-2388.8006	5175.8958	0.0000	1.0000
SYMI	61	-2390.6754	5179.6454	0.0000	1.0000
SYMIG	62	-2388.3792	5181.5824	0.0000	1.0000
K2PG	57	-2406.0046	5184.1860	0.0000	1.0000
K2PI	57	-2408.0788	5188.3344	0.0000	1.0000
K2PIG	58	-2405.5687	5189.8438	0.0000	1.0000
F81G	59	-2418.5928	5222.4213	0.0000	1.0000
GTR	63	-2405.9916	5223.3367	0.0000	1.0000
F81I	59	-2420.9698	5227.1752	0.0000	1.0000
F81IG	60	-2418.1005	5227.9661	0.0000	1.0000
JC69G	56	-2434.6046	5234.8567	0.0000	1.0000
HKY	59	-2426.7635	5238.7627	0.0000	1.0000
JC69I	56	-2436.7737	5239.1949	0.0000	1.0000
JC69IG	57	-2434.1555	5240.4878	0.0000	1.0000
SYM	60	-2430.2527	5252.2706	0.0000	1.0000
K2P	56	-2446.2906	5258.2287	0.0000	1.0000
F81	58	-2457.8810	5294.4684	0.0000	1.0000
JC69	55	-2473.9410	5307.0001	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/7055_hg18.phylip = 1 - 685;
 Partition Dummy = 1:261-aggressive-phylip/7055_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/7055_hg18.phylip = 1 - 685;
 Partition Dummy = 1:261-aggressive-phylip/7055_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/7055_hg18.phylip = 1 - 685;
 Partition Dummy = 1:261-aggressive-phylip/7055_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

