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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-3775.7475	7694.6143	0.4984	0.4984
HKYIG	61	-3781.3298	7696.1551	0.2307	0.7291
GTRG	64	-3777.8825	7696.4674	0.1973	0.9264
HKYG	60	-3783.6661	7698.4399	0.0736	1.0000
HKYI	60	-3804.4024	7739.9126	0.0000	1.0000
GTRI	64	-3800.3643	7741.4309	0.0000	1.0000
K2PIG	58	-3808.0505	7742.4551	0.0000	1.0000
K2PG	57	-3810.0297	7744.0472	0.0000	1.0000
SYMIG	62	-3804.5127	7744.9157	0.0000	1.0000
SYMG	61	-3807.2449	7747.9852	0.0000	1.0000
K2PI	57	-3831.0927	7786.1734	0.0000	1.0000
SYMI	61	-3827.9772	7789.4499	0.0000	1.0000
F81IG	60	-3870.2646	7871.6369	0.0000	1.0000
F81G	59	-3873.3287	7875.3846	0.0000	1.0000
F81I	59	-3892.2264	7913.1801	0.0000	1.0000
JC69IG	57	-3895.6997	7915.3873	0.0000	1.0000
JC69G	56	-3898.3288	7918.2865	0.0000	1.0000
JC69I	56	-3917.5180	7956.6650	0.0000	1.0000
GTR	63	-3917.8347	7973.9622	0.0000	1.0000
HKY	59	-3922.9904	7974.7080	0.0000	1.0000
K2P	56	-3948.7073	8019.0435	0.0000	1.0000
SYM	60	-3947.0271	8025.1618	0.0000	1.0000
F81	58	-4008.2813	8142.9166	0.0000	1.0000
JC69	55	-4033.8368	8186.9506	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3775.7475	7681.4951	0.6325	0.6325
GTRG	64	-3777.8825	7683.7651	0.2033	0.8358
HKYIG	61	-3781.3298	7684.6597	0.1300	0.9658
HKYG	60	-3783.6661	7687.3322	0.0342	1.0000
GTRI	64	-3800.3643	7728.7286	0.0000	1.0000
HKYI	60	-3804.4024	7728.8049	0.0000	1.0000
K2PIG	58	-3808.0505	7732.1011	0.0000	1.0000
SYMIG	62	-3804.5127	7733.0253	0.0000	1.0000
K2PG	57	-3810.0297	7734.0593	0.0000	1.0000
SYMG	61	-3807.2449	7736.4897	0.0000	1.0000
K2PI	57	-3831.0927	7776.1855	0.0000	1.0000
SYMI	61	-3827.9772	7777.9545	0.0000	1.0000
F81IG	60	-3870.2646	7860.5292	0.0000	1.0000
F81G	59	-3873.3287	7864.6574	0.0000	1.0000
F81I	59	-3892.2264	7902.4529	0.0000	1.0000
JC69IG	57	-3895.6997	7905.3994	0.0000	1.0000
JC69G	56	-3898.3288	7908.6575	0.0000	1.0000
JC69I	56	-3917.5180	7947.0360	0.0000	1.0000
GTR	63	-3917.8347	7961.6695	0.0000	1.0000
HKY	59	-3922.9904	7963.9808	0.0000	1.0000
K2P	56	-3948.7073	8009.4145	0.0000	1.0000
SYM	60	-3947.0271	8014.0541	0.0000	1.0000
F81	58	-4008.2813	8132.5626	0.0000	1.0000
JC69	55	-4033.8368	8177.6735	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3783.6661	7962.0872	0.7214	0.7214
HKYIG	61	-3781.3298	7963.9940	0.2780	0.9994
GTRG	64	-3777.8825	7976.8372	0.0005	0.9999
GTRIG	65	-3775.7475	7979.1464	0.0001	1.0000
K2PG	57	-3810.0297	7995.0766	0.0000	1.0000
K2PIG	58	-3808.0505	7997.6977	0.0000	1.0000
HKYI	60	-3804.4024	8003.5600	0.0000	1.0000
SYMG	61	-3807.2449	8015.8241	0.0000	1.0000
SYMIG	62	-3804.5127	8016.9389	0.0000	1.0000
GTRI	64	-3800.3643	8021.8007	0.0000	1.0000
K2PI	57	-3831.0927	8037.2028	0.0000	1.0000
SYMI	61	-3827.9772	8057.2888	0.0000	1.0000
F81G	59	-3873.3287	8134.8332	0.0000	1.0000
F81IG	60	-3870.2646	8135.2843	0.0000	1.0000
JC69G	56	-3898.3288	8165.0956	0.0000	1.0000
JC69IG	57	-3895.6997	8166.4167	0.0000	1.0000
F81I	59	-3892.2264	8172.6287	0.0000	1.0000
JC69I	56	-3917.5180	8203.4741	0.0000	1.0000
HKY	59	-3922.9904	8234.1566	0.0000	1.0000
GTR	63	-3917.8347	8250.1623	0.0000	1.0000
K2P	56	-3948.7073	8265.8526	0.0000	1.0000
SYM	60	-3947.0271	8288.8092	0.0000	1.0000
F81	58	-4008.2813	8398.1592	0.0000	1.0000
JC69	55	-4033.8368	8429.5324	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/12691_hg18.phylip = 1 - 720;
 Partition Dummy = 1:261-aggressive-phylip/12691_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/12691_hg18.phylip = 1 - 720;
 Partition Dummy = 1:261-aggressive-phylip/12691_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 (HKYG)]

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

