
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:56 04/30/2011

Input data from file "591-tomodel-temp-phylip/chr10_4435.phylip" (7 taxa, 694 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1759.2929	3551.3894	0.5325	0.5325
HKYIG	17	-1759.1128	3553.1309	0.2229	0.7554
HKYI	16	-1760.5168	3553.8372	0.1566	0.9120
GTRG	20	-1757.4104	3556.0689	0.0513	0.9633
GTRIG	21	-1757.2054	3557.7857	0.0217	0.9851
GTRI	20	-1758.6481	3558.5444	0.0149	1.0000
HKY	15	-1770.1337	3570.9753	0.0000	1.0000
GTR	19	-1768.7250	3576.5777	0.0000	1.0000
K2PG	13	-1786.9990	3600.5333	0.0000	1.0000
K2PIG	14	-1786.7188	3602.0562	0.0000	1.0000
K2PI	13	-1788.2338	3603.0029	0.0000	1.0000
F81G	15	-1786.8475	3604.4030	0.0000	1.0000
F81IG	16	-1786.7232	3606.2500	0.0000	1.0000
SYMG	17	-1785.8211	3606.5475	0.0000	1.0000
F81I	15	-1788.1108	3606.9295	0.0000	1.0000
SYMIG	18	-1785.4967	3608.0068	0.0000	1.0000
SYMI	17	-1787.1051	3609.1155	0.0000	1.0000
K2P	12	-1799.0408	3622.5398	0.0000	1.0000
F81	14	-1799.7767	3628.1719	0.0000	1.0000
SYM	16	-1798.3294	3629.4624	0.0000	1.0000
JC69G	12	-1812.1554	3648.7690	0.0000	1.0000
JC69IG	13	-1811.9901	3650.5156	0.0000	1.0000
JC69I	12	-1816.1817	3656.8216	0.0000	1.0000
JC69	11	-1825.5737	3673.5344	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1759.2929	3550.5858	0.5143	0.5143
HKYIG	17	-1759.1128	3552.2256	0.2266	0.7409
HKYI	16	-1760.5168	3553.0336	0.1513	0.8922
GTRG	20	-1757.4104	3554.8208	0.0619	0.9541
GTRIG	21	-1757.2054	3556.4107	0.0280	0.9820
GTRI	20	-1758.6481	3557.2963	0.0180	1.0000
HKY	15	-1770.1337	3570.2673	0.0000	1.0000
GTR	19	-1768.7250	3575.4501	0.0000	1.0000
K2PG	13	-1786.9990	3599.9980	0.0000	1.0000
K2PIG	14	-1786.7188	3601.4377	0.0000	1.0000
K2PI	13	-1788.2338	3602.4676	0.0000	1.0000
F81G	15	-1786.8475	3603.6951	0.0000	1.0000
F81IG	16	-1786.7232	3605.4465	0.0000	1.0000
SYMG	17	-1785.8211	3605.6422	0.0000	1.0000
F81I	15	-1788.1108	3606.2215	0.0000	1.0000
SYMIG	18	-1785.4967	3606.9934	0.0000	1.0000
SYMI	17	-1787.1051	3608.2102	0.0000	1.0000
K2P	12	-1799.0408	3622.0816	0.0000	1.0000
F81	14	-1799.7767	3627.5533	0.0000	1.0000
SYM	16	-1798.3294	3628.6589	0.0000	1.0000
JC69G	12	-1812.1554	3648.3109	0.0000	1.0000
JC69IG	13	-1811.9901	3649.9803	0.0000	1.0000
JC69I	12	-1816.1817	3656.3635	0.0000	1.0000
JC69	11	-1825.5737	3673.1473	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1759.2929	3623.2654	0.7462	0.7462
HKYI	16	-1760.5168	3625.7132	0.2195	0.9657
HKYIG	17	-1759.1128	3629.4476	0.0339	0.9996
HKY	15	-1770.1337	3638.4044	0.0004	1.0000
GTRG	20	-1757.4104	3645.6702	0.0000	1.0000
GTRI	20	-1758.6481	3648.1457	0.0000	1.0000
GTRIG	21	-1757.2054	3651.8026	0.0000	1.0000
K2PG	13	-1786.9990	3659.0502	0.0000	1.0000
K2PI	13	-1788.2338	3661.5197	0.0000	1.0000
GTR	19	-1768.7250	3661.7570	0.0000	1.0000
K2PIG	14	-1786.7188	3665.0323	0.0000	1.0000
F81G	15	-1786.8475	3671.8322	0.0000	1.0000
F81I	15	-1788.1108	3674.3586	0.0000	1.0000
K2P	12	-1799.0408	3676.5913	0.0000	1.0000
F81IG	16	-1786.7232	3678.1260	0.0000	1.0000
SYMG	17	-1785.8211	3682.8642	0.0000	1.0000
SYMI	17	-1787.1051	3685.4322	0.0000	1.0000
SYMIG	18	-1785.4967	3688.7579	0.0000	1.0000
F81	14	-1799.7767	3691.1479	0.0000	1.0000
SYM	16	-1798.3294	3701.3384	0.0000	1.0000
JC69G	12	-1812.1554	3702.8205	0.0000	1.0000
JC69IG	13	-1811.9901	3709.0324	0.0000	1.0000
JC69I	12	-1816.1817	3710.8732	0.0000	1.0000
JC69	11	-1825.5737	3723.1145	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-temp-phylip/chr10_4435.phylip = 1 - 694;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4435.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-temp-phylip/chr10_4435.phylip = 1 - 694;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4435.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-temp-phylip/chr10_4435.phylip = 1 - 694;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4435.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

