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

Input data from file "591-tomodel-phylip/chr6_3036.phylip" (7 taxa, 913 characters)

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-3130.2490	6294.4980	0.3359	0.3359
GTRG	20	-3127.4165	6294.8329	0.2841	0.6201
SYMIG	18	-3130.2254	6296.4509	0.1265	0.7466
GTRIG	21	-3127.3653	6296.7305	0.1100	0.8566
K2PG	13	-3136.3936	6298.7871	0.0393	0.8960
SYMI	17	-3132.6174	6299.2348	0.0315	0.9274
GTRI	20	-3129.8513	6299.7026	0.0249	0.9523
HKYG	16	-3134.1111	6300.2221	0.0192	0.9715
K2PIG	14	-3136.3798	6300.7596	0.0147	0.9862
HKYIG	17	-3134.0990	6302.1979	0.0071	0.9933
K2PI	13	-3138.5543	6303.1087	0.0045	0.9979
HKYI	16	-3136.3070	6304.6139	0.0021	1.0000
GTR	19	-3183.2721	6404.5442	0.0000	1.0000
JC69G	12	-3190.8919	6405.7838	0.0000	1.0000
SYM	16	-3187.2007	6406.4014	0.0000	1.0000
F81G	15	-3188.4579	6406.9158	0.0000	1.0000
JC69IG	13	-3190.8966	6407.7932	0.0000	1.0000
JC69I	12	-3192.4097	6408.8193	0.0000	1.0000
F81IG	16	-3188.4579	6408.9157	0.0000	1.0000
F81I	15	-3190.0091	6410.0182	0.0000	1.0000
HKY	15	-3190.4749	6410.9497	0.0000	1.0000
K2P	12	-3193.6141	6411.2282	0.0000	1.0000
JC69	11	-3245.8947	6513.7894	0.0000	1.0000
F81	14	-3243.1839	6514.3679	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3136.3936	6361.4047	0.8697	0.8697
K2PI	13	-3138.5543	6365.7262	0.1002	0.9699
K2PIG	14	-3136.3798	6368.1939	0.0292	0.9991
SYMG	17	-3130.2490	6376.3825	0.0005	0.9996
HKYG	16	-3134.1111	6377.2899	0.0003	0.9999
SYMI	17	-3132.6174	6381.1193	0.0000	0.9999
HKYI	16	-3136.3070	6381.6817	0.0000	1.0000
SYMIG	18	-3130.2254	6383.1521	0.0000	1.0000
HKYIG	17	-3134.0990	6384.0824	0.0000	1.0000
GTRG	20	-3127.4165	6391.1677	0.0000	1.0000
GTRI	20	-3129.8513	6396.0373	0.0000	1.0000
GTRIG	21	-3127.3653	6397.8820	0.0000	1.0000
JC69G	12	-3190.8919	6463.5846	0.0000	1.0000
JC69I	12	-3192.4097	6466.6202	0.0000	1.0000
K2P	12	-3193.6141	6469.0291	0.0000	1.0000
JC69IG	13	-3190.8966	6470.4108	0.0000	1.0000
F81G	15	-3188.4579	6479.1668	0.0000	1.0000
F81I	15	-3190.0091	6482.2693	0.0000	1.0000
HKY	15	-3190.4749	6483.2007	0.0000	1.0000
SYM	16	-3187.2007	6483.4692	0.0000	1.0000
F81IG	16	-3188.4579	6485.9835	0.0000	1.0000
GTR	19	-3183.2721	6496.0622	0.0000	1.0000
JC69	11	-3245.8947	6566.7735	0.0000	1.0000
F81	14	-3243.1839	6581.8022	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (SYMG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_3036.phylip = 1 - 913;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3036.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=Fixed(Equal);
END;


[Mrbayes block for the best AICc model (SYMG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_3036.phylip = 1 - 913;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3036.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=Fixed(Equal);
END;


[Mrbayes block for the best BIC model (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_3036.phylip = 1 - 913;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3036.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

