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

Input data from file "591-tomodel-phylip/chrZ_4787.phylip" (7 taxa, 1128 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	-3710.9499	7461.8999	0.6654	0.6654
GTRIG	21	-3710.8141	7463.6283	0.2804	0.9458
GTRI	20	-3713.8629	7467.7259	0.0361	0.9819
HKYG	16	-3718.9294	7469.8589	0.0124	0.9943
HKYIG	17	-3718.8460	7471.6919	0.0050	0.9993
HKYI	16	-3721.8222	7475.6445	0.0007	1.0000
SYMG	17	-3730.9707	7495.9413	0.0000	1.0000
SYMIG	18	-3730.8374	7497.6748	0.0000	1.0000
K2PG	13	-3735.9147	7497.8294	0.0000	1.0000
K2PIG	14	-3735.8255	7499.6509	0.0000	1.0000
SYMI	17	-3734.0727	7502.1454	0.0000	1.0000
K2PI	13	-3738.7773	7503.5545	0.0000	1.0000
GTR	19	-3755.1420	7548.2840	0.0000	1.0000
HKY	15	-3763.2604	7556.5208	0.0000	1.0000
SYM	16	-3775.5408	7583.0816	0.0000	1.0000
K2P	12	-3781.1419	7586.2838	0.0000	1.0000
F81G	15	-3804.6116	7639.2232	0.0000	1.0000
F81IG	16	-3804.6232	7641.2465	0.0000	1.0000
F81I	15	-3806.4005	7642.8010	0.0000	1.0000
JC69G	12	-3819.3354	7662.6709	0.0000	1.0000
JC69IG	13	-3819.3441	7664.6882	0.0000	1.0000
JC69I	12	-3821.3395	7666.6790	0.0000	1.0000
F81	14	-3841.5732	7711.1463	0.0000	1.0000
JC69	11	-3856.8601	7735.7203	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3718.9294	7550.3101	0.9157	0.9157
HKYI	16	-3721.8222	7556.0957	0.0508	0.9665
HKYIG	17	-3718.8460	7557.1713	0.0296	0.9961
GTRG	20	-3710.9499	7562.4639	0.0021	0.9982
K2PG	13	-3735.9147	7563.1960	0.0015	0.9997
GTRI	20	-3713.8629	7568.2899	0.0001	0.9998
K2PI	13	-3738.7773	7568.9212	0.0001	0.9999
GTRIG	21	-3710.8141	7569.2205	0.0001	1.0000
K2PIG	14	-3735.8255	7570.0457	0.0000	1.0000
SYMG	17	-3730.9707	7581.4208	0.0000	1.0000
SYMI	17	-3734.0727	7587.6248	0.0000	1.0000
SYMIG	18	-3730.8374	7588.1825	0.0000	1.0000
HKY	15	-3763.2604	7631.9438	0.0000	1.0000
GTR	19	-3755.1420	7643.8198	0.0000	1.0000
K2P	12	-3781.1419	7646.6223	0.0000	1.0000
SYM	16	-3775.5408	7663.5328	0.0000	1.0000
F81G	15	-3804.6116	7714.6462	0.0000	1.0000
F81I	15	-3806.4005	7718.2240	0.0000	1.0000
F81IG	16	-3804.6232	7721.6977	0.0000	1.0000
JC69G	12	-3819.3354	7723.0093	0.0000	1.0000
JC69I	12	-3821.3395	7727.0174	0.0000	1.0000
JC69IG	13	-3819.3441	7730.0548	0.0000	1.0000
F81	14	-3841.5732	7781.5411	0.0000	1.0000
JC69	11	-3856.8601	7791.0305	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/chrZ_4787.phylip = 1 - 1128;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4787.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/chrZ_4787.phylip = 1 - 1128;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4787.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/chrZ_4787.phylip = 1 - 1128;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4787.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

