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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3923.0109	7886.0217	0.6572	0.6572
GTRIG	21	-3923.0302	7888.0604	0.2371	0.8943
GTRI	20	-3924.8388	7889.6777	0.1056	1.0000
HKYG	16	-3938.1809	7908.3618	0.0000	1.0000
HKYIG	17	-3938.2486	7910.4972	0.0000	1.0000
HKYI	16	-3939.5206	7911.0413	0.0000	1.0000
SYMG	17	-3945.4448	7924.8896	0.0000	1.0000
SYMIG	18	-3945.4832	7926.9665	0.0000	1.0000
SYMI	17	-3946.6888	7927.3775	0.0000	1.0000
K2PG	13	-3952.2404	7930.4807	0.0000	1.0000
K2PIG	14	-3952.3029	7932.6059	0.0000	1.0000
K2PI	13	-3953.4151	7932.8301	0.0000	1.0000
GTR	19	-3964.1150	7966.2301	0.0000	1.0000
HKY	15	-3981.3684	7992.7368	0.0000	1.0000
SYM	16	-3988.2026	8008.4052	0.0000	1.0000
K2P	12	-3997.4870	8018.9740	0.0000	1.0000
F81G	15	-4025.2302	8080.4604	0.0000	1.0000
F81I	15	-4026.1944	8082.3888	0.0000	1.0000
F81IG	16	-4025.2557	8082.5114	0.0000	1.0000
JC69G	12	-4036.7872	8097.5744	0.0000	1.0000
JC69I	12	-4037.7369	8099.4739	0.0000	1.0000
JC69IG	13	-4036.8153	8099.6305	0.0000	1.0000
F81	14	-4066.6651	8161.3302	0.0000	1.0000
JC69	11	-4079.0773	8180.1547	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3923.0109	7987.1797	0.6049	0.6049
HKYG	16	-3938.1809	7989.2882	0.2108	0.8156
GTRI	20	-3924.8388	7990.8356	0.0972	0.9129
HKYI	16	-3939.5206	7991.9676	0.0552	0.9681
GTRIG	21	-3923.0302	7994.2763	0.0174	0.9855
K2PG	13	-3952.2404	7996.2334	0.0065	0.9920
HKYIG	17	-3938.2486	7996.4814	0.0058	0.9978
K2PI	13	-3953.4151	7998.5828	0.0020	0.9998
K2PIG	14	-3952.3029	8003.4165	0.0002	1.0000
SYMG	17	-3945.4448	8010.8738	0.0000	1.0000
SYMI	17	-3946.6888	8013.3618	0.0000	1.0000
SYMIG	18	-3945.4832	8018.0086	0.0000	1.0000
GTR	19	-3964.1150	8062.3302	0.0000	1.0000
HKY	15	-3981.3684	8068.6053	0.0000	1.0000
K2P	12	-3997.4870	8079.6688	0.0000	1.0000
SYM	16	-3988.2026	8089.3315	0.0000	1.0000
F81G	15	-4025.2302	8156.3289	0.0000	1.0000
F81I	15	-4026.1944	8158.2572	0.0000	1.0000
JC69G	12	-4036.7872	8158.2692	0.0000	1.0000
JC69I	12	-4037.7369	8160.1686	0.0000	1.0000
F81IG	16	-4025.2557	8163.4377	0.0000	1.0000
JC69IG	13	-4036.8153	8165.3832	0.0000	1.0000
F81	14	-4066.6651	8232.1407	0.0000	1.0000
JC69	11	-4079.0773	8235.7915	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/chr3_1300.phylip = 1 - 1162;
 Partition Dummy = 1:591-tomodel-phylip/chr3_1300.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/chr3_1300.phylip = 1 - 1162;
 Partition Dummy = 1:591-tomodel-phylip/chr3_1300.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 (GTRG)]

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


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

