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

Input data from file "591-tomodel-phylip/chr6_8057.phylip" (7 taxa, 1124 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	-3742.7420	7525.4841	0.6861	0.6861
GTRIG	21	-3742.6765	7527.3530	0.2695	0.9556
HKYG	16	-3750.1512	7532.3025	0.0227	0.9782
GTRI	20	-3746.7769	7533.5537	0.0121	0.9904
HKYIG	17	-3750.0732	7534.1465	0.0090	0.9994
HKYI	16	-3753.7955	7539.5910	0.0006	1.0000
GTR	19	-3791.2957	7620.5914	0.0000	1.0000
HKY	15	-3798.6527	7627.3053	0.0000	1.0000
SYMG	17	-3797.2224	7628.4448	0.0000	1.0000
K2PG	13	-3801.7096	7629.4191	0.0000	1.0000
SYMIG	18	-3797.1349	7630.2698	0.0000	1.0000
K2PIG	14	-3801.6333	7631.2665	0.0000	1.0000
SYMI	17	-3799.7743	7633.5485	0.0000	1.0000
K2PI	13	-3804.3851	7634.7702	0.0000	1.0000
F81G	15	-3836.6064	7703.2128	0.0000	1.0000
F81IG	16	-3836.5407	7705.0813	0.0000	1.0000
F81I	15	-3839.0402	7708.0804	0.0000	1.0000
SYM	16	-3839.5587	7711.1174	0.0000	1.0000
K2P	12	-3845.7763	7715.5526	0.0000	1.0000
JC69G	12	-3874.4027	7772.8054	0.0000	1.0000
JC69IG	13	-3874.3201	7774.6403	0.0000	1.0000
JC69I	12	-3876.6025	7777.2050	0.0000	1.0000
F81	14	-3879.2958	7786.5915	0.0000	1.0000
JC69	11	-3916.0554	7854.1108	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3750.1512	7612.6968	0.9436	0.9436
HKYIG	17	-3750.0732	7619.5655	0.0304	0.9740
HKYI	16	-3753.7955	7619.9853	0.0247	0.9987
GTRG	20	-3742.7420	7625.9771	0.0012	0.9999
GTRIG	21	-3742.6765	7632.8707	0.0000	1.0000
GTRI	20	-3746.7769	7634.0467	0.0000	1.0000
K2PG	13	-3801.7096	7694.7396	0.0000	1.0000
K2PI	13	-3804.3851	7700.0906	0.0000	1.0000
K2PIG	14	-3801.6333	7701.6116	0.0000	1.0000
HKY	15	-3798.6527	7702.6751	0.0000	1.0000
SYMG	17	-3797.2224	7713.8638	0.0000	1.0000
GTR	19	-3791.2957	7716.0598	0.0000	1.0000
SYMI	17	-3799.7743	7718.9676	0.0000	1.0000
SYMIG	18	-3797.1349	7720.7134	0.0000	1.0000
K2P	12	-3845.7763	7775.8484	0.0000	1.0000
F81G	15	-3836.6064	7778.5825	0.0000	1.0000
F81I	15	-3839.0402	7783.4502	0.0000	1.0000
F81IG	16	-3836.5407	7785.4757	0.0000	1.0000
SYM	16	-3839.5587	7791.5118	0.0000	1.0000
JC69G	12	-3874.4027	7833.1012	0.0000	1.0000
JC69I	12	-3876.6025	7837.5008	0.0000	1.0000
JC69IG	13	-3874.3201	7839.9607	0.0000	1.0000
F81	14	-3879.2958	7856.9366	0.0000	1.0000
JC69	11	-3916.0554	7909.3820	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/chr6_8057.phylip = 1 - 1124;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8057.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/chr6_8057.phylip = 1 - 1124;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8057.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/chr6_8057.phylip = 1 - 1124;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8057.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

