
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:30 04/04/2011

Input data from file "591-tomodel-phylip/chr11_3419.phylip" (7 taxa, 667 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	21	-2342.9179	4729.2684	0.9915	0.9915
HKYG	16	-2353.7344	4740.3058	0.0040	0.9955
HKYIG	17	-2353.5356	4742.0141	0.0017	0.9972
GTRI	20	-2350.3716	4742.0434	0.0017	0.9989
GTRG	20	-2350.8546	4743.0094	0.0010	0.9999
HKYI	16	-2357.9325	4748.7019	0.0001	1.0000
K2PG	13	-2362.7387	4752.0348	0.0000	1.0000
SYMG	17	-2358.6273	4752.1976	0.0000	1.0000
SYMI	17	-2358.6447	4752.2325	0.0000	1.0000
K2PI	13	-2363.0574	4752.6723	0.0000	1.0000
SYMIG	18	-2358.3194	4753.6945	0.0000	1.0000
K2PIG	14	-2362.6277	4753.8997	0.0000	1.0000
GTR	19	-2393.8253	4826.8253	0.0000	1.0000
HKY	15	-2402.1764	4835.0901	0.0000	1.0000
SYM	16	-2403.8947	4840.6263	0.0000	1.0000
K2P	12	-2409.1689	4842.8148	0.0000	1.0000
F81I	15	-2420.4641	4871.6654	0.0000	1.0000
F81IG	16	-2420.4808	4873.7986	0.0000	1.0000
F81G	15	-2422.1617	4875.0608	0.0000	1.0000
JC69IG	13	-2424.7092	4875.9758	0.0000	1.0000
JC69I	12	-2425.8142	4876.1055	0.0000	1.0000
JC69G	12	-2425.8954	4876.2679	0.0000	1.0000
F81	14	-2461.6750	4951.9941	0.0000	1.0000
JC69	11	-2467.9568	4958.3167	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-2342.9179	4727.8358	0.9931	0.9931
HKYG	16	-2353.7344	4739.4689	0.0030	0.9961
GTRI	20	-2350.3716	4740.7431	0.0016	0.9976
HKYIG	17	-2353.5356	4741.0711	0.0013	0.9990
GTRG	20	-2350.8546	4741.7091	0.0010	0.9999
HKYI	16	-2357.9325	4747.8650	0.0000	1.0000
SYMG	17	-2358.6273	4751.2546	0.0000	1.0000
SYMI	17	-2358.6447	4751.2895	0.0000	1.0000
K2PG	13	-2362.7387	4751.4774	0.0000	1.0000
K2PI	13	-2363.0574	4752.1149	0.0000	1.0000
SYMIG	18	-2358.3194	4752.6389	0.0000	1.0000
K2PIG	14	-2362.6277	4753.2555	0.0000	1.0000
GTR	19	-2393.8253	4825.6507	0.0000	1.0000
HKY	15	-2402.1764	4834.3528	0.0000	1.0000
SYM	16	-2403.8947	4839.7894	0.0000	1.0000
K2P	12	-2409.1689	4842.3377	0.0000	1.0000
F81I	15	-2420.4641	4870.9281	0.0000	1.0000
F81IG	16	-2420.4808	4872.9617	0.0000	1.0000
F81G	15	-2422.1617	4874.3235	0.0000	1.0000
JC69IG	13	-2424.7092	4875.4184	0.0000	1.0000
JC69I	12	-2425.8142	4875.6285	0.0000	1.0000
JC69G	12	-2425.8954	4875.7909	0.0000	1.0000
F81	14	-2461.6750	4951.3499	0.0000	1.0000
JC69	11	-2467.9568	4957.9136	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2362.7387	4810.0137	0.4397	0.4397
K2PI	13	-2363.0574	4810.6512	0.3197	0.7593
HKYG	16	-2353.7344	4811.5135	0.2077	0.9670
K2PIG	14	-2362.6277	4816.2945	0.0190	0.9860
HKYIG	17	-2353.5356	4817.6185	0.0098	0.9958
HKYI	16	-2357.9325	4819.9097	0.0031	0.9990
GTRIG	21	-2342.9179	4822.3944	0.0009	0.9999
SYMG	17	-2358.6273	4827.8021	0.0001	0.9999
SYMI	17	-2358.6447	4827.8369	0.0001	1.0000
GTRI	20	-2350.3716	4830.7989	0.0000	1.0000
GTRG	20	-2350.8546	4831.7649	0.0000	1.0000
SYMIG	18	-2358.3194	4833.6891	0.0000	1.0000
K2P	12	-2409.1689	4896.3712	0.0000	1.0000
HKY	15	-2402.1764	4901.8946	0.0000	1.0000
GTR	19	-2393.8253	4911.2037	0.0000	1.0000
SYM	16	-2403.8947	4911.8340	0.0000	1.0000
JC69I	12	-2425.8142	4929.6620	0.0000	1.0000
JC69G	12	-2425.8954	4929.8243	0.0000	1.0000
JC69IG	13	-2424.7092	4933.9547	0.0000	1.0000
F81I	15	-2420.4641	4938.4700	0.0000	1.0000
F81G	15	-2422.1617	4941.8653	0.0000	1.0000
F81IG	16	-2420.4808	4945.0063	0.0000	1.0000
JC69	11	-2467.9568	5007.4443	0.0000	1.0000
F81	14	-2461.6750	5014.3890	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 591-tomodel-phylip/chr11_3419.phylip = 1 - 667;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3419.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 591-tomodel-phylip/chr11_3419.phylip = 1 - 667;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3419.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3419.phylip = 1 - 667;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3419.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

