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

Input data from file "591-tomodel-phylip/chr3_764.phylip" (7 taxa, 831 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	-1821.0185	3683.0741	0.2763	0.2763
HKYG	16	-1825.3388	3683.3460	0.2412	0.5175
GTRI	20	-1821.6935	3684.4241	0.1407	0.6582
HKYI	16	-1825.9912	3684.6507	0.1256	0.7838
GTRIG	21	-1820.8547	3684.8515	0.1136	0.8974
HKYIG	17	-1825.1576	3685.0680	0.1020	0.9993
GTR	19	-1828.6102	3696.1576	0.0004	0.9997
HKY	15	-1833.1517	3696.8924	0.0003	1.0000
F81G	15	-1840.1424	3710.8737	0.0000	1.0000
F81I	15	-1840.7931	3712.1752	0.0000	1.0000
F81IG	16	-1840.0142	3712.6967	0.0000	1.0000
F81	14	-1847.2137	3722.9422	0.0000	1.0000
K2PG	13	-1872.1531	3770.7517	0.0000	1.0000
K2PI	13	-1872.8162	3772.0779	0.0000	1.0000
K2PIG	14	-1872.0013	3772.5172	0.0000	1.0000
SYMG	17	-1868.9571	3772.6670	0.0000	1.0000
SYMI	17	-1869.6860	3774.1248	0.0000	1.0000
SYMIG	18	-1868.8205	3774.4834	0.0000	1.0000
K2P	12	-1879.4723	3783.3260	0.0000	1.0000
SYM	16	-1875.8187	3784.3058	0.0000	1.0000
JC69G	12	-1883.7042	3791.7898	0.0000	1.0000
JC69I	12	-1884.3475	3793.0764	0.0000	1.0000
JC69IG	13	-1883.5840	3793.6135	0.0000	1.0000
JC69	11	-1890.6576	3803.6375	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1821.0185	3682.0370	0.2968	0.2968
HKYG	16	-1825.3388	3682.6777	0.2155	0.5123
GTRI	20	-1821.6935	3683.3871	0.1511	0.6635
GTRIG	21	-1820.8547	3683.7093	0.1286	0.7921
HKYI	16	-1825.9912	3683.9824	0.1122	0.9043
HKYIG	17	-1825.1576	3684.3153	0.0950	0.9994
GTR	19	-1828.6102	3695.2205	0.0004	0.9998
HKY	15	-1833.1517	3696.3035	0.0002	1.0000
F81G	15	-1840.1424	3710.2848	0.0000	1.0000
F81I	15	-1840.7931	3711.5863	0.0000	1.0000
F81IG	16	-1840.0142	3712.0284	0.0000	1.0000
F81	14	-1847.2137	3722.4275	0.0000	1.0000
K2PG	13	-1872.1531	3770.3061	0.0000	1.0000
K2PI	13	-1872.8162	3771.6323	0.0000	1.0000
SYMG	17	-1868.9571	3771.9142	0.0000	1.0000
K2PIG	14	-1872.0013	3772.0025	0.0000	1.0000
SYMI	17	-1869.6860	3773.3720	0.0000	1.0000
SYMIG	18	-1868.8205	3773.6410	0.0000	1.0000
K2P	12	-1879.4723	3782.9446	0.0000	1.0000
SYM	16	-1875.8187	3783.6375	0.0000	1.0000
JC69G	12	-1883.7042	3791.4084	0.0000	1.0000
JC69I	12	-1884.3475	3792.6950	0.0000	1.0000
JC69IG	13	-1883.5840	3793.1679	0.0000	1.0000
JC69	11	-1890.6576	3803.3151	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1825.3388	3758.2398	0.6352	0.6352
HKYI	16	-1825.9912	3759.5445	0.3308	0.9661
HKYIG	17	-1825.1576	3764.6000	0.0264	0.9925
HKY	15	-1833.1517	3767.1429	0.0074	0.9999
GTRG	20	-1821.0185	3776.4896	0.0001	1.0000
GTRI	20	-1821.6935	3777.8397	0.0000	1.0000
F81G	15	-1840.1424	3781.1242	0.0000	1.0000
F81I	15	-1840.7931	3782.4257	0.0000	1.0000
GTRIG	21	-1820.8547	3782.8845	0.0000	1.0000
GTR	19	-1828.6102	3784.9504	0.0000	1.0000
F81IG	16	-1840.0142	3787.5905	0.0000	1.0000
F81	14	-1847.2137	3788.5443	0.0000	1.0000
K2PG	13	-1872.1531	3831.7003	0.0000	1.0000
K2PI	13	-1872.8162	3833.0265	0.0000	1.0000
K2PIG	14	-1872.0013	3838.1193	0.0000	1.0000
K2P	12	-1879.4723	3839.6161	0.0000	1.0000
JC69G	12	-1883.7042	3848.0800	0.0000	1.0000
JC69I	12	-1884.3475	3849.3666	0.0000	1.0000
SYMG	17	-1868.9571	3852.1989	0.0000	1.0000
SYMI	17	-1869.6860	3853.6567	0.0000	1.0000
JC69IG	13	-1883.5840	3854.5621	0.0000	1.0000
JC69	11	-1890.6576	3855.2641	0.0000	1.0000
SYMIG	18	-1868.8205	3858.6484	0.0000	1.0000
SYM	16	-1875.8187	3859.1996	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/chr3_764.phylip = 1 - 831;
 Partition Dummy = 1:591-tomodel-phylip/chr3_764.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/chr3_764.phylip = 1 - 831;
 Partition Dummy = 1:591-tomodel-phylip/chr3_764.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 591-tomodel-phylip/chr3_764.phylip = 1 - 831;
 Partition Dummy = 1:591-tomodel-phylip/chr3_764.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

