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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1903.1297	3847.6186	0.4492	0.4492
GTRG	20	-1903.7233	3848.8058	0.2481	0.6973
GTRIG	21	-1903.1354	3849.7683	0.1533	0.8506
HKYI	16	-1909.1448	3851.1642	0.0763	0.9269
HKYG	16	-1909.6359	3852.1464	0.0467	0.9736
HKYIG	17	-1909.1490	3853.2834	0.0264	1.0000
GTR	19	-1919.5739	3878.3757	0.0000	1.0000
HKY	15	-1924.7588	3880.2881	0.0000	1.0000
F81I	15	-1939.0277	3908.8258	0.0000	1.0000
F81G	15	-1939.3040	3909.3785	0.0000	1.0000
F81IG	16	-1939.0232	3910.9209	0.0000	1.0000
K2PI	13	-1942.4046	3911.3916	0.0000	1.0000
K2PG	13	-1942.7237	3912.0297	0.0000	1.0000
K2PIG	14	-1942.4073	3913.4877	0.0000	1.0000
SYMI	17	-1941.5780	3918.1416	0.0000	1.0000
SYMG	17	-1941.9225	3918.8306	0.0000	1.0000
SYMIG	18	-1941.5811	3920.2654	0.0000	1.0000
F81	14	-1953.8471	3936.3672	0.0000	1.0000
K2P	12	-1957.7370	3939.9725	0.0000	1.0000
SYM	16	-1956.7289	3946.3324	0.0000	1.0000
JC69I	12	-1967.1212	3958.7408	0.0000	1.0000
JC69G	12	-1967.5435	3959.5855	0.0000	1.0000
JC69IG	13	-1967.1277	3960.8379	0.0000	1.0000
JC69	11	-1982.0445	3986.5101	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1903.1297	3846.2593	0.4583	0.4583
GTRG	20	-1903.7233	3847.4466	0.2531	0.7115
GTRIG	21	-1903.1354	3848.2707	0.1677	0.8791
HKYI	16	-1909.1448	3850.2896	0.0611	0.9402
HKYG	16	-1909.6359	3851.2718	0.0374	0.9776
HKYIG	17	-1909.1490	3852.2979	0.0224	1.0000
GTR	19	-1919.5739	3877.1479	0.0000	1.0000
HKY	15	-1924.7588	3879.5177	0.0000	1.0000
F81I	15	-1939.0277	3908.0553	0.0000	1.0000
F81G	15	-1939.3040	3908.6080	0.0000	1.0000
F81IG	16	-1939.0232	3910.0463	0.0000	1.0000
K2PI	13	-1942.4046	3910.8092	0.0000	1.0000
K2PG	13	-1942.7237	3911.4473	0.0000	1.0000
K2PIG	14	-1942.4073	3912.8146	0.0000	1.0000
SYMI	17	-1941.5780	3917.1561	0.0000	1.0000
SYMG	17	-1941.9225	3917.8451	0.0000	1.0000
SYMIG	18	-1941.5811	3919.1622	0.0000	1.0000
F81	14	-1953.8471	3935.6941	0.0000	1.0000
K2P	12	-1957.7370	3939.4741	0.0000	1.0000
SYM	16	-1956.7289	3945.4578	0.0000	1.0000
JC69I	12	-1967.1212	3958.2424	0.0000	1.0000
JC69G	12	-1967.5435	3959.0871	0.0000	1.0000
JC69IG	13	-1967.1277	3960.2555	0.0000	1.0000
JC69	11	-1982.0445	3986.0891	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1909.1448	3921.6481	0.6050	0.6050
HKYG	16	-1909.6359	3922.6303	0.3702	0.9752
HKYIG	17	-1909.1490	3928.1163	0.0238	0.9990
GTRI	20	-1903.1297	3935.4574	0.0006	0.9996
GTRG	20	-1903.7233	3936.6447	0.0003	1.0000
GTRIG	21	-1903.1354	3941.9287	0.0000	1.0000
HKY	15	-1924.7588	3946.4162	0.0000	1.0000
GTR	19	-1919.5739	3961.8861	0.0000	1.0000
K2PI	13	-1942.4046	3968.7880	0.0000	1.0000
K2PG	13	-1942.7237	3969.4261	0.0000	1.0000
F81I	15	-1939.0277	3974.9539	0.0000	1.0000
K2PIG	14	-1942.4073	3975.2533	0.0000	1.0000
F81G	15	-1939.3040	3975.5066	0.0000	1.0000
F81IG	16	-1939.0232	3981.4048	0.0000	1.0000
SYMI	17	-1941.5780	3992.9744	0.0000	1.0000
K2P	12	-1957.7370	3992.9929	0.0000	1.0000
SYMG	17	-1941.9225	3993.6635	0.0000	1.0000
F81	14	-1953.8471	3998.1328	0.0000	1.0000
SYMIG	18	-1941.5811	3999.4405	0.0000	1.0000
JC69I	12	-1967.1212	4011.7613	0.0000	1.0000
JC69G	12	-1967.5435	4012.6059	0.0000	1.0000
SYM	16	-1956.7289	4016.8162	0.0000	1.0000
JC69IG	13	-1967.1277	4018.2342	0.0000	1.0000
JC69	11	-1982.0445	4035.1480	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 591-tomodel-phylip/chr3_18256.phylip = 1 - 639;
 Partition Dummy = 1:591-tomodel-phylip/chr3_18256.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 (GTRI)]

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

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

