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

Input data from file "591-tomodel-phylip/chrUn_random_9893.phylip" (7 taxa, 869 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	-2713.3756	5466.7512	0.4922	0.4922
GTRI	20	-2713.9655	5467.9311	0.2728	0.7650
GTRIG	21	-2713.3756	5468.7512	0.1811	0.9460
HKYG	16	-2720.1368	5472.2737	0.0311	0.9771
HKYIG	17	-2720.1368	5474.2737	0.0114	0.9886
HKYI	16	-2721.1401	5474.2803	0.0114	1.0000
GTR	19	-2766.5195	5571.0391	0.0000	1.0000
HKY	15	-2773.5587	5577.1174	0.0000	1.0000
F81G	15	-2784.9914	5599.9829	0.0000	1.0000
F81I	15	-2785.6928	5601.3856	0.0000	1.0000
F81IG	16	-2784.9915	5601.9829	0.0000	1.0000
K2PG	13	-2790.2788	5606.5576	0.0000	1.0000
SYMG	17	-2786.3599	5606.7198	0.0000	1.0000
K2PIG	14	-2790.2699	5608.5397	0.0000	1.0000
SYMIG	18	-2786.3539	5608.7078	0.0000	1.0000
K2PI	13	-2791.8897	5609.7793	0.0000	1.0000
SYMI	17	-2787.9191	5609.8382	0.0000	1.0000
F81	14	-2834.0868	5696.1736	0.0000	1.0000
JC69G	12	-2843.5024	5711.0049	0.0000	1.0000
JC69IG	13	-2843.5055	5713.0110	0.0000	1.0000
JC69I	12	-2844.6530	5713.3060	0.0000	1.0000
SYM	16	-2842.1444	5716.2889	0.0000	1.0000
K2P	12	-2846.3484	5716.6969	0.0000	1.0000
JC69	11	-2894.4929	5810.9858	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2720.1368	5548.5512	0.7131	0.7131
HKYI	16	-2721.1401	5550.5578	0.2615	0.9745
HKYIG	17	-2720.1368	5555.3185	0.0242	0.9987
GTRG	20	-2713.3756	5562.0981	0.0008	0.9995
GTRI	20	-2713.9655	5563.2779	0.0005	1.0000
GTRIG	21	-2713.3756	5568.8654	0.0000	1.0000
HKY	15	-2773.5587	5648.6276	0.0000	1.0000
GTR	19	-2766.5195	5661.6186	0.0000	1.0000
K2PG	13	-2790.2788	5668.5330	0.0000	1.0000
F81G	15	-2784.9914	5671.4930	0.0000	1.0000
K2PI	13	-2791.8897	5671.7548	0.0000	1.0000
F81I	15	-2785.6928	5672.8957	0.0000	1.0000
K2PIG	14	-2790.2699	5675.2825	0.0000	1.0000
F81IG	16	-2784.9915	5678.2604	0.0000	1.0000
SYMG	17	-2786.3599	5687.7646	0.0000	1.0000
SYMI	17	-2787.9191	5690.8831	0.0000	1.0000
SYMIG	18	-2786.3539	5694.5200	0.0000	1.0000
F81	14	-2834.0868	5762.9164	0.0000	1.0000
JC69G	12	-2843.5024	5768.2130	0.0000	1.0000
JC69I	12	-2844.6530	5770.5141	0.0000	1.0000
K2P	12	-2846.3484	5773.9050	0.0000	1.0000
JC69IG	13	-2843.5055	5774.9865	0.0000	1.0000
SYM	16	-2842.1444	5792.5664	0.0000	1.0000
JC69	11	-2894.4929	5863.4266	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/chrUn_random_9893.phylip = 1 - 869;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_9893.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/chrUn_random_9893.phylip = 1 - 869;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_9893.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/chrUn_random_9893.phylip = 1 - 869;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_9893.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

