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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1908.6333	3858.8193	0.3817	0.3817
GTRI	20	-1909.4315	3860.4157	0.1718	0.5535
SYMG	17	-1912.6756	3860.4763	0.1667	0.7202
GTRIG	21	-1908.6568	3861.0246	0.1267	0.8469
SYMI	17	-1913.2658	3861.6566	0.0924	0.9392
SYMIG	18	-1912.7134	3862.6865	0.0552	0.9944
K2PG	13	-1921.4492	3869.5627	0.0018	0.9962
K2PI	13	-1921.8572	3870.3786	0.0012	0.9974
HKYG	16	-1918.8564	3870.7109	0.0010	0.9984
HKYI	16	-1919.2705	3871.5392	0.0007	0.9991
K2PIG	14	-1921.4702	3871.7082	0.0006	0.9997
HKYIG	17	-1918.8706	3872.8662	0.0003	1.0000
GTR	19	-1940.2801	3919.9624	0.0000	1.0000
SYM	16	-1943.7097	3920.4175	0.0000	1.0000
K2P	12	-1951.6428	3927.8539	0.0000	1.0000
JC69G	12	-1951.8650	3928.2984	0.0000	1.0000
HKY	15	-1949.0687	3929.0165	0.0000	1.0000
JC69I	12	-1952.4417	3929.4516	0.0000	1.0000
F81G	15	-1949.6117	3930.1026	0.0000	1.0000
JC69IG	13	-1951.8901	3930.4444	0.0000	1.0000
F81I	15	-1950.1217	3931.1225	0.0000	1.0000
F81IG	16	-1949.6358	3932.2697	0.0000	1.0000
JC69	11	-1979.8334	3982.1468	0.0000	1.0000
F81	14	-1977.5505	3983.8688	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1908.6333	3857.2666	0.4012	0.4012
GTRI	20	-1909.4315	3858.8630	0.1806	0.5818
GTRIG	21	-1908.6568	3859.3135	0.1442	0.7260
SYMG	17	-1912.6756	3859.3513	0.1415	0.8674
SYMI	17	-1913.2658	3860.5316	0.0784	0.9459
SYMIG	18	-1912.7134	3861.4268	0.0501	0.9960
K2PG	13	-1921.4492	3868.8984	0.0012	0.9972
HKYG	16	-1918.8564	3869.7128	0.0008	0.9980
K2PI	13	-1921.8572	3869.7144	0.0008	0.9988
HKYI	16	-1919.2705	3870.5411	0.0005	0.9993
K2PIG	14	-1921.4702	3870.9403	0.0004	0.9997
HKYIG	17	-1918.8706	3871.7412	0.0003	1.0000
GTR	19	-1940.2801	3918.5602	0.0000	1.0000
SYM	16	-1943.7097	3919.4193	0.0000	1.0000
K2P	12	-1951.6428	3927.2856	0.0000	1.0000
JC69G	12	-1951.8650	3927.7301	0.0000	1.0000
HKY	15	-1949.0687	3928.1374	0.0000	1.0000
JC69I	12	-1952.4417	3928.8833	0.0000	1.0000
F81G	15	-1949.6117	3929.2234	0.0000	1.0000
JC69IG	13	-1951.8901	3929.7802	0.0000	1.0000
F81I	15	-1950.1217	3930.2434	0.0000	1.0000
F81IG	16	-1949.6358	3931.2716	0.0000	1.0000
JC69	11	-1979.8334	3981.6668	0.0000	1.0000
F81	14	-1977.5505	3983.1010	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1921.4492	3925.2080	0.5745	0.5745
K2PI	13	-1921.8572	3926.0239	0.3820	0.9565
K2PIG	14	-1921.4702	3931.5814	0.0237	0.9802
SYMG	17	-1912.6756	3932.9868	0.0118	0.9920
SYMI	17	-1913.2658	3934.1672	0.0065	0.9985
HKYG	16	-1918.8564	3939.0168	0.0006	0.9990
SYMIG	18	-1912.7134	3939.3939	0.0005	0.9995
HKYI	16	-1919.2705	3939.8451	0.0004	0.9999
GTRG	20	-1908.6333	3943.8967	0.0001	1.0000
HKYIG	17	-1918.8706	3945.3767	0.0000	1.0000
GTRI	20	-1909.4315	3945.4931	0.0000	1.0000
GTRIG	21	-1908.6568	3950.2750	0.0000	1.0000
K2P	12	-1951.6428	3979.2637	0.0000	1.0000
JC69G	12	-1951.8650	3979.7081	0.0000	1.0000
JC69I	12	-1952.4417	3980.8613	0.0000	1.0000
JC69IG	13	-1951.8901	3986.0897	0.0000	1.0000
SYM	16	-1943.7097	3988.7233	0.0000	1.0000
HKY	15	-1949.0687	3993.1099	0.0000	1.0000
F81G	15	-1949.6117	3994.1960	0.0000	1.0000
F81I	15	-1950.1217	3995.2159	0.0000	1.0000
F81IG	16	-1949.6358	4000.5756	0.0000	1.0000
GTR	19	-1940.2801	4000.8587	0.0000	1.0000
JC69	11	-1979.8334	4029.3133	0.0000	1.0000
F81	14	-1977.5505	4043.7420	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_5495.phylip = 1 - 562;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_5495.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_5495.phylip = 1 - 562;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_5495.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

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

