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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-2947.1409	5928.2818	0.5279	0.5279
SYMIG	18	-2947.1049	5930.2098	0.2013	0.7292
GTRG	20	-2945.9337	5931.8674	0.0879	0.8171
K2PG	13	-2953.0586	5932.1171	0.0776	0.8947
GTRIG	21	-2945.9139	5933.8279	0.0330	0.9277
K2PIG	14	-2952.9803	5933.9606	0.0309	0.9585
SYMI	17	-2950.2755	5934.5510	0.0230	0.9815
HKYG	16	-2952.4841	5936.9681	0.0069	0.9884
K2PI	13	-2955.7577	5937.5154	0.0052	0.9936
GTRI	20	-2949.2530	5938.5060	0.0032	0.9968
HKYIG	17	-2952.3982	5938.7963	0.0028	0.9995
HKYI	16	-2955.1305	5942.2610	0.0005	1.0000
SYM	16	-2981.7584	5995.5167	0.0000	1.0000
GTR	19	-2980.5957	5999.1915	0.0000	1.0000
K2P	12	-2988.9663	6001.9325	0.0000	1.0000
HKY	15	-2988.3400	6006.6801	0.0000	1.0000
JC69G	12	-3000.5513	6025.1026	0.0000	1.0000
JC69IG	13	-3000.5303	6027.0605	0.0000	1.0000
F81G	15	-2999.5490	6029.0979	0.0000	1.0000
JC69I	12	-3002.6383	6029.2767	0.0000	1.0000
F81IG	16	-2999.5294	6031.0588	0.0000	1.0000
F81I	15	-3001.6605	6033.3209	0.0000	1.0000
JC69	11	-3032.3252	6086.6505	0.0000	1.0000
F81	14	-3031.2447	6090.4893	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2953.0586	5993.9723	0.9052	0.9052
K2PI	13	-2955.7577	5999.3706	0.0609	0.9661
K2PIG	14	-2952.9803	6000.5739	0.0334	0.9994
SYMG	17	-2947.1409	6009.1694	0.0005	0.9999
HKYG	16	-2952.4841	6013.0976	0.0001	1.0000
SYMI	17	-2950.2755	6015.4386	0.0000	1.0000
SYMIG	18	-2947.1049	6015.8555	0.0000	1.0000
HKYI	16	-2955.1305	6018.3906	0.0000	1.0000
HKYIG	17	-2952.3982	6019.6839	0.0000	1.0000
GTRG	20	-2945.9337	6027.0293	0.0000	1.0000
GTRI	20	-2949.2530	6033.6679	0.0000	1.0000
GTRIG	21	-2945.9139	6033.7479	0.0000	1.0000
K2P	12	-2988.9663	6059.0297	0.0000	1.0000
SYM	16	-2981.7584	6071.6463	0.0000	1.0000
HKY	15	-2988.3400	6078.0515	0.0000	1.0000
JC69G	12	-3000.5513	6082.1998	0.0000	1.0000
JC69I	12	-3002.6383	6086.3738	0.0000	1.0000
JC69IG	13	-3000.5303	6088.9157	0.0000	1.0000
GTR	19	-2980.5957	6089.5953	0.0000	1.0000
F81G	15	-2999.5490	6100.4693	0.0000	1.0000
F81I	15	-3001.6605	6104.6924	0.0000	1.0000
F81IG	16	-2999.5294	6107.1883	0.0000	1.0000
JC69	11	-3032.3252	6138.9895	0.0000	1.0000
F81	14	-3031.2447	6157.1026	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_9806.phylip = 1 - 861;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9806.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_9806.phylip = 1 - 861;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9806.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=Fixed(Equal);
END;


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

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

