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

Input data from file "591-tomodel-phylip/chr8_8812.phylip" (7 taxa, 677 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	-2431.4984	4904.2773	0.6784	0.6784
GTRIG	21	-2431.4680	4906.3467	0.2410	0.9194
GTRI	20	-2433.7889	4908.8582	0.0687	0.9881
HKYG	16	-2440.2911	4913.4064	0.0071	0.9952
HKYIG	17	-2440.1430	4915.2146	0.0029	0.9980
HKYI	16	-2441.5591	4915.9425	0.0020	1.0000
SYMI	17	-2451.3673	4937.6632	0.0000	1.0000
SYMG	17	-2451.4413	4937.8113	0.0000	1.0000
K2PG	13	-2455.9803	4938.5097	0.0000	1.0000
SYMIG	18	-2450.8109	4938.6614	0.0000	1.0000
K2PI	13	-2456.6047	4939.7583	0.0000	1.0000
K2PIG	14	-2455.7197	4940.0738	0.0000	1.0000
GTR	19	-2476.1766	4991.5101	0.0000	1.0000
HKY	15	-2482.6695	4996.0652	0.0000	1.0000
F81G	15	-2489.1018	5008.9297	0.0000	1.0000
F81I	15	-2489.7958	5010.3177	0.0000	1.0000
F81IG	16	-2489.0142	5010.8527	0.0000	1.0000
SYM	16	-2493.3781	5019.5804	0.0000	1.0000
K2P	12	-2497.7733	5020.0165	0.0000	1.0000
JC69G	12	-2502.9007	5030.2714	0.0000	1.0000
JC69IG	13	-2501.9959	5030.5409	0.0000	1.0000
JC69I	12	-2503.0553	5030.5805	0.0000	1.0000
F81	14	-2528.8259	5086.2862	0.0000	1.0000
JC69	11	-2542.0081	5106.4131	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2431.4984	4902.9969	0.6691	0.6691
GTRIG	21	-2431.4680	4904.9360	0.2537	0.9228
GTRI	20	-2433.7889	4907.5777	0.0677	0.9905
HKYG	16	-2440.2911	4912.5822	0.0055	0.9961
HKYIG	17	-2440.1430	4914.2860	0.0024	0.9984
HKYI	16	-2441.5591	4915.1182	0.0016	1.0000
SYMI	17	-2451.3673	4936.7345	0.0000	1.0000
SYMG	17	-2451.4413	4936.8827	0.0000	1.0000
SYMIG	18	-2450.8109	4937.6219	0.0000	1.0000
K2PG	13	-2455.9803	4937.9607	0.0000	1.0000
K2PI	13	-2456.6047	4939.2093	0.0000	1.0000
K2PIG	14	-2455.7197	4939.4393	0.0000	1.0000
GTR	19	-2476.1766	4990.3533	0.0000	1.0000
HKY	15	-2482.6695	4995.3390	0.0000	1.0000
F81G	15	-2489.1018	5008.2035	0.0000	1.0000
F81I	15	-2489.7958	5009.5915	0.0000	1.0000
F81IG	16	-2489.0142	5010.0285	0.0000	1.0000
SYM	16	-2493.3781	5018.7562	0.0000	1.0000
K2P	12	-2497.7733	5019.5466	0.0000	1.0000
JC69G	12	-2502.9007	5029.8015	0.0000	1.0000
JC69IG	13	-2501.9959	5029.9919	0.0000	1.0000
JC69I	12	-2503.0553	5030.1106	0.0000	1.0000
F81	14	-2528.8259	5085.6518	0.0000	1.0000
JC69	11	-2542.0081	5106.0161	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2440.2911	4984.8649	0.7426	0.7426
HKYI	16	-2441.5591	4987.4010	0.2090	0.9515
HKYIG	17	-2440.1430	4991.0864	0.0331	0.9846
GTRG	20	-2431.4984	4993.3503	0.0107	0.9953
K2PG	13	-2455.9803	4996.6904	0.0020	0.9973
GTRI	20	-2433.7889	4997.9312	0.0011	0.9984
K2PI	13	-2456.6047	4997.9390	0.0011	0.9995
GTRIG	21	-2431.4680	4999.8071	0.0004	0.9999
K2PIG	14	-2455.7197	5002.6867	0.0001	1.0000
SYMI	17	-2451.3673	5013.5349	0.0000	1.0000
SYMG	17	-2451.4413	5013.6831	0.0000	1.0000
SYMIG	18	-2450.8109	5018.9400	0.0000	1.0000
HKY	15	-2482.6695	5063.1040	0.0000	1.0000
K2P	12	-2497.7733	5073.7586	0.0000	1.0000
F81G	15	-2489.1018	5075.9686	0.0000	1.0000
GTR	19	-2476.1766	5076.1890	0.0000	1.0000
F81I	15	-2489.7958	5077.3566	0.0000	1.0000
F81IG	16	-2489.0142	5082.3112	0.0000	1.0000
JC69G	12	-2502.9007	5084.0135	0.0000	1.0000
JC69I	12	-2503.0553	5084.3227	0.0000	1.0000
JC69IG	13	-2501.9959	5088.7216	0.0000	1.0000
SYM	16	-2493.3781	5091.0389	0.0000	1.0000
F81	14	-2528.8259	5148.8992	0.0000	1.0000
JC69	11	-2542.0081	5155.7105	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/chr8_8812.phylip = 1 - 677;
 Partition Dummy = 1:591-tomodel-phylip/chr8_8812.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/chr8_8812.phylip = 1 - 677;
 Partition Dummy = 1:591-tomodel-phylip/chr8_8812.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/chr8_8812.phylip = 1 - 677;
 Partition Dummy = 1:591-tomodel-phylip/chr8_8812.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

