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

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

Minimum AIC  model: HKYI
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	-1929.0300	3890.9277	0.4126	0.4126
HKYG	16	-1929.4001	3891.6678	0.2850	0.6976
HKYIG	17	-1928.9156	3892.8088	0.1611	0.8587
GTRI	20	-1926.6688	3894.6859	0.0630	0.9218
GTRG	20	-1926.8518	3895.0520	0.0525	0.9742
GTRIG	21	-1926.4948	3896.4751	0.0258	1.0000
K2PI	13	-1950.2719	3927.1215	0.0000	1.0000
SYMI	17	-1946.1362	3927.2500	0.0000	1.0000
SYMG	17	-1946.1525	3927.2827	0.0000	1.0000
K2PG	13	-1950.6375	3927.8528	0.0000	1.0000
SYMIG	18	-1945.8529	3928.8001	0.0000	1.0000
K2PIG	14	-1950.1524	3928.9726	0.0000	1.0000
HKY	15	-1957.4113	3945.5869	0.0000	1.0000
GTR	19	-1955.1847	3949.5873	0.0000	1.0000
F81I	15	-1960.9541	3952.6725	0.0000	1.0000
F81G	15	-1961.6247	3954.0138	0.0000	1.0000
F81IG	16	-1960.9382	3954.7440	0.0000	1.0000
SYM	16	-1974.0279	3980.9233	0.0000	1.0000
K2P	12	-1978.8604	3982.2152	0.0000	1.0000
JC69I	12	-1980.5243	3985.5430	0.0000	1.0000
JC69G	12	-1981.1267	3986.7479	0.0000	1.0000
JC69IG	13	-1980.4979	3987.5736	0.0000	1.0000
F81	14	-1989.4793	4007.6264	0.0000	1.0000
JC69	11	-2008.5733	4039.5644	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1929.0300	3890.0601	0.3931	0.3931
HKYG	16	-1929.4001	3890.8002	0.2715	0.6645
HKYIG	17	-1928.9156	3891.8311	0.1621	0.8267
GTRI	20	-1926.6688	3893.3376	0.0763	0.9030
GTRG	20	-1926.8518	3893.7036	0.0636	0.9666
GTRIG	21	-1926.4948	3894.9896	0.0334	1.0000
SYMI	17	-1946.1362	3926.2723	0.0000	1.0000
SYMG	17	-1946.1525	3926.3050	0.0000	1.0000
K2PI	13	-1950.2719	3926.5437	0.0000	1.0000
K2PG	13	-1950.6375	3927.2750	0.0000	1.0000
SYMIG	18	-1945.8529	3927.7057	0.0000	1.0000
K2PIG	14	-1950.1524	3928.3049	0.0000	1.0000
HKY	15	-1957.4113	3944.8225	0.0000	1.0000
GTR	19	-1955.1847	3948.3693	0.0000	1.0000
F81I	15	-1960.9541	3951.9082	0.0000	1.0000
F81G	15	-1961.6247	3953.2495	0.0000	1.0000
F81IG	16	-1960.9382	3953.8764	0.0000	1.0000
SYM	16	-1974.0279	3980.0557	0.0000	1.0000
K2P	12	-1978.8604	3981.7207	0.0000	1.0000
JC69I	12	-1980.5243	3985.0485	0.0000	1.0000
JC69G	12	-1981.1267	3986.2535	0.0000	1.0000
JC69IG	13	-1980.4979	3986.9958	0.0000	1.0000
F81	14	-1989.4793	4006.9587	0.0000	1.0000
JC69	11	-2008.5733	4039.1466	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1929.0300	3961.5432	0.5764	0.5764
HKYG	16	-1929.4001	3962.2834	0.3981	0.9745
HKYIG	17	-1928.9156	3967.7820	0.0255	1.0000
GTRI	20	-1926.6688	3982.6916	0.0000	1.0000
GTRG	20	-1926.8518	3983.0576	0.0000	1.0000
K2PI	13	-1950.2719	3984.6238	0.0000	1.0000
K2PG	13	-1950.6375	3985.3551	0.0000	1.0000
GTRIG	21	-1926.4948	3988.8113	0.0000	1.0000
K2PIG	14	-1950.1524	3990.8527	0.0000	1.0000
SYMI	17	-1946.1362	4002.2232	0.0000	1.0000
SYMG	17	-1946.1525	4002.2559	0.0000	1.0000
SYMIG	18	-1945.8529	4008.1243	0.0000	1.0000
HKY	15	-1957.4113	4011.8380	0.0000	1.0000
F81I	15	-1960.9541	4018.9236	0.0000	1.0000
F81G	15	-1961.6247	4020.2649	0.0000	1.0000
F81IG	16	-1960.9382	4025.3596	0.0000	1.0000
GTR	19	-1955.1847	4033.2556	0.0000	1.0000
K2P	12	-1978.8604	4035.3331	0.0000	1.0000
JC69I	12	-1980.5243	4038.6609	0.0000	1.0000
JC69G	12	-1981.1267	4039.8659	0.0000	1.0000
JC69IG	13	-1980.4979	4045.0759	0.0000	1.0000
SYM	16	-1974.0279	4051.5389	0.0000	1.0000
F81	14	-1989.4793	4069.5064	0.0000	1.0000
JC69	11	-2008.5733	4088.2913	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_5476.phylip = 1 - 644;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5476.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 AICc model (HKYI)]

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

