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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1909.9963	3852.8182	0.2978	0.2978
GTRG	20	-1906.1789	3853.6402	0.1974	0.4952
HKYI	16	-1910.4153	3853.6560	0.1959	0.6910
GTRI	20	-1906.5860	3854.4545	0.1314	0.8224
HKYIG	17	-1909.9650	3854.8602	0.1073	0.9297
GTRIG	21	-1906.1460	3855.7049	0.0703	1.0000
HKY	15	-1932.9256	3896.5785	0.0000	1.0000
GTR	19	-1930.0338	3899.2261	0.0000	1.0000
K2PG	13	-1942.4726	3911.4950	0.0000	1.0000
F81G	15	-1940.3863	3911.4998	0.0000	1.0000
F81I	15	-1940.4353	3911.5980	0.0000	1.0000
K2PI	13	-1942.7965	3912.1428	0.0000	1.0000
K2PIG	14	-1942.4614	3913.5581	0.0000	1.0000
F81IG	16	-1940.4200	3913.6654	0.0000	1.0000
SYMG	17	-1940.8052	3916.5406	0.0000	1.0000
SYMI	17	-1941.1804	3917.2910	0.0000	1.0000
SYMIG	18	-1940.7760	3918.5931	0.0000	1.0000
F81	14	-1961.1233	3950.8820	0.0000	1.0000
K2P	12	-1963.6847	3951.8400	0.0000	1.0000
SYM	16	-1961.6614	3956.1483	0.0000	1.0000
JC69G	12	-1967.6971	3959.8648	0.0000	1.0000
JC69I	12	-1967.7516	3959.9737	0.0000	1.0000
JC69IG	13	-1967.7428	3962.0354	0.0000	1.0000
JC69	11	-1987.7736	3997.9448	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1909.9963	3851.9927	0.2673	0.2673
GTRG	20	-1906.1789	3852.3577	0.2227	0.4899
HKYI	16	-1910.4153	3852.8305	0.1758	0.6657
GTRI	20	-1906.5860	3853.1721	0.1482	0.8139
HKYIG	17	-1909.9650	3853.9301	0.1014	0.9154
GTRIG	21	-1906.1460	3854.2921	0.0846	1.0000
HKY	15	-1932.9256	3895.8512	0.0000	1.0000
GTR	19	-1930.0338	3898.0675	0.0000	1.0000
F81G	15	-1940.3863	3910.7725	0.0000	1.0000
F81I	15	-1940.4353	3910.8707	0.0000	1.0000
K2PG	13	-1942.4726	3910.9451	0.0000	1.0000
K2PI	13	-1942.7965	3911.5930	0.0000	1.0000
F81IG	16	-1940.4200	3912.8399	0.0000	1.0000
K2PIG	14	-1942.4614	3912.9227	0.0000	1.0000
SYMG	17	-1940.8052	3915.6105	0.0000	1.0000
SYMI	17	-1941.1804	3916.3609	0.0000	1.0000
SYMIG	18	-1940.7760	3917.5520	0.0000	1.0000
F81	14	-1961.1233	3950.2466	0.0000	1.0000
K2P	12	-1963.6847	3951.3694	0.0000	1.0000
SYM	16	-1961.6614	3955.3228	0.0000	1.0000
JC69G	12	-1967.6971	3959.3942	0.0000	1.0000
JC69I	12	-1967.7516	3959.5032	0.0000	1.0000
JC69IG	13	-1967.7428	3961.4855	0.0000	1.0000
JC69	11	-1987.7736	3997.5472	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1909.9963	3924.2518	0.5891	0.5891
HKYI	16	-1910.4153	3925.0896	0.3875	0.9765
HKYIG	17	-1909.9650	3930.7054	0.0234	0.9999
GTRG	20	-1906.1789	3942.6816	0.0001	1.0000
GTRI	20	-1906.5860	3943.4959	0.0000	1.0000
GTRIG	21	-1906.1460	3949.1321	0.0000	1.0000
HKY	15	-1932.9256	3963.5941	0.0000	1.0000
K2PG	13	-1942.4726	3969.6556	0.0000	1.0000
K2PI	13	-1942.7965	3970.3035	0.0000	1.0000
K2PIG	14	-1942.4614	3976.1494	0.0000	1.0000
F81G	15	-1940.3863	3978.5154	0.0000	1.0000
F81I	15	-1940.4353	3978.6136	0.0000	1.0000
GTR	19	-1930.0338	3983.8752	0.0000	1.0000
F81IG	16	-1940.4200	3985.0990	0.0000	1.0000
SYMG	17	-1940.8052	3992.3857	0.0000	1.0000
SYMI	17	-1941.1804	3993.1361	0.0000	1.0000
SYMIG	18	-1940.7760	3998.8434	0.0000	1.0000
K2P	12	-1963.6847	4005.5637	0.0000	1.0000
F81	14	-1961.1233	4013.4733	0.0000	1.0000
JC69G	12	-1967.6971	4013.5885	0.0000	1.0000
JC69I	12	-1967.7516	4013.6975	0.0000	1.0000
JC69IG	13	-1967.7428	4020.1960	0.0000	1.0000
SYM	16	-1961.6614	4027.5819	0.0000	1.0000
JC69	11	-1987.7736	4047.2254	0.0000	1.0000



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

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

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


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

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


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

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

