
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:58 02/17/2011

Input data from file "183-phylip/chr7_10687.phylip" (29 taxa, 623 characters)

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


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	64	-2725.7986	5594.5077	0.6905	0.6905
GTRIG	65	-2725.3803	5596.1645	0.3016	0.9921
HKYG	60	-2735.5858	5604.1966	0.0054	0.9975
HKYIG	61	-2735.1450	5605.7731	0.0025	1.0000
GTRI	64	-2737.1763	5617.2630	0.0000	1.0000
K2PG	57	-2749.5739	5624.8504	0.0000	1.0000
HKYI	60	-2747.1717	5627.3683	0.0000	1.0000
K2PIG	58	-2749.7357	5627.6062	0.0000	1.0000
SYMG	61	-2746.6240	5628.7310	0.0000	1.0000
SYMIG	62	-2745.7512	5629.4525	0.0000	1.0000
K2PI	57	-2763.1333	5651.9693	0.0000	1.0000
SYMI	61	-2758.6097	5652.7025	0.0000	1.0000
F81G	59	-2791.1354	5712.8463	0.0000	1.0000
F81IG	60	-2790.8306	5714.6860	0.0000	1.0000
GTR	63	-2791.9947	5724.4152	0.0000	1.0000
JC69G	56	-2801.4696	5726.2183	0.0000	1.0000
JC69IG	57	-2800.5749	5726.8524	0.0000	1.0000
HKY	59	-2799.6645	5729.9044	0.0000	1.0000
F81I	59	-2801.9755	5734.5265	0.0000	1.0000
JC69I	56	-2812.6403	5748.5597	0.0000	1.0000
K2P	56	-2814.9501	5753.1794	0.0000	1.0000
SYM	60	-2811.0197	5755.0642	0.0000	1.0000
F81	58	-2851.8856	5831.9060	0.0000	1.0000
JC69	55	-2862.6236	5846.1113	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2725.7986	5579.5973	0.6395	0.6395
GTRIG	65	-2725.3803	5580.7605	0.3574	0.9969
HKYG	60	-2735.5858	5591.1717	0.0020	0.9989
HKYIG	61	-2735.1450	5592.2900	0.0011	1.0000
GTRI	64	-2737.1763	5602.3526	0.0000	1.0000
K2PG	57	-2749.5739	5613.1478	0.0000	1.0000
HKYI	60	-2747.1717	5614.3434	0.0000	1.0000
SYMG	61	-2746.6240	5615.2479	0.0000	1.0000
K2PIG	58	-2749.7357	5615.4714	0.0000	1.0000
SYMIG	62	-2745.7512	5615.5025	0.0000	1.0000
SYMI	61	-2758.6097	5639.2194	0.0000	1.0000
K2PI	57	-2763.1333	5640.2667	0.0000	1.0000
F81G	59	-2791.1354	5700.2708	0.0000	1.0000
F81IG	60	-2790.8306	5701.6611	0.0000	1.0000
GTR	63	-2791.9947	5709.9895	0.0000	1.0000
JC69G	56	-2801.4696	5714.9392	0.0000	1.0000
JC69IG	57	-2800.5749	5715.1498	0.0000	1.0000
HKY	59	-2799.6645	5717.3289	0.0000	1.0000
F81I	59	-2801.9755	5721.9510	0.0000	1.0000
JC69I	56	-2812.6403	5737.2805	0.0000	1.0000
K2P	56	-2814.9501	5741.9003	0.0000	1.0000
SYM	60	-2811.0197	5742.0393	0.0000	1.0000
F81	58	-2851.8856	5819.7712	0.0000	1.0000
JC69	55	-2862.6236	5835.2471	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2735.5858	5857.2445	0.8893	0.8893
HKYIG	61	-2735.1450	5862.7973	0.0554	0.9447
GTRG	64	-2725.7986	5863.4083	0.0408	0.9855
K2PG	57	-2749.5739	5865.9169	0.0116	0.9971
GTRIG	65	-2725.3803	5869.0061	0.0025	0.9996
K2PIG	58	-2749.7357	5872.6751	0.0004	1.0000
HKYI	60	-2747.1717	5880.4162	0.0000	1.0000
SYMG	61	-2746.6240	5885.7553	0.0000	1.0000
GTRI	64	-2737.1763	5886.1636	0.0000	1.0000
SYMIG	62	-2745.7512	5890.4443	0.0000	1.0000
K2PI	57	-2763.1333	5893.0358	0.0000	1.0000
SYMI	61	-2758.6097	5909.7268	0.0000	1.0000
F81G	59	-2791.1354	5961.9090	0.0000	1.0000
JC69G	56	-2801.4696	5963.2738	0.0000	1.0000
F81IG	60	-2790.8306	5967.7339	0.0000	1.0000
JC69IG	57	-2800.5749	5967.9189	0.0000	1.0000
HKY	59	-2799.6645	5978.9672	0.0000	1.0000
F81I	59	-2801.9755	5983.5892	0.0000	1.0000
JC69I	56	-2812.6403	5985.6151	0.0000	1.0000
GTR	63	-2791.9947	5989.3659	0.0000	1.0000
K2P	56	-2814.9501	5990.2349	0.0000	1.0000
SYM	60	-2811.0197	6008.1121	0.0000	1.0000
F81	58	-2851.8856	6076.9749	0.0000	1.0000
JC69	55	-2862.6236	6079.1472	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 183-phylip/chr7_10687.phylip = 1 - 623;
 Partition Dummy = 1:183-phylip/chr7_10687.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 183-phylip/chr7_10687.phylip = 1 - 623;
 Partition Dummy = 1:183-phylip/chr7_10687.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 183-phylip/chr7_10687.phylip = 1 - 623;
 Partition Dummy = 1:183-phylip/chr7_10687.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

