
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:41 04/04/2011

Input data from file "591-tomodel-phylip/chr6_9803.phylip" (7 taxa, 660 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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.8456	3645.0058	0.4899	0.4899
GTRI	20	-1802.2893	3645.8931	0.3144	0.8043
GTRIG	21	-1801.7239	3646.8960	0.1904	0.9947
HKYG	16	-1811.9343	3656.7147	0.0014	0.9961
HKYI	16	-1812.2857	3657.4174	0.0010	0.9971
SYMG	17	-1811.2341	3657.4216	0.0010	0.9981
SYMI	17	-1811.5082	3657.9697	0.0007	0.9988
HKYIG	17	-1811.5617	3658.0766	0.0007	0.9995
SYMIG	18	-1811.1693	3659.4056	0.0004	0.9999
K2PG	13	-1818.3840	3663.3316	0.0001	0.9999
K2PI	13	-1818.6779	3663.9193	0.0000	1.0000
K2PIG	14	-1817.9517	3664.5546	0.0000	1.0000
F81G	15	-1833.1666	3697.0786	0.0000	1.0000
F81I	15	-1833.5340	3697.8133	0.0000	1.0000
F81IG	16	-1833.0132	3698.8725	0.0000	1.0000
JC69G	12	-1839.4412	3703.3646	0.0000	1.0000
JC69I	12	-1839.8326	3704.1475	0.0000	1.0000
JC69IG	13	-1839.0487	3704.6608	0.0000	1.0000
GTR	19	-1848.2048	3735.5972	0.0000	1.0000
HKY	15	-1858.4599	3747.6651	0.0000	1.0000
SYM	16	-1858.0532	3748.9524	0.0000	1.0000
K2P	12	-1864.4541	3753.3903	0.0000	1.0000
F81	14	-1877.1921	3783.0354	0.0000	1.0000
JC69	11	-1884.1313	3790.6699	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1801.8456	3643.6913	0.4840	0.4840
GTRI	20	-1802.2893	3644.5786	0.3106	0.7946
GTRIG	21	-1801.7239	3645.4477	0.2011	0.9957
HKYG	16	-1811.9343	3655.8687	0.0011	0.9968
SYMG	17	-1811.2341	3656.4683	0.0008	0.9976
HKYI	16	-1812.2857	3656.5713	0.0008	0.9984
SYMI	17	-1811.5082	3657.0164	0.0006	0.9990
HKYIG	17	-1811.5617	3657.1233	0.0006	0.9996
SYMIG	18	-1811.1693	3658.3385	0.0003	0.9999
K2PG	13	-1818.3840	3662.7681	0.0000	1.0000
K2PI	13	-1818.6779	3663.3559	0.0000	1.0000
K2PIG	14	-1817.9517	3663.9034	0.0000	1.0000
F81G	15	-1833.1666	3696.3333	0.0000	1.0000
F81I	15	-1833.5340	3697.0680	0.0000	1.0000
F81IG	16	-1833.0132	3698.0265	0.0000	1.0000
JC69G	12	-1839.4412	3702.8824	0.0000	1.0000
JC69I	12	-1839.8326	3703.6652	0.0000	1.0000
JC69IG	13	-1839.0487	3704.0974	0.0000	1.0000
GTR	19	-1848.2048	3734.4097	0.0000	1.0000
HKY	15	-1858.4599	3746.9197	0.0000	1.0000
SYM	16	-1858.0532	3748.1063	0.0000	1.0000
K2P	12	-1864.4541	3752.9081	0.0000	1.0000
F81	14	-1877.1921	3782.3842	0.0000	1.0000
JC69	11	-1884.1313	3790.2625	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1818.3840	3721.1672	0.5320	0.5320
K2PI	13	-1818.6779	3721.7550	0.3965	0.9285
K2PIG	14	-1817.9517	3726.7948	0.0319	0.9604
HKYG	16	-1811.9343	3727.7445	0.0198	0.9803
HKYI	16	-1812.2857	3728.4472	0.0140	0.9942
SYMG	17	-1811.2341	3732.8364	0.0016	0.9958
SYMI	17	-1811.5082	3733.3845	0.0012	0.9970
HKYIG	17	-1811.5617	3733.4914	0.0011	0.9981
GTRG	20	-1801.8456	3733.5361	0.0011	0.9992
GTRI	20	-1802.2893	3734.4234	0.0007	0.9999
SYMIG	18	-1811.1693	3739.1989	0.0001	1.0000
GTRIG	21	-1801.7239	3739.7848	0.0000	1.0000
JC69G	12	-1839.4412	3756.7893	0.0000	1.0000
JC69I	12	-1839.8326	3757.5721	0.0000	1.0000
JC69IG	13	-1839.0487	3762.4965	0.0000	1.0000
F81G	15	-1833.1666	3763.7169	0.0000	1.0000
F81I	15	-1833.5340	3764.4516	0.0000	1.0000
F81IG	16	-1833.0132	3769.9023	0.0000	1.0000
K2P	12	-1864.4541	3806.8150	0.0000	1.0000
HKY	15	-1858.4599	3814.3033	0.0000	1.0000
GTR	19	-1848.2048	3819.7622	0.0000	1.0000
SYM	16	-1858.0532	3819.9822	0.0000	1.0000
JC69	11	-1884.1313	3839.6772	0.0000	1.0000
F81	14	-1877.1921	3845.2756	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 591-tomodel-phylip/chr6_9803.phylip = 1 - 660;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9803.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 591-tomodel-phylip/chr6_9803.phylip = 1 - 660;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9803.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_9803.phylip = 1 - 660;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9803.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

