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

Input data from file "591-tomodel-phylip/chr3_16833.phylip" (7 taxa, 672 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	-1882.5881	3806.4666	0.5988	0.5988
GTRG	20	-1883.2813	3807.8530	0.2994	0.8982
GTRIG	21	-1883.3931	3810.2077	0.0922	0.9904
HKYI	16	-1891.7422	3816.3149	0.0044	0.9948
HKYG	16	-1891.8502	3816.5309	0.0039	0.9987
HKYIG	17	-1891.8874	3818.7107	0.0013	1.0000
GTR	19	-1901.9969	3843.1595	0.0000	1.0000
HKY	15	-1909.0124	3848.7565	0.0000	1.0000
F81I	15	-1920.9457	3872.6231	0.0000	1.0000
F81G	15	-1921.1869	3873.1056	0.0000	1.0000
F81IG	16	-1920.9407	3874.7119	0.0000	1.0000
SYMI	17	-1922.7906	3880.5169	0.0000	1.0000
SYMG	17	-1923.3714	3881.6786	0.0000	1.0000
SYMIG	18	-1922.7959	3882.6392	0.0000	1.0000
K2PI	13	-1930.1996	3886.9524	0.0000	1.0000
K2PG	13	-1930.5670	3887.6873	0.0000	1.0000
K2PIG	14	-1930.2013	3889.0419	0.0000	1.0000
F81	14	-1935.4370	3899.5133	0.0000	1.0000
SYM	16	-1937.5151	3907.8608	0.0000	1.0000
K2P	12	-1944.8991	3914.2716	0.0000	1.0000
JC69I	12	-1952.8463	3930.1660	0.0000	1.0000
JC69G	12	-1953.2231	3930.9196	0.0000	1.0000
JC69IG	13	-1952.8490	3932.2512	0.0000	1.0000
JC69	11	-1966.3365	3955.0729	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1882.5881	3805.1763	0.5962	0.5962
GTRG	20	-1883.2813	3806.5627	0.2981	0.8943
GTRIG	21	-1883.3931	3808.7862	0.0981	0.9924
HKYI	16	-1891.7422	3815.4844	0.0034	0.9958
HKYG	16	-1891.8502	3815.7004	0.0031	0.9989
HKYIG	17	-1891.8874	3817.7749	0.0011	1.0000
GTR	19	-1901.9969	3841.9938	0.0000	1.0000
HKY	15	-1909.0124	3848.0248	0.0000	1.0000
F81I	15	-1920.9457	3871.8914	0.0000	1.0000
F81G	15	-1921.1869	3872.3739	0.0000	1.0000
F81IG	16	-1920.9407	3873.8814	0.0000	1.0000
SYMI	17	-1922.7906	3879.5812	0.0000	1.0000
SYMG	17	-1923.3714	3880.7428	0.0000	1.0000
SYMIG	18	-1922.7959	3881.5917	0.0000	1.0000
K2PI	13	-1930.1996	3886.3992	0.0000	1.0000
K2PG	13	-1930.5670	3887.1341	0.0000	1.0000
K2PIG	14	-1930.2013	3888.4026	0.0000	1.0000
F81	14	-1935.4370	3898.8741	0.0000	1.0000
SYM	16	-1937.5151	3907.0303	0.0000	1.0000
K2P	12	-1944.8991	3913.7982	0.0000	1.0000
JC69I	12	-1952.8463	3929.6925	0.0000	1.0000
JC69G	12	-1953.2231	3930.4462	0.0000	1.0000
JC69IG	13	-1952.8490	3931.6980	0.0000	1.0000
JC69	11	-1966.3365	3954.6729	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1891.7422	3887.6485	0.5095	0.5095
HKYG	16	-1891.8502	3887.8645	0.4573	0.9668
HKYIG	17	-1891.8874	3894.4493	0.0170	0.9838
GTRI	20	-1882.5881	3895.3814	0.0107	0.9945
GTRG	20	-1883.2813	3896.7678	0.0053	0.9998
GTRIG	21	-1883.3931	3903.5016	0.0002	1.0000
HKY	15	-1909.0124	3915.6787	0.0000	1.0000
GTR	19	-1901.9969	3927.6887	0.0000	1.0000
F81I	15	-1920.9457	3939.5453	0.0000	1.0000
F81G	15	-1921.1869	3940.0277	0.0000	1.0000
K2PI	13	-1930.1996	3945.0326	0.0000	1.0000
K2PG	13	-1930.5670	3945.7674	0.0000	1.0000
F81IG	16	-1920.9407	3946.0455	0.0000	1.0000
K2PIG	14	-1930.2013	3951.5462	0.0000	1.0000
SYMI	17	-1922.7906	3956.2556	0.0000	1.0000
SYMG	17	-1923.3714	3957.4172	0.0000	1.0000
F81	14	-1935.4370	3962.0177	0.0000	1.0000
SYMIG	18	-1922.7959	3962.7764	0.0000	1.0000
K2P	12	-1944.8991	3967.9213	0.0000	1.0000
SYM	16	-1937.5151	3979.1944	0.0000	1.0000
JC69I	12	-1952.8463	3983.8156	0.0000	1.0000
JC69G	12	-1953.2231	3984.5693	0.0000	1.0000
JC69IG	13	-1952.8490	3990.3314	0.0000	1.0000
JC69	11	-1966.3365	4004.2857	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_16833.phylip = 1 - 672;
 Partition Dummy = 1:591-tomodel-phylip/chr3_16833.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_16833.phylip = 1 - 672;
 Partition Dummy = 1:591-tomodel-phylip/chr3_16833.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_16833.phylip = 1 - 672;
 Partition Dummy = 1:591-tomodel-phylip/chr3_16833.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

