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

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

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	13	-1992.5892	4011.8295	0.3256	0.3256
SYMG	17	-1988.9195	4012.9418	0.1867	0.5123
K2PIG	14	-1992.5957	4013.9441	0.1131	0.6254
K2PI	13	-1993.7448	4014.1407	0.1025	0.7280
SYMI	17	-1989.9733	4015.0493	0.0651	0.7931
SYMIG	18	-1988.9082	4015.0511	0.0650	0.8581
HKYG	16	-1991.0675	4015.1135	0.0630	0.9211
GTRG	20	-1987.8355	4017.1927	0.0223	0.9434
HKYIG	17	-1991.0624	4017.2276	0.0219	0.9653
HKYI	16	-1992.2574	4017.4933	0.0192	0.9845
GTRI	20	-1988.8735	4019.2687	0.0079	0.9924
GTRIG	21	-1987.8326	4019.3421	0.0076	1.0000
JC69G	12	-2030.5502	4085.6576	0.0000	1.0000
JC69I	12	-2030.9773	4086.5118	0.0000	1.0000
F81G	15	-2028.3908	4087.6433	0.0000	1.0000
JC69IG	13	-2030.6265	4087.9042	0.0000	1.0000
F81I	15	-2028.8188	4088.4993	0.0000	1.0000
F81IG	16	-2028.4664	4089.9113	0.0000	1.0000
K2P	12	-2035.2812	4095.1196	0.0000	1.0000
SYM	16	-2031.8702	4096.7188	0.0000	1.0000
HKY	15	-2033.1553	4097.1724	0.0000	1.0000
GTR	19	-2030.1620	4099.6982	0.0000	1.0000
JC69	11	-2070.4592	4163.3890	0.0000	1.0000
F81	14	-2068.1743	4165.1014	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	13	-1992.5892	4011.1783	0.2873	0.2873
SYMG	17	-1988.9195	4011.8391	0.2065	0.4937
K2PIG	14	-1992.5957	4013.1914	0.1050	0.5987
K2PI	13	-1993.7448	4013.4895	0.0905	0.6892
SYMIG	18	-1988.9082	4013.8164	0.0768	0.7660
SYMI	17	-1989.9733	4013.9466	0.0720	0.8380
HKYG	16	-1991.0675	4014.1350	0.0655	0.9035
GTRG	20	-1987.8355	4015.6710	0.0304	0.9339
HKYIG	17	-1991.0624	4016.1249	0.0242	0.9581
HKYI	16	-1992.2574	4016.5148	0.0199	0.9780
GTRIG	21	-1987.8326	4017.6652	0.0112	0.9892
GTRI	20	-1988.8735	4017.7469	0.0108	1.0000
JC69G	12	-2030.5502	4085.1004	0.0000	1.0000
JC69I	12	-2030.9773	4085.9547	0.0000	1.0000
F81G	15	-2028.3908	4086.7815	0.0000	1.0000
JC69IG	13	-2030.6265	4087.2531	0.0000	1.0000
F81I	15	-2028.8188	4087.6376	0.0000	1.0000
F81IG	16	-2028.4664	4088.9329	0.0000	1.0000
K2P	12	-2035.2812	4094.5625	0.0000	1.0000
SYM	16	-2031.8702	4095.7404	0.0000	1.0000
HKY	15	-2033.1553	4096.3106	0.0000	1.0000
GTR	19	-2030.1620	4098.3239	0.0000	1.0000
JC69	11	-2070.4592	4162.9184	0.0000	1.0000
F81	14	-2068.1743	4164.3487	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1992.5892	4067.7398	0.7369	0.7369
K2PI	13	-1993.7448	4070.0511	0.2320	0.9690
K2PIG	14	-1992.5957	4074.1038	0.0306	0.9995
HKYG	16	-1991.0675	4083.7492	0.0002	0.9998
SYMG	17	-1988.9195	4085.8041	0.0001	0.9999
HKYI	16	-1992.2574	4086.1290	0.0001	1.0000
SYMI	17	-1989.9733	4087.9116	0.0000	1.0000
HKYIG	17	-1991.0624	4090.0899	0.0000	1.0000
SYMIG	18	-1988.9082	4092.1324	0.0000	1.0000
GTRG	20	-1987.8355	4102.6887	0.0000	1.0000
GTRI	20	-1988.8735	4104.7646	0.0000	1.0000
GTRIG	21	-1987.8326	4109.0338	0.0000	1.0000
JC69G	12	-2030.5502	4137.3110	0.0000	1.0000
JC69I	12	-2030.9773	4138.1653	0.0000	1.0000
JC69IG	13	-2030.6265	4143.8146	0.0000	1.0000
K2P	12	-2035.2812	4146.7731	0.0000	1.0000
F81G	15	-2028.3908	4152.0448	0.0000	1.0000
F81I	15	-2028.8188	4152.9008	0.0000	1.0000
F81IG	16	-2028.4664	4158.5470	0.0000	1.0000
HKY	15	-2033.1553	4161.5739	0.0000	1.0000
SYM	16	-2031.8702	4165.3546	0.0000	1.0000
GTR	19	-2030.1620	4180.9907	0.0000	1.0000
JC69	11	-2070.4592	4210.7782	0.0000	1.0000
F81	14	-2068.1743	4225.2611	0.0000	1.0000



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

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

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


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

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


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

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

