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

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

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-2557.3815	5155.8051	0.6792	0.6792
GTRIG	21	-2557.5961	5158.3400	0.1912	0.8704
GTRI	20	-2559.0383	5159.1189	0.1295	0.9999
SYMG	17	-2570.4878	5175.7321	0.0000	1.0000
SYMIG	18	-2570.6791	5178.2046	0.0000	1.0000
SYMI	17	-2572.0665	5178.8894	0.0000	1.0000
HKYG	16	-2573.9406	5180.5528	0.0000	1.0000
HKYIG	17	-2573.9475	5182.6515	0.0000	1.0000
HKYI	16	-2575.3984	5183.4684	0.0000	1.0000
GTR	19	-2575.0224	5188.9866	0.0000	1.0000
K2PG	13	-2582.3988	5191.2453	0.0000	1.0000
K2PI	13	-2583.8357	5194.1191	0.0000	1.0000
K2PIG	14	-2583.1285	5194.7743	0.0000	1.0000
SYM	16	-2587.8180	5208.3076	0.0000	1.0000
HKY	15	-2590.9003	5212.3924	0.0000	1.0000
K2P	12	-2599.5096	5223.4025	0.0000	1.0000
F81G	15	-2623.3711	5277.3342	0.0000	1.0000
F81I	15	-2624.3229	5279.2376	0.0000	1.0000
F81IG	16	-2623.3713	5279.4143	0.0000	1.0000
JC69G	12	-2630.1908	5284.7648	0.0000	1.0000
JC69I	12	-2631.1432	5286.6697	0.0000	1.0000
JC69IG	13	-2630.2030	5286.8537	0.0000	1.0000
F81	14	-2638.1630	5304.8432	0.0000	1.0000
JC69	11	-2645.1142	5312.5524	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2557.3815	5154.7630	0.6722	0.6722
GTRIG	21	-2557.5961	5157.1922	0.1995	0.8717
GTRI	20	-2559.0383	5158.0767	0.1282	1.0000
SYMG	17	-2570.4878	5174.9757	0.0000	1.0000
SYMIG	18	-2570.6791	5177.3581	0.0000	1.0000
SYMI	17	-2572.0665	5178.1329	0.0000	1.0000
HKYG	16	-2573.9406	5179.8812	0.0000	1.0000
HKYIG	17	-2573.9475	5181.8951	0.0000	1.0000
HKYI	16	-2575.3984	5182.7968	0.0000	1.0000
GTR	19	-2575.0224	5188.0448	0.0000	1.0000
K2PG	13	-2582.3988	5190.7976	0.0000	1.0000
K2PI	13	-2583.8357	5193.6714	0.0000	1.0000
K2PIG	14	-2583.1285	5194.2571	0.0000	1.0000
SYM	16	-2587.8180	5207.6360	0.0000	1.0000
HKY	15	-2590.9003	5211.8005	0.0000	1.0000
K2P	12	-2599.5096	5223.0192	0.0000	1.0000
F81G	15	-2623.3711	5276.7423	0.0000	1.0000
F81I	15	-2624.3229	5278.6457	0.0000	1.0000
F81IG	16	-2623.3713	5278.7427	0.0000	1.0000
JC69G	12	-2630.1908	5284.3815	0.0000	1.0000
JC69I	12	-2631.1432	5286.2864	0.0000	1.0000
JC69IG	13	-2630.2030	5286.4060	0.0000	1.0000
F81	14	-2638.1630	5304.3260	0.0000	1.0000
JC69	11	-2645.1142	5312.2285	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2557.3815	5249.1191	0.6201	0.6201
K2PG	13	-2582.3988	5252.1290	0.1377	0.7577
GTRI	20	-2559.0383	5252.4328	0.1183	0.8760
K2PI	13	-2583.8357	5255.0029	0.0327	0.9087
SYMG	17	-2570.4878	5255.1783	0.0300	0.9387
HKYG	16	-2573.9406	5255.3661	0.0273	0.9660
GTRIG	21	-2557.5961	5256.2661	0.0174	0.9834
HKYI	16	-2575.3984	5258.2816	0.0064	0.9897
SYMI	17	-2572.0665	5258.3356	0.0062	0.9959
K2PIG	14	-2583.1285	5260.3063	0.0023	0.9982
HKYIG	17	-2573.9475	5262.0977	0.0009	0.9991
SYMIG	18	-2570.6791	5262.2786	0.0009	1.0000
GTR	19	-2575.0224	5277.6831	0.0000	1.0000
K2P	12	-2599.5096	5279.6329	0.0000	1.0000
HKY	15	-2590.9003	5282.5676	0.0000	1.0000
SYM	16	-2587.8180	5283.1209	0.0000	1.0000
JC69G	12	-2630.1908	5340.9952	0.0000	1.0000
JC69I	12	-2631.1432	5342.9000	0.0000	1.0000
F81G	15	-2623.3711	5347.5094	0.0000	1.0000
JC69IG	13	-2630.2030	5347.7374	0.0000	1.0000
F81I	15	-2624.3229	5349.4128	0.0000	1.0000
F81IG	16	-2623.3713	5354.2276	0.0000	1.0000
JC69	11	-2645.1142	5364.1243	0.0000	1.0000
F81	14	-2638.1630	5370.3753	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/chr6_7982.phylip = 1 - 827;
 Partition Dummy = 1:683-tomodel-phylip/chr6_7982.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/chr6_7982.phylip = 1 - 827;
 Partition Dummy = 1:683-tomodel-phylip/chr6_7982.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 (GTRG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_7982.phylip = 1 - 827;
 Partition Dummy = 1:683-tomodel-phylip/chr6_7982.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;


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

