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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: SYMG


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
GTRI	20	-1763.5998	3568.7668	0.3813	0.3813
GTRG	20	-1763.7642	3569.0956	0.3235	0.7048
GTRIG	21	-1763.5852	3570.8976	0.1314	0.8362
SYMG	17	-1768.4958	3572.1269	0.0711	0.9073
SYMI	17	-1768.5562	3572.2479	0.0669	0.9742
SYMIG	18	-1768.4411	3574.1537	0.0258	1.0000
GTR	19	-1776.6785	3592.7722	0.0000	1.0000
SYM	16	-1780.0412	3593.0898	0.0000	1.0000
K2PG	13	-1785.7068	3598.0839	0.0000	1.0000
K2PI	13	-1786.0835	3598.8374	0.0000	1.0000
HKYG	16	-1782.9645	3598.9365	0.0000	1.0000
HKYI	16	-1783.3420	3599.6914	0.0000	1.0000
K2PIG	14	-1785.7117	3600.1982	0.0000	1.0000
HKYIG	17	-1782.9734	3601.0822	0.0000	1.0000
K2P	12	-1798.2164	3621.0063	0.0000	1.0000
HKY	15	-1796.0876	3623.0625	0.0000	1.0000
JC69G	12	-1827.2731	3679.1198	0.0000	1.0000
JC69I	12	-1827.5686	3679.7108	0.0000	1.0000
JC69IG	13	-1827.2824	3681.2352	0.0000	1.0000
F81G	15	-1825.7499	3682.3870	0.0000	1.0000
F81I	15	-1826.0706	3683.0285	0.0000	1.0000
F81IG	16	-1825.7563	3684.5201	0.0000	1.0000
JC69	11	-1836.8090	3696.1023	0.0000	1.0000
F81	14	-1835.4847	3699.7443	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1763.5998	3567.1997	0.3889	0.3889
GTRG	20	-1763.7642	3567.5285	0.3299	0.7188
GTRIG	21	-1763.5852	3569.1705	0.1452	0.8639
SYMG	17	-1768.4958	3570.9915	0.0584	0.9223
SYMI	17	-1768.5562	3571.1125	0.0550	0.9773
SYMIG	18	-1768.4411	3572.8823	0.0227	1.0000
GTR	19	-1776.6785	3591.3570	0.0000	1.0000
SYM	16	-1780.0412	3592.0824	0.0000	1.0000
K2PG	13	-1785.7068	3597.4136	0.0000	1.0000
HKYG	16	-1782.9645	3597.9291	0.0000	1.0000
K2PI	13	-1786.0835	3598.1670	0.0000	1.0000
HKYI	16	-1783.3420	3598.6840	0.0000	1.0000
K2PIG	14	-1785.7117	3599.4233	0.0000	1.0000
HKYIG	17	-1782.9734	3599.9468	0.0000	1.0000
K2P	12	-1798.2164	3620.4328	0.0000	1.0000
HKY	15	-1796.0876	3622.1753	0.0000	1.0000
JC69G	12	-1827.2731	3678.5463	0.0000	1.0000
JC69I	12	-1827.5686	3679.1372	0.0000	1.0000
JC69IG	13	-1827.2824	3680.5649	0.0000	1.0000
F81G	15	-1825.7499	3681.4998	0.0000	1.0000
F81I	15	-1826.0706	3682.1412	0.0000	1.0000
F81IG	16	-1825.7563	3683.5127	0.0000	1.0000
JC69	11	-1836.8090	3695.6179	0.0000	1.0000
F81	14	-1835.4847	3698.9694	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	17	-1768.4958	3644.4751	0.4942	0.4942
SYMI	17	-1768.5562	3644.5961	0.4652	0.9594
SYMIG	18	-1768.4411	3650.6885	0.0221	0.9815
K2PG	13	-1785.7068	3653.6069	0.0051	0.9866
GTRI	20	-1763.5998	3653.6510	0.0050	0.9917
GTRG	20	-1763.7642	3653.9798	0.0043	0.9959
K2PI	13	-1786.0835	3654.3604	0.0035	0.9994
K2PIG	14	-1785.7117	3659.9393	0.0002	0.9997
GTRIG	21	-1763.5852	3659.9444	0.0002	0.9999
SYM	16	-1780.0412	3661.2434	0.0001	1.0000
HKYG	16	-1782.9645	3667.0901	0.0000	1.0000
HKYI	16	-1783.3420	3667.8451	0.0000	1.0000
K2P	12	-1798.2164	3672.3035	0.0000	1.0000
HKYIG	17	-1782.9734	3673.4304	0.0000	1.0000
GTR	19	-1776.6785	3673.4857	0.0000	1.0000
HKY	15	-1796.0876	3687.0138	0.0000	1.0000
JC69G	12	-1827.2731	3730.4171	0.0000	1.0000
JC69I	12	-1827.5686	3731.0080	0.0000	1.0000
JC69IG	13	-1827.2824	3736.7582	0.0000	1.0000
JC69	11	-1836.8090	3743.1661	0.0000	1.0000
F81G	15	-1825.7499	3746.3383	0.0000	1.0000
F81I	15	-1826.0706	3746.9797	0.0000	1.0000
F81IG	16	-1825.7563	3752.6737	0.0000	1.0000
F81	14	-1835.4847	3759.4854	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr11_3549.phylip = 1 - 557;
 Partition Dummy = 1:683-tomodel-phylip/chr11_3549.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr11_3549.phylip = 1 - 557;
 Partition Dummy = 1:683-tomodel-phylip/chr11_3549.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

