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

Input data from file "183-phylip/chr8_675.phylip" (29 taxa, 760 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	64	-5229.2398	10598.4507	0.6097	0.6097
GTRIG	65	-5228.4900	10599.3432	0.3903	1.0000
HKYG	60	-5249.0087	10628.4896	0.0000	1.0000
HKYIG	61	-5248.4267	10629.6901	0.0000	1.0000
GTRI	64	-5268.5059	10676.9830	0.0000	1.0000
HKYI	60	-5283.8482	10698.1686	0.0000	1.0000
SYMG	61	-5303.6775	10740.1917	0.0000	1.0000
SYMIG	62	-5302.8832	10740.9744	0.0000	1.0000
K2PG	57	-5310.9405	10745.2999	0.0000	1.0000
K2PIG	58	-5310.1734	10746.1100	0.0000	1.0000
SYMI	61	-5335.6922	10804.2211	0.0000	1.0000
F81G	59	-5338.2359	10804.5860	0.0000	1.0000
F81IG	60	-5337.4460	10805.3641	0.0000	1.0000
K2PI	57	-5342.6571	10808.7330	0.0000	1.0000
F81I	59	-5368.0351	10864.1845	0.0000	1.0000
JC69G	56	-5376.7538	10874.5887	0.0000	1.0000
JC69IG	57	-5375.7253	10874.8695	0.0000	1.0000
GTR	63	-5389.0022	10915.5907	0.0000	1.0000
JC69I	56	-5405.2178	10931.5167	0.0000	1.0000
HKY	59	-5403.0814	10934.2771	0.0000	1.0000
SYM	60	-5452.6694	11035.8110	0.0000	1.0000
K2P	56	-5459.2918	11039.6648	0.0000	1.0000
F81	58	-5483.2600	11092.2832	0.0000	1.0000
JC69	55	-5519.5926	11157.9352	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-5229.2398	10586.4795	0.5622	0.5622
GTRIG	65	-5228.4900	10586.9801	0.4378	1.0000
HKYG	60	-5249.0087	10618.0175	0.0000	1.0000
HKYIG	61	-5248.4267	10618.8534	0.0000	1.0000
GTRI	64	-5268.5059	10665.0118	0.0000	1.0000
HKYI	60	-5283.8482	10687.6965	0.0000	1.0000
SYMG	61	-5303.6775	10729.3550	0.0000	1.0000
SYMIG	62	-5302.8832	10729.7663	0.0000	1.0000
K2PG	57	-5310.9405	10735.8811	0.0000	1.0000
K2PIG	58	-5310.1734	10736.3468	0.0000	1.0000
SYMI	61	-5335.6922	10793.3844	0.0000	1.0000
F81G	59	-5338.2359	10794.4717	0.0000	1.0000
F81IG	60	-5337.4460	10794.8920	0.0000	1.0000
K2PI	57	-5342.6571	10799.3142	0.0000	1.0000
F81I	59	-5368.0351	10854.0702	0.0000	1.0000
JC69IG	57	-5375.7253	10865.4507	0.0000	1.0000
JC69G	56	-5376.7538	10865.5076	0.0000	1.0000
GTR	63	-5389.0022	10904.0045	0.0000	1.0000
JC69I	56	-5405.2178	10922.4357	0.0000	1.0000
HKY	59	-5403.0814	10924.1628	0.0000	1.0000
SYM	60	-5452.6694	11025.3389	0.0000	1.0000
K2P	56	-5459.2918	11030.5837	0.0000	1.0000
F81	58	-5483.2600	11082.5200	0.0000	1.0000
JC69	55	-5519.5926	11149.1852	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-5229.2398	10883.0119	0.9273	0.9273
GTRIG	65	-5228.4900	10888.1458	0.0712	0.9985
HKYG	60	-5249.0087	10896.0166	0.0014	0.9999
HKYIG	61	-5248.4267	10901.4859	0.0001	1.0000
GTRI	64	-5268.5059	10961.5442	0.0000	1.0000
HKYI	60	-5283.8482	10965.6956	0.0000	1.0000
K2PG	57	-5310.9405	10999.9802	0.0000	1.0000
K2PIG	58	-5310.1734	11005.0793	0.0000	1.0000
SYMG	61	-5303.6775	11011.9875	0.0000	1.0000
SYMIG	62	-5302.8832	11017.0321	0.0000	1.0000
K2PI	57	-5342.6571	11063.4134	0.0000	1.0000
F81G	59	-5338.2359	11067.8375	0.0000	1.0000
F81IG	60	-5337.4460	11072.8911	0.0000	1.0000
SYMI	61	-5335.6922	11076.0169	0.0000	1.0000
JC69G	56	-5376.7538	11124.9734	0.0000	1.0000
F81I	59	-5368.0351	11127.4360	0.0000	1.0000
JC69IG	57	-5375.7253	11129.5498	0.0000	1.0000
JC69I	56	-5405.2178	11181.9015	0.0000	1.0000
GTR	63	-5389.0022	11195.9035	0.0000	1.0000
HKY	59	-5403.0814	11197.5286	0.0000	1.0000
K2P	56	-5459.2918	11290.0495	0.0000	1.0000
SYM	60	-5452.6694	11303.3380	0.0000	1.0000
F81	58	-5483.2600	11351.2525	0.0000	1.0000
JC69	55	-5519.5926	11404.0177	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/chr8_675.phylip = 1 - 760;
 Partition Dummy = 1:183-phylip/chr8_675.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/chr8_675.phylip = 1 - 760;
 Partition Dummy = 1:183-phylip/chr8_675.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 183-phylip/chr8_675.phylip = 1 - 760;
 Partition Dummy = 1:183-phylip/chr8_675.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

