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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3043.7769	6127.5537	0.5273	0.5273
GTRI	20	-3044.4078	6128.8156	0.2806	0.8079
GTRIG	21	-3043.7920	6129.5840	0.1911	0.9990
HKYI	16	-3054.8824	6141.7648	0.0004	0.9994
HKYG	16	-3054.9153	6141.8305	0.0004	0.9998
HKYIG	17	-3054.7436	6143.4872	0.0002	1.0000
GTR	19	-3062.8822	6163.7643	0.0000	1.0000
HKY	15	-3072.1647	6174.3294	0.0000	1.0000
K2PI	13	-3120.7303	6267.4606	0.0000	1.0000
SYMI	17	-3117.7606	6269.5212	0.0000	1.0000
K2PIG	14	-3120.8029	6269.6057	0.0000	1.0000
K2PG	13	-3121.8101	6269.6202	0.0000	1.0000
SYMG	17	-3119.1255	6272.2511	0.0000	1.0000
SYMIG	18	-3118.2652	6272.5304	0.0000	1.0000
K2P	12	-3137.0537	6298.1074	0.0000	1.0000
SYM	16	-3134.2088	6300.4176	0.0000	1.0000
F81G	15	-3156.7097	6343.4194	0.0000	1.0000
F81I	15	-3157.0013	6344.0027	0.0000	1.0000
F81IG	16	-3157.0055	6346.0110	0.0000	1.0000
F81	14	-3170.3520	6368.7040	0.0000	1.0000
JC69G	12	-3203.8753	6431.7507	0.0000	1.0000
JC69I	12	-3203.9654	6431.9308	0.0000	1.0000
JC69IG	13	-3203.9765	6433.9529	0.0000	1.0000
JC69	11	-3217.3466	6456.6932	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3054.8824	6219.1277	0.4704	0.4704
HKYG	16	-3054.9153	6219.1935	0.4552	0.9257
GTRG	20	-3043.7769	6224.2574	0.0362	0.9619
GTRI	20	-3044.4078	6225.5193	0.0193	0.9811
HKYIG	17	-3054.7436	6225.6853	0.0177	0.9988
GTRIG	21	-3043.7920	6231.1229	0.0012	1.0000
HKY	15	-3072.1647	6246.8571	0.0000	1.0000
GTR	19	-3062.8822	6255.6328	0.0000	1.0000
K2PI	13	-3120.7303	6330.3180	0.0000	1.0000
K2PG	13	-3121.8101	6332.4776	0.0000	1.0000
K2PIG	14	-3120.8029	6337.2983	0.0000	1.0000
SYMI	17	-3117.7606	6351.7193	0.0000	1.0000
SYMG	17	-3119.1255	6354.4492	0.0000	1.0000
K2P	12	-3137.0537	6356.1297	0.0000	1.0000
SYMIG	18	-3118.2652	6359.5637	0.0000	1.0000
SYM	16	-3134.2088	6377.7805	0.0000	1.0000
F81G	15	-3156.7097	6415.9472	0.0000	1.0000
F81I	15	-3157.0013	6416.5304	0.0000	1.0000
F81IG	16	-3157.0055	6423.3739	0.0000	1.0000
F81	14	-3170.3520	6436.3965	0.0000	1.0000
JC69G	12	-3203.8753	6489.7729	0.0000	1.0000
JC69I	12	-3203.9654	6489.9530	0.0000	1.0000
JC69IG	13	-3203.9765	6496.8103	0.0000	1.0000
JC69	11	-3217.3466	6509.8803	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/chr7_10487.phylip = 1 - 930;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10487.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/chr7_10487.phylip = 1 - 930;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10487.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 (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_10487.phylip = 1 - 930;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10487.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

