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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-1473.7512	2980.3927	0.2660	0.2660
HKYG	16	-1473.8794	2980.6491	0.2340	0.5001
GTRI	20	-1469.8868	2981.1575	0.1815	0.6816
GTRG	20	-1469.9123	2981.2085	0.1769	0.8585
HKYIG	17	-1473.8794	2982.7621	0.0814	0.9399
GTRIG	21	-1469.9211	2983.3669	0.0601	1.0000
HKY	15	-1488.2002	3007.1848	0.0000	1.0000
GTR	19	-1484.4089	3008.0679	0.0000	1.0000
F81I	15	-1502.4626	3035.7095	0.0000	1.0000
F81G	15	-1502.5976	3035.9796	0.0000	1.0000
F81IG	16	-1502.5977	3038.0857	0.0000	1.0000
K2PI	13	-1512.5523	3051.6975	0.0000	1.0000
K2PG	13	-1512.6199	3051.8327	0.0000	1.0000
K2PIG	14	-1512.6199	3053.9250	0.0000	1.0000
SYMI	17	-1511.0496	3057.1025	0.0000	1.0000
SYMG	17	-1511.0975	3057.1982	0.0000	1.0000
SYMIG	18	-1511.0975	3059.3181	0.0000	1.0000
F81	14	-1516.4628	3061.6107	0.0000	1.0000
K2P	12	-1528.1805	3080.8682	0.0000	1.0000
SYM	16	-1526.3823	3085.6550	0.0000	1.0000
JC69I	12	-1536.9494	3098.4061	0.0000	1.0000
JC69G	12	-1537.0259	3098.5592	0.0000	1.0000
JC69IG	13	-1537.0259	3100.6447	0.0000	1.0000
JC69	11	-1551.6851	3125.7988	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1473.7512	2979.5023	0.2360	0.2360
HKYG	16	-1473.8794	2979.7588	0.2076	0.4435
GTRI	20	-1469.8868	2979.7736	0.2060	0.6496
GTRG	20	-1469.9123	2979.8246	0.2008	0.8504
HKYIG	17	-1473.8794	2981.7588	0.0764	0.9268
GTRIG	21	-1469.9211	2981.8421	0.0732	1.0000
HKY	15	-1488.2002	3006.4005	0.0000	1.0000
GTR	19	-1484.4089	3006.8179	0.0000	1.0000
F81I	15	-1502.4626	3034.9252	0.0000	1.0000
F81G	15	-1502.5976	3035.1953	0.0000	1.0000
F81IG	16	-1502.5977	3037.1953	0.0000	1.0000
K2PI	13	-1512.5523	3051.1047	0.0000	1.0000
K2PG	13	-1512.6199	3051.2398	0.0000	1.0000
K2PIG	14	-1512.6199	3053.2398	0.0000	1.0000
SYMI	17	-1511.0496	3056.0992	0.0000	1.0000
SYMG	17	-1511.0975	3056.1949	0.0000	1.0000
SYMIG	18	-1511.0975	3058.1949	0.0000	1.0000
F81	14	-1516.4628	3060.9256	0.0000	1.0000
K2P	12	-1528.1805	3080.3609	0.0000	1.0000
SYM	16	-1526.3823	3084.7647	0.0000	1.0000
JC69I	12	-1536.9494	3097.8988	0.0000	1.0000
JC69G	12	-1537.0259	3098.0519	0.0000	1.0000
JC69IG	13	-1537.0259	3100.0519	0.0000	1.0000
JC69	11	-1551.6851	3125.3702	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1473.7512	3050.5830	0.5222	0.5222
HKYG	16	-1473.8794	3050.8394	0.4593	0.9815
HKYIG	17	-1473.8794	3057.2820	0.0183	0.9999
GTRI	20	-1469.8868	3068.6244	0.0001	0.9999
GTRG	20	-1469.9123	3068.6754	0.0001	1.0000
HKY	15	-1488.2002	3073.0386	0.0000	1.0000
GTRIG	21	-1469.9211	3075.1355	0.0000	1.0000
GTR	19	-1484.4089	3091.2261	0.0000	1.0000
F81I	15	-1502.4626	3101.5633	0.0000	1.0000
F81G	15	-1502.5976	3101.8334	0.0000	1.0000
F81IG	16	-1502.5977	3108.2760	0.0000	1.0000
K2PI	13	-1512.5523	3108.8577	0.0000	1.0000
K2PG	13	-1512.6199	3108.9928	0.0000	1.0000
K2PIG	14	-1512.6199	3115.4354	0.0000	1.0000
F81	14	-1516.4628	3123.1211	0.0000	1.0000
SYMI	17	-1511.0496	3131.6224	0.0000	1.0000
SYMG	17	-1511.0975	3131.7181	0.0000	1.0000
K2P	12	-1528.1805	3133.6714	0.0000	1.0000
SYMIG	18	-1511.0975	3138.1606	0.0000	1.0000
JC69I	12	-1536.9494	3151.2092	0.0000	1.0000
JC69G	12	-1537.0259	3151.3623	0.0000	1.0000
SYM	16	-1526.3823	3155.8453	0.0000	1.0000
JC69IG	13	-1537.0259	3157.8049	0.0000	1.0000
JC69	11	-1551.6851	3174.2381	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_299.phylip = 1 - 628;
 Partition Dummy = 1:591-tomodel-phylip/chr3_299.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_299.phylip = 1 - 628;
 Partition Dummy = 1:591-tomodel-phylip/chr3_299.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_299.phylip = 1 - 628;
 Partition Dummy = 1:591-tomodel-phylip/chr3_299.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

