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

Input data from file "183-phylip/chrZ_6727.phylip" (29 taxa, 632 characters)

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: K2PIG


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
HKYIG	61	-1911.9346	3959.1394	0.9438	0.9438
GTRIG	65	-1910.0635	3965.2861	0.0437	0.9874
K2PIG	58	-1919.9820	3967.9082	0.0118	0.9992
SYMIG	62	-1918.1756	3974.0806	0.0005	0.9997
HKYI	60	-1921.4773	3975.7742	0.0002	1.0000
HKYG	60	-1923.5975	3980.0147	0.0000	1.0000
GTRI	64	-1920.1079	3982.8895	0.0000	1.0000
K2PI	57	-1929.2071	3983.9333	0.0000	1.0000
GTRG	64	-1922.2149	3987.1036	0.0000	1.0000
K2PG	57	-1931.3261	3988.1714	0.0000	1.0000
SYMI	61	-1927.9930	3991.2563	0.0000	1.0000
SYMG	61	-1930.2911	3995.8525	0.0000	1.0000
F81IG	60	-1947.7307	4028.2810	0.0000	1.0000
JC69IG	57	-1954.2837	4034.0865	0.0000	1.0000
F81I	59	-1956.9930	4044.3635	0.0000	1.0000
JC69I	56	-1963.5209	4050.1443	0.0000	1.0000
F81G	59	-1960.3546	4051.0869	0.0000	1.0000
JC69G	56	-1966.8277	4056.7580	0.0000	1.0000
HKY	59	-2019.9521	4170.2818	0.0000	1.0000
K2P	56	-2026.9322	4176.9669	0.0000	1.0000
GTR	63	-2019.0900	4178.3772	0.0000	1.0000
SYM	60	-2026.8279	4186.4754	0.0000	1.0000
F81	58	-2056.0409	4240.0260	0.0000	1.0000
JC69	55	-2062.0871	4244.8687	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1911.9346	3945.8692	0.8878	0.8878
GTRIG	65	-1910.0635	3950.1270	0.1056	0.9934
K2PIG	58	-1919.9820	3955.9640	0.0057	0.9992
SYMIG	62	-1918.1756	3960.3513	0.0006	0.9998
HKYI	60	-1921.4773	3962.9546	0.0002	1.0000
HKYG	60	-1923.5975	3967.1950	0.0000	1.0000
GTRI	64	-1920.1079	3968.2158	0.0000	1.0000
K2PI	57	-1929.2071	3972.4142	0.0000	1.0000
GTRG	64	-1922.2149	3972.4299	0.0000	1.0000
K2PG	57	-1931.3261	3976.6522	0.0000	1.0000
SYMI	61	-1927.9930	3977.9861	0.0000	1.0000
SYMG	61	-1930.2911	3982.5823	0.0000	1.0000
F81IG	60	-1947.7307	4015.4614	0.0000	1.0000
JC69IG	57	-1954.2837	4022.5674	0.0000	1.0000
F81I	59	-1956.9930	4031.9859	0.0000	1.0000
F81G	59	-1960.3546	4038.7093	0.0000	1.0000
JC69I	56	-1963.5209	4039.0417	0.0000	1.0000
JC69G	56	-1966.8277	4045.6554	0.0000	1.0000
HKY	59	-2019.9521	4157.9042	0.0000	1.0000
GTR	63	-2019.0900	4164.1801	0.0000	1.0000
K2P	56	-2026.9322	4165.8643	0.0000	1.0000
SYM	60	-2026.8279	4173.6557	0.0000	1.0000
F81	58	-2056.0409	4228.0818	0.0000	1.0000
JC69	55	-2062.0871	4234.1743	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	58	-1919.9820	4213.9996	0.8334	0.8334
HKYIG	61	-1911.9346	4217.2515	0.1639	0.9973
K2PI	57	-1929.2071	4226.0009	0.0021	0.9994
HKYI	60	-1921.4773	4229.8879	0.0003	0.9997
K2PG	57	-1931.3261	4230.2389	0.0002	0.9999
HKYG	60	-1923.5975	4234.1284	0.0000	1.0000
SYMIG	62	-1918.1756	4236.1824	0.0000	1.0000
GTRIG	65	-1910.0635	4239.3049	0.0000	1.0000
SYMI	61	-1927.9930	4249.3683	0.0000	1.0000
GTRI	64	-1920.1079	4252.9447	0.0000	1.0000
SYMG	61	-1930.2911	4253.9645	0.0000	1.0000
GTRG	64	-1922.2149	4257.1588	0.0000	1.0000
JC69IG	57	-1954.2837	4276.1541	0.0000	1.0000
F81IG	60	-1947.7307	4282.3948	0.0000	1.0000
JC69I	56	-1963.5209	4288.1795	0.0000	1.0000
F81I	59	-1956.9930	4294.4704	0.0000	1.0000
JC69G	56	-1966.8277	4294.7932	0.0000	1.0000
F81G	59	-1960.3546	4301.1937	0.0000	1.0000
K2P	56	-2026.9322	4415.0021	0.0000	1.0000
HKY	59	-2019.9521	4420.3887	0.0000	1.0000
SYM	60	-2026.8279	4440.5891	0.0000	1.0000
GTR	63	-2019.0900	4444.4601	0.0000	1.0000
JC69	55	-2062.0871	4478.8632	0.0000	1.0000
F81	58	-2056.0409	4486.1174	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_6727.phylip = 1 - 632;
 Partition Dummy = 1:183-phylip/chrZ_6727.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_6727.phylip = 1 - 632;
 Partition Dummy = 1:183-phylip/chrZ_6727.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_6727.phylip = 1 - 632;
 Partition Dummy = 1:183-phylip/chrZ_6727.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

