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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-2471.9120	4985.3596	0.6094	0.6094
GTRIG	21	-2471.9148	4987.5220	0.2067	0.8162
HKYG	16	-2478.4355	4989.8583	0.0643	0.8804
SYMG	17	-2477.6276	4990.3679	0.0498	0.9302
HKYIG	17	-2478.4357	4991.9841	0.0222	0.9525
GTRI	20	-2475.2681	4992.0718	0.0213	0.9737
SYMIG	18	-2477.6102	4992.4664	0.0174	0.9911
HKYI	16	-2481.4835	4995.9544	0.0030	0.9942
SYMI	17	-2480.4669	4996.0465	0.0029	0.9971
K2PG	13	-2485.0476	4996.7522	0.0020	0.9992
K2PIG	14	-2485.0471	4998.8538	0.0007	0.9999
K2PI	13	-2487.8386	5002.3342	0.0001	1.0000
GTR	19	-2505.1011	5049.5891	0.0000	1.0000
HKY	15	-2511.5066	5053.8828	0.0000	1.0000
SYM	16	-2511.3496	5055.6866	0.0000	1.0000
K2P	12	-2518.2593	5061.0807	0.0000	1.0000
F81G	15	-2536.5487	5103.9669	0.0000	1.0000
F81IG	16	-2536.5470	5106.0813	0.0000	1.0000
F81I	15	-2539.4640	5109.7976	0.0000	1.0000
JC69G	12	-2543.2968	5111.1558	0.0000	1.0000
JC69IG	13	-2543.2964	5113.2497	0.0000	1.0000
JC69I	12	-2546.4181	5117.3983	0.0000	1.0000
F81	14	-2567.3663	5163.4920	0.0000	1.0000
JC69	11	-2574.5294	5171.5336	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2471.9120	4983.8239	0.6199	0.6199
GTRIG	21	-2471.9148	4985.8297	0.2274	0.8474
HKYG	16	-2478.4355	4988.8710	0.0497	0.8971
SYMG	17	-2477.6276	4989.2552	0.0410	0.9381
GTRI	20	-2475.2681	4990.5362	0.0216	0.9597
HKYIG	17	-2478.4357	4990.8714	0.0183	0.9780
SYMIG	18	-2477.6102	4991.2205	0.0154	0.9933
SYMI	17	-2480.4669	4994.9338	0.0024	0.9957
HKYI	16	-2481.4835	4994.9671	0.0024	0.9981
K2PG	13	-2485.0476	4996.0951	0.0013	0.9994
K2PIG	14	-2485.0471	4998.0943	0.0005	0.9999
K2PI	13	-2487.8386	5001.6772	0.0001	1.0000
GTR	19	-2505.1011	5048.2023	0.0000	1.0000
HKY	15	-2511.5066	5053.0133	0.0000	1.0000
SYM	16	-2511.3496	5054.6993	0.0000	1.0000
K2P	12	-2518.2593	5060.5185	0.0000	1.0000
F81G	15	-2536.5487	5103.0973	0.0000	1.0000
F81IG	16	-2536.5470	5105.0940	0.0000	1.0000
F81I	15	-2539.4640	5108.9280	0.0000	1.0000
JC69G	12	-2543.2968	5110.5936	0.0000	1.0000
JC69IG	13	-2543.2964	5112.5927	0.0000	1.0000
JC69I	12	-2546.4181	5116.8362	0.0000	1.0000
F81	14	-2567.3663	5162.7325	0.0000	1.0000
JC69	11	-2574.5294	5171.0588	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2485.0476	5052.5427	0.8554	0.8554
K2PI	13	-2487.8386	5058.1248	0.0525	0.9079
HKYG	16	-2478.4355	5058.3450	0.0470	0.9549
K2PIG	14	-2485.0471	5058.8840	0.0359	0.9908
SYMG	17	-2477.6276	5063.0712	0.0044	0.9952
HKYI	16	-2481.4835	5064.4410	0.0022	0.9975
HKYIG	17	-2478.4357	5064.6874	0.0020	0.9994
SYMI	17	-2480.4669	5068.7498	0.0003	0.9997
SYMIG	18	-2477.6102	5069.3787	0.0002	0.9999
GTRG	20	-2471.9120	5070.6663	0.0001	1.0000
GTRIG	21	-2471.9148	5077.0142	0.0000	1.0000
GTRI	20	-2475.2681	5077.3786	0.0000	1.0000
K2P	12	-2518.2593	5112.6240	0.0000	1.0000
HKY	15	-2511.5066	5118.1451	0.0000	1.0000
SYM	16	-2511.3496	5124.1732	0.0000	1.0000
GTR	19	-2505.1011	5130.7026	0.0000	1.0000
JC69G	12	-2543.2968	5162.6991	0.0000	1.0000
F81G	15	-2536.5487	5168.2291	0.0000	1.0000
JC69I	12	-2546.4181	5168.9416	0.0000	1.0000
JC69IG	13	-2543.2964	5169.0403	0.0000	1.0000
F81I	15	-2539.4640	5174.0598	0.0000	1.0000
F81IG	16	-2536.5470	5174.5680	0.0000	1.0000
JC69	11	-2574.5294	5218.8221	0.0000	1.0000
F81	14	-2567.3663	5223.5222	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 591-tomodel-phylip/chr11_3545.phylip = 1 - 568;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3545.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 591-tomodel-phylip/chr11_3545.phylip = 1 - 568;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3545.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 (K2PG)]

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

