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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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	-1887.0149	3815.1179	0.5428	0.5428
GTRIG	21	-1886.6731	3816.5447	0.2660	0.8088
GTRI	20	-1888.2032	3817.4945	0.1654	0.9742
HKYG	16	-1894.8834	3822.4679	0.0138	0.9880
HKYIG	17	-1894.6454	3824.0805	0.0061	0.9941
HKYI	16	-1895.7378	3824.1767	0.0059	1.0000
GTR	19	-1900.4536	3839.8903	0.0000	1.0000
HKY	15	-1907.6184	3845.8546	0.0000	1.0000
F81G	15	-1930.3495	3891.3168	0.0000	1.0000
F81I	15	-1930.9542	3892.5262	0.0000	1.0000
F81IG	16	-1930.2327	3893.1665	0.0000	1.0000
K2PG	13	-1936.6942	3899.8556	0.0000	1.0000
K2PI	13	-1937.3417	3901.1506	0.0000	1.0000
K2PIG	14	-1936.4667	3901.4733	0.0000	1.0000
SYMG	17	-1933.8755	3902.5406	0.0000	1.0000
SYMIG	18	-1933.5873	3904.0583	0.0000	1.0000
SYMI	17	-1934.6722	3904.1341	0.0000	1.0000
F81	14	-1941.9365	3912.4128	0.0000	1.0000
K2P	12	-1950.1403	3924.6806	0.0000	1.0000
SYM	16	-1947.3513	3927.4036	0.0000	1.0000
JC69G	12	-1967.4438	3959.2877	0.0000	1.0000
JC69I	12	-1967.9691	3960.3382	0.0000	1.0000
JC69IG	13	-1967.3312	3961.1296	0.0000	1.0000
JC69	11	-1979.5132	3981.3645	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1887.0149	3814.0298	0.5371	0.5371
GTRIG	21	-1886.6731	3815.3463	0.2781	0.8151
GTRI	20	-1888.2032	3816.4065	0.1637	0.9788
HKYG	16	-1894.8834	3821.7669	0.0112	0.9900
HKYIG	17	-1894.6454	3823.2909	0.0052	0.9952
HKYI	16	-1895.7378	3823.4757	0.0048	1.0000
GTR	19	-1900.4536	3838.9071	0.0000	1.0000
HKY	15	-1907.6184	3845.2368	0.0000	1.0000
F81G	15	-1930.3495	3890.6990	0.0000	1.0000
F81I	15	-1930.9542	3891.9084	0.0000	1.0000
F81IG	16	-1930.2327	3892.4655	0.0000	1.0000
K2PG	13	-1936.6942	3899.3883	0.0000	1.0000
K2PI	13	-1937.3417	3900.6834	0.0000	1.0000
K2PIG	14	-1936.4667	3900.9335	0.0000	1.0000
SYMG	17	-1933.8755	3901.7510	0.0000	1.0000
SYMIG	18	-1933.5873	3903.1746	0.0000	1.0000
SYMI	17	-1934.6722	3903.3444	0.0000	1.0000
F81	14	-1941.9365	3911.8730	0.0000	1.0000
K2P	12	-1950.1403	3924.2806	0.0000	1.0000
SYM	16	-1947.3513	3926.7025	0.0000	1.0000
JC69G	12	-1967.4438	3958.8877	0.0000	1.0000
JC69I	12	-1967.9691	3959.9382	0.0000	1.0000
JC69IG	13	-1967.3312	3960.6623	0.0000	1.0000
JC69	11	-1979.5132	3981.0264	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1894.8834	3896.5801	0.6774	0.6774
HKYI	16	-1895.7378	3898.2889	0.2882	0.9656
HKYIG	17	-1894.6454	3902.7799	0.0305	0.9961
GTRG	20	-1887.0149	3907.5462	0.0028	0.9989
GTRI	20	-1888.2032	3909.9229	0.0009	0.9998
GTRIG	21	-1886.6731	3913.5385	0.0001	0.9999
HKY	15	-1907.6184	3915.3742	0.0001	1.0000
GTR	19	-1900.4536	3927.7478	0.0000	1.0000
K2PG	13	-1936.6942	3960.1740	0.0000	1.0000
F81G	15	-1930.3495	3960.8364	0.0000	1.0000
K2PI	13	-1937.3417	3961.4691	0.0000	1.0000
F81I	15	-1930.9542	3962.0458	0.0000	1.0000
K2PIG	14	-1936.4667	3966.3950	0.0000	1.0000
F81IG	16	-1930.2327	3967.2787	0.0000	1.0000
F81	14	-1941.9365	3977.3345	0.0000	1.0000
K2P	12	-1950.1403	3980.3904	0.0000	1.0000
SYMG	17	-1933.8755	3981.2400	0.0000	1.0000
SYMI	17	-1934.6722	3982.8334	0.0000	1.0000
SYMIG	18	-1933.5873	3987.3394	0.0000	1.0000
SYM	16	-1947.3513	4001.5157	0.0000	1.0000
JC69G	12	-1967.4438	4014.9976	0.0000	1.0000
JC69I	12	-1967.9691	4016.0481	0.0000	1.0000
JC69IG	13	-1967.3312	4021.4480	0.0000	1.0000
JC69	11	-1979.5132	4032.4605	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/chr7_9094.phylip = 1 - 793;
 Partition Dummy = 1:591-tomodel-phylip/chr7_9094.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/chr7_9094.phylip = 1 - 793;
 Partition Dummy = 1:591-tomodel-phylip/chr7_9094.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_9094.phylip = 1 - 793;
 Partition Dummy = 1:591-tomodel-phylip/chr7_9094.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

