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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1930.1230	3892.9871	0.4036	0.4036
HKYI	16	-1930.2365	3893.2141	0.3603	0.7639
HKYIG	17	-1930.1284	3895.0917	0.1409	0.9048
HKY	15	-1933.8848	3898.4226	0.0266	0.9314
K2PG	13	-1936.4767	3899.4474	0.0160	0.9474
K2PI	13	-1936.5796	3899.6531	0.0144	0.9618
GTRG	20	-1929.3522	3899.8551	0.0130	0.9748
GTRI	20	-1929.4771	3900.1049	0.0115	0.9863
K2PIG	14	-1936.4816	3901.5339	0.0056	0.9919
GTRIG	21	-1929.3590	3901.9854	0.0045	0.9964
K2P	12	-1940.0791	3904.5810	0.0012	0.9976
GTR	19	-1933.0393	3905.1182	0.0009	0.9985
SYMG	17	-1935.5421	3905.9192	0.0006	0.9992
SYMI	17	-1935.6638	3906.1625	0.0006	0.9997
SYMIG	18	-1935.5492	3908.0329	0.0002	0.9999
SYM	16	-1939.0539	3910.8490	0.0001	1.0000
F81G	15	-1970.7128	3972.0788	0.0000	1.0000
F81I	15	-1970.7292	3972.1114	0.0000	1.0000
F81IG	16	-1970.7094	3974.1599	0.0000	1.0000
F81	14	-1973.3809	3975.3324	0.0000	1.0000
JC69G	12	-1976.3574	3977.1376	0.0000	1.0000
JC69I	12	-1976.3737	3977.1701	0.0000	1.0000
JC69IG	13	-1976.3543	3979.2025	0.0000	1.0000
JC69	11	-1978.9195	3980.1962	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1930.1230	3892.2459	0.4001	0.4001
HKYI	16	-1930.2365	3892.4730	0.3572	0.7573
HKYIG	17	-1930.1284	3894.2568	0.1464	0.9037
HKY	15	-1933.8848	3897.7696	0.0253	0.9290
GTRG	20	-1929.3522	3898.7044	0.0158	0.9449
K2PG	13	-1936.4767	3898.9535	0.0140	0.9589
GTRI	20	-1929.4771	3898.9543	0.0140	0.9728
K2PI	13	-1936.5796	3899.1592	0.0126	0.9855
GTRIG	21	-1929.3590	3900.7179	0.0058	0.9912
K2PIG	14	-1936.4816	3900.9633	0.0051	0.9964
GTR	19	-1933.0393	3904.0786	0.0011	0.9974
K2P	12	-1940.0791	3904.1582	0.0010	0.9985
SYMG	17	-1935.5421	3905.0842	0.0007	0.9991
SYMI	17	-1935.6638	3905.3276	0.0006	0.9997
SYMIG	18	-1935.5492	3907.0985	0.0002	0.9999
SYM	16	-1939.0539	3910.1078	0.0001	1.0000
F81G	15	-1970.7128	3971.4257	0.0000	1.0000
F81I	15	-1970.7292	3971.4584	0.0000	1.0000
F81IG	16	-1970.7094	3973.4187	0.0000	1.0000
F81	14	-1973.3809	3974.7617	0.0000	1.0000
JC69G	12	-1976.3574	3976.7148	0.0000	1.0000
JC69I	12	-1976.3737	3976.7474	0.0000	1.0000
JC69IG	13	-1976.3543	3978.7086	0.0000	1.0000
JC69	11	-1978.9195	3979.8389	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1936.4767	3959.0318	0.3627	0.3627
K2PI	13	-1936.5796	3959.2374	0.3272	0.6899
K2P	12	-1940.0791	3959.6151	0.2709	0.9608
K2PIG	14	-1936.4816	3965.6630	0.0132	0.9740
HKYG	16	-1930.1230	3966.1884	0.0101	0.9841
HKYI	16	-1930.2365	3966.4155	0.0090	0.9932
HKY	15	-1933.8848	3967.0906	0.0064	0.9996
HKYIG	17	-1930.1284	3972.8207	0.0004	1.0000
SYMG	17	-1935.5421	3983.6481	0.0000	1.0000
SYMI	17	-1935.6638	3983.8915	0.0000	1.0000
SYM	16	-1939.0539	3984.0503	0.0000	1.0000
SYMIG	18	-1935.5492	3990.2838	0.0000	1.0000
GTRG	20	-1929.3522	3991.1325	0.0000	1.0000
GTRI	20	-1929.4771	3991.3824	0.0000	1.0000
GTR	19	-1933.0393	3991.8853	0.0000	1.0000
GTRIG	21	-1929.3590	3997.7674	0.0000	1.0000
JC69	11	-1978.9195	4030.6744	0.0000	1.0000
JC69G	12	-1976.3574	4032.1717	0.0000	1.0000
JC69I	12	-1976.3737	4032.2042	0.0000	1.0000
JC69IG	13	-1976.3543	4038.7869	0.0000	1.0000
F81	14	-1973.3809	4039.4614	0.0000	1.0000
F81G	15	-1970.7128	4040.7468	0.0000	1.0000
F81I	15	-1970.7292	4040.7795	0.0000	1.0000
F81IG	16	-1970.7094	4047.3612	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/chrZ_10758.phylip = 1 - 751;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_10758.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/chrZ_10758.phylip = 1 - 751;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_10758.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 (K2PG)]

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

