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

Input data from file "683-tomodel-phylip/chr2_23507.phylip" (7 taxa, 687 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	20	-1801.8573	3644.9759	0.2362	0.2362
HKYG	16	-1806.1804	3645.1727	0.2141	0.4502
GTRI	20	-1801.9954	3645.2522	0.2057	0.6560
HKYI	16	-1806.3371	3645.4860	0.1830	0.8390
GTRIG	21	-1801.8627	3647.1148	0.0811	0.9200
HKYIG	17	-1806.1900	3647.2949	0.0741	0.9941
HKY	15	-1811.2963	3653.3079	0.0037	0.9978
GTR	19	-1807.5882	3654.3159	0.0022	1.0000
K2PG	13	-1832.7024	3691.9457	0.0000	1.0000
K2PI	13	-1832.8653	3692.2714	0.0000	1.0000
K2PIG	14	-1832.7149	3694.0547	0.0000	1.0000
SYMG	17	-1831.7201	3698.3549	0.0000	1.0000
SYMI	17	-1831.8525	3698.6197	0.0000	1.0000
K2P	12	-1837.3875	3699.2378	0.0000	1.0000
SYMIG	18	-1831.7306	3700.4851	0.0000	1.0000
SYM	16	-1836.5533	3705.9186	0.0000	1.0000
F81I	15	-1846.0040	3722.7233	0.0000	1.0000
F81G	15	-1846.0475	3722.8103	0.0000	1.0000
F81IG	16	-1846.0406	3724.8931	0.0000	1.0000
F81	14	-1849.8171	3728.2592	0.0000	1.0000
JC69I	12	-1868.0631	3760.5890	0.0000	1.0000
JC69G	12	-1868.0776	3760.6180	0.0000	1.0000
JC69IG	13	-1868.0741	3762.6892	0.0000	1.0000
JC69	11	-1871.8434	3766.0779	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1801.8573	3643.7147	0.2588	0.2588
GTRI	20	-1801.9954	3643.9909	0.2254	0.4842
HKYG	16	-1806.1804	3644.3608	0.1873	0.6715
HKYI	16	-1806.3371	3644.6741	0.1602	0.8317
GTRIG	21	-1801.8627	3645.7254	0.0947	0.9264
HKYIG	17	-1806.1900	3646.3801	0.0683	0.9947
HKY	15	-1811.2963	3652.5925	0.0031	0.9977
GTR	19	-1807.5882	3653.1765	0.0023	1.0000
K2PG	13	-1832.7024	3691.4049	0.0000	1.0000
K2PI	13	-1832.8653	3691.7305	0.0000	1.0000
K2PIG	14	-1832.7149	3693.4297	0.0000	1.0000
SYMG	17	-1831.7201	3697.4401	0.0000	1.0000
SYMI	17	-1831.8525	3697.7049	0.0000	1.0000
K2P	12	-1837.3875	3698.7749	0.0000	1.0000
SYMIG	18	-1831.7306	3699.4612	0.0000	1.0000
SYM	16	-1836.5533	3705.1066	0.0000	1.0000
F81I	15	-1846.0040	3722.0080	0.0000	1.0000
F81G	15	-1846.0475	3722.0949	0.0000	1.0000
F81IG	16	-1846.0406	3724.0812	0.0000	1.0000
F81	14	-1849.8171	3727.6342	0.0000	1.0000
JC69I	12	-1868.0631	3760.1261	0.0000	1.0000
JC69G	12	-1868.0776	3760.1551	0.0000	1.0000
JC69IG	13	-1868.0741	3762.1483	0.0000	1.0000
JC69	11	-1871.8434	3765.6868	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1806.1804	3716.8781	0.4877	0.4877
HKYI	16	-1806.3371	3717.1914	0.4170	0.9047
HKY	15	-1811.2963	3720.5775	0.0767	0.9814
HKYIG	17	-1806.1900	3723.4298	0.0184	0.9998
GTRG	20	-1801.8573	3734.3614	0.0001	0.9999
GTRI	20	-1801.9954	3734.6376	0.0001	1.0000
GTR	19	-1807.5882	3739.2908	0.0000	1.0000
GTRIG	21	-1801.8627	3740.9044	0.0000	1.0000
K2PG	13	-1832.7024	3750.3252	0.0000	1.0000
K2PI	13	-1832.8653	3750.6508	0.0000	1.0000
K2P	12	-1837.3875	3753.1630	0.0000	1.0000
K2PIG	14	-1832.7149	3756.8824	0.0000	1.0000
SYMG	17	-1831.7201	3774.4898	0.0000	1.0000
SYMI	17	-1831.8525	3774.7546	0.0000	1.0000
SYM	16	-1836.5533	3777.6240	0.0000	1.0000
SYMIG	18	-1831.7306	3781.0432	0.0000	1.0000
F81I	15	-1846.0040	3789.9930	0.0000	1.0000
F81G	15	-1846.0475	3790.0800	0.0000	1.0000
F81	14	-1849.8171	3791.0869	0.0000	1.0000
F81IG	16	-1846.0406	3796.5985	0.0000	1.0000
JC69I	12	-1868.0631	3814.5141	0.0000	1.0000
JC69G	12	-1868.0776	3814.5431	0.0000	1.0000
JC69	11	-1871.8434	3815.5425	0.0000	1.0000
JC69IG	13	-1868.0741	3821.0686	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 683-tomodel-phylip/chr2_23507.phylip = 1 - 687;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23507.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 683-tomodel-phylip/chr2_23507.phylip = 1 - 687;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23507.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 683-tomodel-phylip/chr2_23507.phylip = 1 - 687;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23507.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

