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

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

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


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
SYMG	17	-2614.0679	5263.0687	0.6394	0.6394
SYMIG	18	-2614.0678	5265.1798	0.2225	0.8619
GTRG	20	-2613.3885	5268.0633	0.0526	0.9146
K2PG	13	-2620.9583	5268.4682	0.0430	0.9575
GTRIG	21	-2613.3883	5270.1938	0.0181	0.9757
K2PIG	14	-2620.9584	5270.5541	0.0151	0.9908
HKYG	16	-2619.7623	5272.3526	0.0062	0.9970
HKYIG	17	-2619.7624	5274.4577	0.0022	0.9991
SYMI	17	-2620.8267	5276.5862	0.0007	0.9999
GTRI	20	-2620.1588	5281.6041	0.0001	1.0000
K2PI	13	-2627.8440	5282.2396	0.0000	1.0000
HKYI	16	-2626.7335	5286.2950	0.0000	1.0000
GTR	19	-2651.8562	5342.8745	0.0000	1.0000
K2P	12	-2659.9823	5344.4366	0.0000	1.0000
SYM	16	-2655.8559	5344.5399	0.0000	1.0000
HKY	15	-2657.6234	5345.9763	0.0000	1.0000
JC69G	12	-2661.2685	5347.0090	0.0000	1.0000
JC69IG	13	-2661.2684	5349.0884	0.0000	1.0000
F81G	15	-2659.3108	5349.3510	0.0000	1.0000
F81IG	16	-2659.3108	5351.4495	0.0000	1.0000
JC69I	12	-2667.3454	5359.1627	0.0000	1.0000
F81I	15	-2665.4411	5361.6117	0.0000	1.0000
JC69	11	-2697.9852	5418.3691	0.0000	1.0000
F81	14	-2695.2678	5419.1728	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-2614.0679	5262.1357	0.6281	0.6281
SYMIG	18	-2614.0678	5264.1355	0.2311	0.8592
GTRG	20	-2613.3885	5266.7770	0.0617	0.9209
K2PG	13	-2620.9583	5267.9167	0.0349	0.9558
GTRIG	21	-2613.3883	5268.7766	0.0227	0.9785
K2PIG	14	-2620.9584	5269.9168	0.0128	0.9913
HKYG	16	-2619.7623	5271.5246	0.0057	0.9970
HKYIG	17	-2619.7624	5273.5248	0.0021	0.9992
SYMI	17	-2620.8267	5275.6533	0.0007	0.9999
GTRI	20	-2620.1588	5280.3177	0.0001	1.0000
K2PI	13	-2627.8440	5281.6880	0.0000	1.0000
HKYI	16	-2626.7335	5285.4670	0.0000	1.0000
GTR	19	-2651.8562	5341.7124	0.0000	1.0000
SYM	16	-2655.8559	5343.7119	0.0000	1.0000
K2P	12	-2659.9823	5343.9646	0.0000	1.0000
HKY	15	-2657.6234	5345.2469	0.0000	1.0000
JC69G	12	-2661.2685	5346.5370	0.0000	1.0000
JC69IG	13	-2661.2684	5348.5369	0.0000	1.0000
F81G	15	-2659.3108	5348.6216	0.0000	1.0000
F81IG	16	-2659.3108	5350.6215	0.0000	1.0000
JC69I	12	-2667.3454	5358.6907	0.0000	1.0000
F81I	15	-2665.4411	5360.8822	0.0000	1.0000
JC69	11	-2697.9852	5417.9703	0.0000	1.0000
F81	14	-2695.2678	5418.5355	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2620.9583	5326.5887	0.9597	0.9597
K2PIG	14	-2620.9584	5333.1020	0.0370	0.9967
SYMG	17	-2614.0679	5338.8607	0.0021	0.9987
K2PI	13	-2627.8440	5340.3600	0.0010	0.9997
HKYG	16	-2619.7623	5343.7363	0.0002	0.9999
SYMIG	18	-2614.0678	5345.3737	0.0001	1.0000
HKYIG	17	-2619.7624	5350.2497	0.0000	1.0000
SYMI	17	-2620.8267	5352.3782	0.0000	1.0000
GTRG	20	-2613.3885	5357.0416	0.0000	1.0000
HKYI	16	-2626.7335	5357.6787	0.0000	1.0000
GTRIG	21	-2613.3883	5363.5545	0.0000	1.0000
GTRI	20	-2620.1588	5370.5823	0.0000	1.0000
K2P	12	-2659.9823	5398.1234	0.0000	1.0000
JC69G	12	-2661.2685	5400.6957	0.0000	1.0000
JC69IG	13	-2661.2684	5407.2089	0.0000	1.0000
JC69I	12	-2667.3454	5412.8495	0.0000	1.0000
HKY	15	-2657.6234	5412.9453	0.0000	1.0000
SYM	16	-2655.8559	5415.9236	0.0000	1.0000
F81G	15	-2659.3108	5416.3200	0.0000	1.0000
F81IG	16	-2659.3108	5422.8332	0.0000	1.0000
GTR	19	-2651.8562	5427.4638	0.0000	1.0000
F81I	15	-2665.4411	5428.5807	0.0000	1.0000
JC69	11	-2697.9852	5467.6159	0.0000	1.0000
F81	14	-2695.2678	5481.7207	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/chr7_1379.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr7_1379.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/chr7_1379.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr7_1379.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/chr7_1379.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr7_1379.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

