
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
01:02 02/17/2011

Input data from file "183-phylip/chr7_4941.phylip" (29 taxa, 678 characters)

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


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	-2507.0190	5155.6106	0.8288	0.8288
GTRIG	65	-2507.5625	5159.1446	0.1416	0.9704
HKYG	60	-2515.3903	5162.6445	0.0246	0.9951
HKYIG	61	-2515.8005	5165.8803	0.0049	0.9999
GTRI	64	-2516.5048	5174.5821	0.0001	1.0000
HKYI	60	-2524.2843	5180.4324	0.0000	1.0000
K2PG	57	-2541.9781	5208.6207	0.0000	1.0000
K2PIG	58	-2541.8654	5210.7873	0.0000	1.0000
SYMG	61	-2538.7647	5211.8087	0.0000	1.0000
SYMIG	62	-2538.6655	5214.0335	0.0000	1.0000
K2PI	57	-2550.3555	5225.3754	0.0000	1.0000
SYMI	61	-2546.9588	5228.1968	0.0000	1.0000
GTR	63	-2560.6091	5260.3518	0.0000	1.0000
HKY	59	-2566.2340	5261.9244	0.0000	1.0000
F81G	59	-2576.3751	5282.2065	0.0000	1.0000
F81IG	60	-2576.7446	5285.3530	0.0000	1.0000
F81I	59	-2584.8168	5299.0900	0.0000	1.0000
K2P	56	-2589.2942	5300.8686	0.0000	1.0000
SYM	60	-2586.6143	5305.0924	0.0000	1.0000
JC69G	56	-2595.5583	5313.3967	0.0000	1.0000
JC69IG	57	-2595.5210	5315.7066	0.0000	1.0000
JC69I	56	-2603.7912	5329.8626	0.0000	1.0000
F81	58	-2622.5265	5372.1095	0.0000	1.0000
JC69	55	-2641.1057	5402.1149	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2507.0190	5142.0380	0.8134	0.8134
GTRIG	65	-2507.5625	5145.1250	0.1738	0.9872
HKYG	60	-2515.3903	5150.7807	0.0103	0.9974
HKYIG	61	-2515.8005	5153.6010	0.0025	0.9999
GTRI	64	-2516.5048	5161.0095	0.0001	1.0000
HKYI	60	-2524.2843	5168.5686	0.0000	1.0000
K2PG	57	-2541.9781	5197.9562	0.0000	1.0000
SYMG	61	-2538.7647	5199.5295	0.0000	1.0000
K2PIG	58	-2541.8654	5199.7307	0.0000	1.0000
SYMIG	62	-2538.6655	5201.3310	0.0000	1.0000
K2PI	57	-2550.3555	5214.7109	0.0000	1.0000
SYMI	61	-2546.9588	5215.9176	0.0000	1.0000
GTR	63	-2560.6091	5247.2183	0.0000	1.0000
HKY	59	-2566.2340	5250.4681	0.0000	1.0000
F81G	59	-2576.3751	5270.7502	0.0000	1.0000
F81IG	60	-2576.7446	5273.4892	0.0000	1.0000
F81I	59	-2584.8168	5287.6337	0.0000	1.0000
K2P	56	-2589.2942	5290.5884	0.0000	1.0000
SYM	60	-2586.6143	5293.2286	0.0000	1.0000
JC69G	56	-2595.5583	5303.1165	0.0000	1.0000
JC69IG	57	-2595.5210	5305.0421	0.0000	1.0000
JC69I	56	-2603.7912	5319.5824	0.0000	1.0000
F81	58	-2622.5265	5361.0530	0.0000	1.0000
JC69	55	-2641.1057	5392.2114	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2515.3903	5421.9295	0.9660	0.9660
HKYIG	61	-2515.8005	5429.2690	0.0246	0.9906
GTRG	64	-2507.0190	5431.2634	0.0091	0.9997
GTRIG	65	-2507.5625	5438.8696	0.0002	0.9999
HKYI	60	-2524.2843	5439.7174	0.0001	1.0000
GTRI	64	-2516.5048	5450.2350	0.0000	1.0000
K2PG	57	-2541.9781	5455.5476	0.0000	1.0000
K2PIG	58	-2541.8654	5461.8413	0.0000	1.0000
K2PI	57	-2550.3555	5472.3023	0.0000	1.0000
SYMG	61	-2538.7647	5475.1974	0.0000	1.0000
SYMIG	62	-2538.6655	5481.5182	0.0000	1.0000
SYMI	61	-2546.9588	5491.5856	0.0000	1.0000
HKY	59	-2566.2340	5517.0978	0.0000	1.0000
GTR	63	-2560.6091	5531.9245	0.0000	1.0000
F81G	59	-2576.3751	5537.3799	0.0000	1.0000
K2P	56	-2589.2942	5543.6606	0.0000	1.0000
F81IG	60	-2576.7446	5544.6380	0.0000	1.0000
F81I	59	-2584.8168	5554.2633	0.0000	1.0000
JC69G	56	-2595.5583	5556.1888	0.0000	1.0000
JC69IG	57	-2595.5210	5562.6335	0.0000	1.0000
SYM	60	-2586.6143	5564.3774	0.0000	1.0000
JC69I	56	-2603.7912	5572.6546	0.0000	1.0000
F81	58	-2622.5265	5623.1635	0.0000	1.0000
JC69	55	-2641.1057	5640.7645	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 183-phylip/chr7_4941.phylip = 1 - 678;
 Partition Dummy = 1:183-phylip/chr7_4941.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 183-phylip/chr7_4941.phylip = 1 - 678;
 Partition Dummy = 1:183-phylip/chr7_4941.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr7_4941.phylip = 1 - 678;
 Partition Dummy = 1:183-phylip/chr7_4941.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

