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

Input data from file "683-tomodel-phylip/chr5_11334.phylip" (7 taxa, 593 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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	-1901.0045	3843.4776	0.5242	0.5242
GTRG	20	-1901.5726	3844.6136	0.2970	0.8213
GTRIG	21	-1901.0074	3845.6330	0.1784	0.9997
HKYI	16	-1913.8044	3860.5532	0.0001	0.9998
GTR	19	-1910.6725	3860.6713	0.0001	0.9999
HKYG	16	-1914.3726	3861.6896	0.0001	1.0000
HKYIG	17	-1913.8126	3862.6896	0.0000	1.0000
HKY	15	-1923.6471	3878.1261	0.0000	1.0000
SYMI	17	-1937.9374	3910.9391	0.0000	1.0000
K2PI	13	-1942.5460	3911.7208	0.0000	1.0000
SYMG	17	-1938.6494	3912.3632	0.0000	1.0000
K2PG	13	-1943.0774	3912.7836	0.0000	1.0000
SYMIG	18	-1937.9511	3913.0939	0.0000	1.0000
K2PIG	14	-1942.5551	3913.8369	0.0000	1.0000
SYM	16	-1949.2903	3931.5250	0.0000	1.0000
K2P	12	-1954.3404	3933.2187	0.0000	1.0000
F81I	15	-1969.4386	3969.7092	0.0000	1.0000
F81G	15	-1970.2278	3971.2875	0.0000	1.0000
F81IG	16	-1969.4543	3971.8531	0.0000	1.0000
F81	14	-1979.2400	3987.2066	0.0000	1.0000
JC69I	12	-1993.1667	4010.8713	0.0000	1.0000
JC69G	12	-1993.8406	4012.2191	0.0000	1.0000
JC69IG	13	-1993.1796	4012.9880	0.0000	1.0000
JC69	11	-2003.8695	4030.1934	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1901.0045	3842.0091	0.5171	0.5171
GTRG	20	-1901.5726	3843.1451	0.2930	0.8101
GTRIG	21	-1901.0074	3844.0147	0.1897	0.9998
GTR	19	-1910.6725	3859.3449	0.0001	0.9998
HKYI	16	-1913.8044	3859.6088	0.0001	0.9999
HKYG	16	-1914.3726	3860.7452	0.0000	1.0000
HKYIG	17	-1913.8126	3861.6253	0.0000	1.0000
HKY	15	-1923.6471	3877.2942	0.0000	1.0000
SYMI	17	-1937.9374	3909.8748	0.0000	1.0000
K2PI	13	-1942.5460	3911.0921	0.0000	1.0000
SYMG	17	-1938.6494	3911.2989	0.0000	1.0000
SYMIG	18	-1937.9511	3911.9023	0.0000	1.0000
K2PG	13	-1943.0774	3912.1549	0.0000	1.0000
K2PIG	14	-1942.5551	3913.1102	0.0000	1.0000
SYM	16	-1949.2903	3930.5806	0.0000	1.0000
K2P	12	-1954.3404	3932.6808	0.0000	1.0000
F81I	15	-1969.4386	3968.8773	0.0000	1.0000
F81G	15	-1970.2278	3970.4557	0.0000	1.0000
F81IG	16	-1969.4543	3970.9086	0.0000	1.0000
F81	14	-1979.2400	3986.4799	0.0000	1.0000
JC69I	12	-1993.1667	4010.3333	0.0000	1.0000
JC69G	12	-1993.8406	4011.6812	0.0000	1.0000
JC69IG	13	-1993.1796	4012.3593	0.0000	1.0000
JC69	11	-2003.8695	4029.7390	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	20	-1901.0045	3929.7130	0.3154	0.3154
HKYI	16	-1913.8044	3929.7719	0.3062	0.6216
GTRG	20	-1901.5726	3930.8490	0.1787	0.8003
HKYG	16	-1914.3726	3930.9083	0.1735	0.9737
GTRIG	21	-1901.0074	3936.1038	0.0129	0.9866
HKYIG	17	-1913.8126	3936.1736	0.0125	0.9991
GTR	19	-1910.6725	3942.6636	0.0005	0.9996
HKY	15	-1923.6471	3943.0721	0.0004	1.0000
K2PI	13	-1942.5460	3968.0996	0.0000	1.0000
K2PG	13	-1943.0774	3969.1624	0.0000	1.0000
K2PIG	14	-1942.5551	3974.5029	0.0000	1.0000
SYMI	17	-1937.9374	3984.4231	0.0000	1.0000
K2P	12	-1954.3404	3985.3031	0.0000	1.0000
SYMG	17	-1938.6494	3985.8472	0.0000	1.0000
SYMIG	18	-1937.9511	3990.8358	0.0000	1.0000
SYM	16	-1949.2903	4000.7437	0.0000	1.0000
F81I	15	-1969.4386	4034.6552	0.0000	1.0000
F81G	15	-1970.2278	4036.2336	0.0000	1.0000
F81IG	16	-1969.4543	4041.0718	0.0000	1.0000
F81	14	-1979.2400	4047.8727	0.0000	1.0000
JC69I	12	-1993.1667	4062.9557	0.0000	1.0000
JC69G	12	-1993.8406	4064.3035	0.0000	1.0000
JC69IG	13	-1993.1796	4069.3668	0.0000	1.0000
JC69	11	-2003.8695	4077.9762	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 683-tomodel-phylip/chr5_11334.phylip = 1 - 593;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11334.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 683-tomodel-phylip/chr5_11334.phylip = 1 - 593;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11334.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 (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_11334.phylip = 1 - 593;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11334.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;


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

