
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:12 04/04/2011

Input data from file "683-tomodel-phylip/chr2_9409.phylip" (7 taxa, 960 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2342.3056	4724.6113	0.5240	0.5240
GTRI	20	-2342.9418	4725.8836	0.2774	0.8014
GTRIG	21	-2342.2939	4726.5878	0.1950	0.9964
HKYG	16	-2352.0629	4736.1259	0.0017	0.9981
HKYI	16	-2352.6060	4737.2121	0.0010	0.9990
HKYIG	17	-2352.0600	4738.1201	0.0006	0.9996
GTR	19	-2350.5793	4739.1587	0.0004	1.0000
HKY	15	-2359.3983	4748.7965	0.0000	1.0000
F81G	15	-2400.8443	4831.6885	0.0000	1.0000
F81I	15	-2401.3363	4832.6726	0.0000	1.0000
F81IG	16	-2400.8430	4833.6859	0.0000	1.0000
F81	14	-2407.7356	4843.4712	0.0000	1.0000
K2PG	13	-2437.7319	4901.4638	0.0000	1.0000
K2PI	13	-2438.2814	4902.5629	0.0000	1.0000
K2PIG	14	-2437.7265	4903.4530	0.0000	1.0000
SYMG	17	-2436.0247	4906.0494	0.0000	1.0000
SYMI	17	-2436.5235	4907.0469	0.0000	1.0000
SYMIG	18	-2436.0229	4908.0457	0.0000	1.0000
K2P	12	-2445.3755	4914.7509	0.0000	1.0000
SYM	16	-2443.4077	4918.8154	0.0000	1.0000
JC69G	12	-2475.0754	4974.1509	0.0000	1.0000
JC69I	12	-2475.5940	4975.1880	0.0000	1.0000
JC69IG	13	-2475.0731	4976.1462	0.0000	1.0000
JC69	11	-2482.4390	4986.8781	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2352.0629	4813.9968	0.6013	0.6013
HKYI	16	-2352.6060	4815.0830	0.3493	0.9507
HKYIG	17	-2352.0600	4820.8579	0.0195	0.9701
HKY	15	-2359.3983	4821.8005	0.0121	0.9823
GTRG	20	-2342.3056	4821.9499	0.0113	0.9936
GTRI	20	-2342.9418	4823.2223	0.0060	0.9995
GTRIG	21	-2342.2939	4828.7934	0.0004	0.9999
GTR	19	-2350.5793	4831.6304	0.0001	1.0000
F81G	15	-2400.8443	4904.6925	0.0000	1.0000
F81I	15	-2401.3363	4905.6766	0.0000	1.0000
F81IG	16	-2400.8430	4911.5569	0.0000	1.0000
F81	14	-2407.7356	4911.6083	0.0000	1.0000
K2PG	13	-2437.7319	4964.7339	0.0000	1.0000
K2PI	13	-2438.2814	4965.8330	0.0000	1.0000
K2PIG	14	-2437.7265	4971.5900	0.0000	1.0000
K2P	12	-2445.3755	4973.1541	0.0000	1.0000
SYMG	17	-2436.0247	4988.7873	0.0000	1.0000
SYMI	17	-2436.5235	4989.7848	0.0000	1.0000
SYMIG	18	-2436.0229	4995.6505	0.0000	1.0000
SYM	16	-2443.4077	4996.6863	0.0000	1.0000
JC69G	12	-2475.0754	5032.5541	0.0000	1.0000
JC69I	12	-2475.5940	5033.5912	0.0000	1.0000
JC69IG	13	-2475.0731	5039.4163	0.0000	1.0000
JC69	11	-2482.4390	5040.4143	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 683-tomodel-phylip/chr2_9409.phylip = 1 - 960;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9409.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 683-tomodel-phylip/chr2_9409.phylip = 1 - 960;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9409.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 683-tomodel-phylip/chr2_9409.phylip = 1 - 960;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9409.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

