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

Input data from file "591-tomodel-phylip/chr11_3420.phylip" (7 taxa, 460 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	-1497.7818	3037.4770	0.4015	0.4015
GTRIG	21	-1497.7021	3039.5137	0.1450	0.5465
HKYG	16	-1503.2294	3039.6867	0.1330	0.6795
K2PG	13	-1506.9442	3040.7045	0.0800	0.7595
SYMG	17	-1503.0907	3041.5660	0.0520	0.8115
HKYIG	17	-1503.1617	3041.7080	0.0484	0.8599
GTRI	20	-1499.9370	3041.7874	0.0465	0.9064
K2PIG	14	-1506.8858	3042.7155	0.0293	0.9357
HKYI	16	-1505.0305	3043.2889	0.0220	0.9576
SYMIG	18	-1503.0401	3043.6313	0.0185	0.9761
K2PI	13	-1508.6157	3044.0475	0.0150	0.9912
SYMI	17	-1504.8639	3045.1124	0.0088	1.0000
GTR	19	-1518.2316	3076.1905	0.0000	1.0000
HKY	15	-1523.3059	3077.6929	0.0000	1.0000
K2P	12	-1526.9912	3078.6804	0.0000	1.0000
SYM	16	-1522.9152	3079.0584	0.0000	1.0000
JC69G	12	-1536.2737	3097.2454	0.0000	1.0000
F81G	15	-1533.3196	3097.7202	0.0000	1.0000
JC69IG	13	-1536.2151	3099.2464	0.0000	1.0000
F81IG	16	-1533.2608	3099.7496	0.0000	1.0000
JC69I	12	-1538.0011	3100.7001	0.0000	1.0000
F81I	15	-1535.0865	3101.2541	0.0000	1.0000
JC69	11	-1555.9244	3134.4380	0.0000	1.0000
F81	14	-1552.9842	3134.9123	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1497.7818	3035.5636	0.4510	0.4510
GTRIG	21	-1497.7021	3037.4041	0.1797	0.6307
HKYG	16	-1503.2294	3038.4587	0.1061	0.7368
GTRI	20	-1499.9370	3039.8739	0.0523	0.7891
K2PG	13	-1506.9442	3039.8883	0.0519	0.8410
SYMG	17	-1503.0907	3040.1813	0.0448	0.8858
HKYIG	17	-1503.1617	3040.3234	0.0417	0.9275
K2PIG	14	-1506.8858	3041.7716	0.0202	0.9478
HKYI	16	-1505.0305	3042.0609	0.0175	0.9653
SYMIG	18	-1503.0401	3042.0803	0.0173	0.9826
K2PI	13	-1508.6157	3043.2313	0.0098	0.9924
SYMI	17	-1504.8639	3043.7278	0.0076	1.0000
GTR	19	-1518.2316	3074.4632	0.0000	1.0000
HKY	15	-1523.3059	3076.6118	0.0000	1.0000
SYM	16	-1522.9152	3077.8304	0.0000	1.0000
K2P	12	-1526.9912	3077.9824	0.0000	1.0000
JC69G	12	-1536.2737	3096.5474	0.0000	1.0000
F81G	15	-1533.3196	3096.6391	0.0000	1.0000
JC69IG	13	-1536.2151	3098.4303	0.0000	1.0000
F81IG	16	-1533.2608	3098.5216	0.0000	1.0000
JC69I	12	-1538.0011	3100.0021	0.0000	1.0000
F81I	15	-1535.0865	3100.1730	0.0000	1.0000
JC69	11	-1555.9244	3133.8487	0.0000	1.0000
F81	14	-1552.9842	3133.9685	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1506.9442	3093.5943	0.8047	0.8047
K2PI	13	-1508.6157	3096.9373	0.1513	0.9559
K2PIG	14	-1506.8858	3099.6088	0.0398	0.9957
HKYG	16	-1503.2294	3104.5584	0.0033	0.9991
HKYI	16	-1505.0305	3108.1606	0.0006	0.9996
SYMG	17	-1503.0907	3110.4122	0.0002	0.9998
HKYIG	17	-1503.1617	3110.5542	0.0002	1.0000
SYMI	17	-1504.8639	3113.9587	0.0000	1.0000
SYMIG	18	-1503.0401	3116.4423	0.0000	1.0000
GTRG	20	-1497.7818	3118.1881	0.0000	1.0000
GTRI	20	-1499.9370	3122.4984	0.0000	1.0000
GTRIG	21	-1497.7021	3124.1599	0.0000	1.0000
K2P	12	-1526.9912	3127.5571	0.0000	1.0000
HKY	15	-1523.3059	3138.5802	0.0000	1.0000
SYM	16	-1522.9152	3143.9300	0.0000	1.0000
JC69G	12	-1536.2737	3146.1221	0.0000	1.0000
JC69I	12	-1538.0011	3149.5769	0.0000	1.0000
JC69IG	13	-1536.2151	3152.1362	0.0000	1.0000
GTR	19	-1518.2316	3152.9565	0.0000	1.0000
F81G	15	-1533.3196	3158.6075	0.0000	1.0000
F81I	15	-1535.0865	3162.1414	0.0000	1.0000
F81IG	16	-1533.2608	3164.6212	0.0000	1.0000
JC69	11	-1555.9244	3179.2922	0.0000	1.0000
F81	14	-1552.9842	3191.8056	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/chr11_3420.phylip = 1 - 460;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3420.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/chr11_3420.phylip = 1 - 460;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3420.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/chr11_3420.phylip = 1 - 460;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3420.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

