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

Input data from file "261-aggressive-phylip/4976_hg18.phylip" (29 taxa, 539 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	-2200.2512	4546.0551	0.7360	0.7360
GTRIG	65	-2199.9966	4548.1327	0.2604	0.9964
GTRI	64	-2205.5726	4556.6979	0.0036	1.0000
HKYG	60	-2220.6571	4576.6280	0.0000	1.0000
HKYIG	61	-2220.4298	4578.7170	0.0000	1.0000
HKYI	60	-2225.6286	4586.5710	0.0000	1.0000
GTR	63	-2237.7368	4618.4504	0.0000	1.0000
SYMG	61	-2241.9249	4621.7072	0.0000	1.0000
K2PG	57	-2250.5460	4628.8384	0.0000	1.0000
SYMIG	62	-2244.6724	4629.7565	0.0000	1.0000
SYMI	61	-2247.2160	4632.2895	0.0000	1.0000
K2PIG	58	-2253.4091	4637.0766	0.0000	1.0000
K2PI	57	-2256.1907	4640.1277	0.0000	1.0000
HKY	59	-2256.0599	4644.9006	0.0000	1.0000
SYM	60	-2278.0928	4691.4995	0.0000	1.0000
K2P	56	-2287.2980	4699.8407	0.0000	1.0000
F81G	59	-2294.0380	4720.8567	0.0000	1.0000
F81IG	60	-2293.8054	4722.9245	0.0000	1.0000
F81I	59	-2299.1930	4731.1667	0.0000	1.0000
JC69G	56	-2315.1408	4755.5264	0.0000	1.0000
JC69IG	57	-2318.7441	4765.2346	0.0000	1.0000
JC69I	56	-2320.5723	4766.3894	0.0000	1.0000
F81	58	-2328.8977	4788.0538	0.0000	1.0000
JC69	55	-2350.8482	4824.4501	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2200.2512	4528.5023	0.6759	0.6759
GTRIG	65	-2199.9966	4529.9932	0.3208	0.9967
GTRI	64	-2205.5726	4539.1452	0.0033	1.0000
HKYG	60	-2220.6571	4561.3142	0.0000	1.0000
HKYIG	61	-2220.4298	4562.8596	0.0000	1.0000
HKYI	60	-2225.6286	4571.2572	0.0000	1.0000
GTR	63	-2237.7368	4601.4736	0.0000	1.0000
SYMG	61	-2241.9249	4605.8498	0.0000	1.0000
SYMIG	62	-2244.6724	4613.3448	0.0000	1.0000
K2PG	57	-2250.5460	4615.0920	0.0000	1.0000
SYMI	61	-2247.2160	4616.4320	0.0000	1.0000
K2PIG	58	-2253.4091	4622.8183	0.0000	1.0000
K2PI	57	-2256.1907	4626.3814	0.0000	1.0000
HKY	59	-2256.0599	4630.1198	0.0000	1.0000
SYM	60	-2278.0928	4676.1857	0.0000	1.0000
K2P	56	-2287.2980	4686.5959	0.0000	1.0000
F81G	59	-2294.0380	4706.0759	0.0000	1.0000
F81IG	60	-2293.8054	4707.6107	0.0000	1.0000
F81I	59	-2299.1930	4716.3859	0.0000	1.0000
JC69G	56	-2315.1408	4742.2816	0.0000	1.0000
JC69IG	57	-2318.7441	4751.4883	0.0000	1.0000
JC69I	56	-2320.5723	4753.1446	0.0000	1.0000
F81	58	-2328.8977	4773.7955	0.0000	1.0000
JC69	55	-2350.8482	4811.6965	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2200.2512	4803.0441	0.9426	0.9426
GTRIG	65	-2199.9966	4808.8247	0.0524	0.9950
GTRI	64	-2205.5726	4813.6870	0.0046	0.9996
HKYG	60	-2220.6571	4818.6971	0.0004	1.0000
HKYIG	61	-2220.4298	4824.5322	0.0000	1.0000
HKYI	60	-2225.6286	4828.6402	0.0000	1.0000
K2PG	57	-2250.5460	4859.6058	0.0000	1.0000
SYMG	61	-2241.9249	4867.5224	0.0000	1.0000
K2PI	57	-2256.1907	4870.8951	0.0000	1.0000
K2PIG	58	-2253.4091	4871.6218	0.0000	1.0000
GTR	63	-2237.7368	4871.7257	0.0000	1.0000
SYMI	61	-2247.2160	4878.1047	0.0000	1.0000
SYMIG	62	-2244.6724	4879.3071	0.0000	1.0000
HKY	59	-2256.0599	4883.2131	0.0000	1.0000
K2P	56	-2287.2980	4926.8200	0.0000	1.0000
SYM	60	-2278.0928	4933.5686	0.0000	1.0000
F81G	59	-2294.0380	4959.1691	0.0000	1.0000
F81IG	60	-2293.8054	4964.9937	0.0000	1.0000
F81I	59	-2299.1930	4969.4791	0.0000	1.0000
JC69G	56	-2315.1408	4982.5056	0.0000	1.0000
JC69I	56	-2320.5723	4993.3686	0.0000	1.0000
JC69IG	57	-2318.7441	4996.0021	0.0000	1.0000
F81	58	-2328.8977	5022.5990	0.0000	1.0000
JC69	55	-2350.8482	5047.6308	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/4976_hg18.phylip = 1 - 539;
 Partition Dummy = 1:261-aggressive-phylip/4976_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/4976_hg18.phylip = 1 - 539;
 Partition Dummy = 1:261-aggressive-phylip/4976_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/4976_hg18.phylip = 1 - 539;
 Partition Dummy = 1:261-aggressive-phylip/4976_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

