
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:17 04/04/2011

Input data from file "683-tomodel-phylip/chr6_3276.phylip" (7 taxa, 722 characters)

Minimum AIC  model: GTRI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-2444.4009	4921.5734	0.3461	0.3461
GTRI	20	-2440.3744	4921.9471	0.2871	0.6332
HKYIG	17	-2444.4105	4923.6904	0.1201	0.7532
GTRIG	21	-2440.3859	4924.0919	0.0982	0.8515
HKYG	16	-2445.7271	4924.2258	0.0919	0.9434
GTRG	20	-2441.9977	4925.1937	0.0566	1.0000
K2PI	13	-2459.4172	4945.3484	0.0000	1.0000
SYMI	17	-2455.8881	4946.6454	0.0000	1.0000
K2PIG	14	-2459.4264	4947.4468	0.0000	1.0000
K2PG	13	-2460.6554	4947.8250	0.0000	1.0000
SYMIG	18	-2455.9108	4948.7946	0.0000	1.0000
SYMG	17	-2457.7908	4950.4509	0.0000	1.0000
GTR	19	-2463.7539	4966.5904	0.0000	1.0000
HKY	15	-2469.0443	4968.7684	0.0000	1.0000
K2P	12	-2487.0929	4998.6258	0.0000	1.0000
SYM	16	-2483.8066	5000.3849	0.0000	1.0000
F81I	15	-2512.9968	5056.6735	0.0000	1.0000
F81IG	16	-2513.0208	5058.8133	0.0000	1.0000
F81G	15	-2514.6017	5059.8833	0.0000	1.0000
JC69I	12	-2526.3910	5077.2220	0.0000	1.0000
JC69IG	13	-2526.4164	5079.3470	0.0000	1.0000
JC69G	12	-2528.0976	5080.6352	0.0000	1.0000
F81	14	-2535.1233	5098.8407	0.0000	1.0000
JC69	11	-2550.0082	5122.3882	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2440.3744	4920.7488	0.3176	0.3176
HKYI	16	-2444.4009	4920.8018	0.3093	0.6270
GTRIG	21	-2440.3859	4922.7719	0.1155	0.7425
HKYIG	17	-2444.4105	4922.8211	0.1127	0.8552
HKYG	16	-2445.7271	4923.4542	0.0821	0.9373
GTRG	20	-2441.9977	4923.9955	0.0627	1.0000
K2PI	13	-2459.4172	4944.8343	0.0000	1.0000
SYMI	17	-2455.8881	4945.7761	0.0000	1.0000
K2PIG	14	-2459.4264	4946.8527	0.0000	1.0000
K2PG	13	-2460.6554	4947.3109	0.0000	1.0000
SYMIG	18	-2455.9108	4947.8216	0.0000	1.0000
SYMG	17	-2457.7908	4949.5816	0.0000	1.0000
GTR	19	-2463.7539	4965.5078	0.0000	1.0000
HKY	15	-2469.0443	4968.0886	0.0000	1.0000
K2P	12	-2487.0929	4998.1858	0.0000	1.0000
SYM	16	-2483.8066	4999.6133	0.0000	1.0000
F81I	15	-2512.9968	5055.9936	0.0000	1.0000
F81IG	16	-2513.0208	5058.0417	0.0000	1.0000
F81G	15	-2514.6017	5059.2034	0.0000	1.0000
JC69I	12	-2526.3910	5076.7819	0.0000	1.0000
JC69IG	13	-2526.4164	5078.8328	0.0000	1.0000
JC69G	12	-2528.0976	5080.1952	0.0000	1.0000
F81	14	-2535.1233	5098.2466	0.0000	1.0000
JC69	11	-2550.0082	5122.0164	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2444.4009	4994.1142	0.7632	0.7632
HKYG	16	-2445.7271	4996.7666	0.2026	0.9659
HKYIG	17	-2444.4105	5000.7155	0.0281	0.9940
K2PI	13	-2459.4172	5004.4006	0.0045	0.9984
K2PG	13	-2460.6554	5006.8772	0.0013	0.9997
K2PIG	14	-2459.4264	5011.0011	0.0002	0.9999
GTRI	20	-2440.3744	5012.3893	0.0001	1.0000
GTRG	20	-2441.9977	5015.6360	0.0000	1.0000
GTRIG	21	-2440.3859	5018.9944	0.0000	1.0000
SYMI	17	-2455.8881	5023.6705	0.0000	1.0000
SYMG	17	-2457.7908	5027.4760	0.0000	1.0000
SYMIG	18	-2455.9108	5030.2981	0.0000	1.0000
HKY	15	-2469.0443	5036.8189	0.0000	1.0000
GTR	19	-2463.7539	5052.5663	0.0000	1.0000
K2P	12	-2487.0929	5053.1701	0.0000	1.0000
SYM	16	-2483.8066	5072.9257	0.0000	1.0000
F81I	15	-2512.9968	5124.7240	0.0000	1.0000
F81G	15	-2514.6017	5127.9338	0.0000	1.0000
F81IG	16	-2513.0208	5131.3541	0.0000	1.0000
JC69I	12	-2526.3910	5131.7662	0.0000	1.0000
JC69G	12	-2528.0976	5135.1795	0.0000	1.0000
JC69IG	13	-2526.4164	5138.3992	0.0000	1.0000
F81	14	-2535.1233	5162.3950	0.0000	1.0000
JC69	11	-2550.0082	5172.4186	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_3276.phylip = 1 - 722;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3276.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_3276.phylip = 1 - 722;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3276.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_3276.phylip = 1 - 722;
 Partition Dummy = 1:683-tomodel-phylip/chr6_3276.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

