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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1578.3270	3198.4412	0.3769	0.3769
GTRG	20	-1578.4782	3198.7436	0.3240	0.7009
GTRIG	21	-1578.2831	3200.5363	0.1322	0.8331
HKYI	16	-1584.1968	3201.5413	0.0800	0.9131
HKYG	16	-1584.5280	3202.2037	0.0574	0.9705
HKYIG	17	-1584.1231	3203.5401	0.0294	1.0000
K2PI	13	-1597.4547	3221.6724	0.0000	1.0000
K2PG	13	-1597.9623	3222.6877	0.0000	1.0000
K2PIG	14	-1597.4381	3223.7586	0.0000	1.0000
SYMI	17	-1595.9198	3227.1334	0.0000	1.0000
SYMG	17	-1596.3167	3227.9273	0.0000	1.0000
SYMIG	18	-1595.9723	3229.3938	0.0000	1.0000
GTR	19	-1599.8959	3239.4054	0.0000	1.0000
HKY	15	-1607.2172	3245.4449	0.0000	1.0000
F81I	15	-1612.3356	3255.6818	0.0000	1.0000
F81G	15	-1612.5366	3256.0836	0.0000	1.0000
F81IG	16	-1612.2619	3257.6714	0.0000	1.0000
K2P	12	-1622.9528	3270.5583	0.0000	1.0000
SYM	16	-1619.7032	3272.5540	0.0000	1.0000
JC69I	12	-1624.3937	3273.4402	0.0000	1.0000
JC69G	12	-1624.5837	3273.8200	0.0000	1.0000
JC69IG	13	-1624.3236	3275.4102	0.0000	1.0000
F81	14	-1635.4989	3299.8803	0.0000	1.0000
JC69	11	-1648.7166	3319.9844	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1578.3270	3196.6539	0.3891	0.3891
GTRG	20	-1578.4782	3196.9563	0.3345	0.7236
GTRIG	21	-1578.2831	3198.5662	0.1496	0.8732
HKYI	16	-1584.1968	3200.3936	0.0600	0.9332
HKYG	16	-1584.5280	3201.0560	0.0431	0.9762
HKYIG	17	-1584.1231	3202.2462	0.0238	1.0000
K2PI	13	-1597.4547	3220.9093	0.0000	1.0000
K2PG	13	-1597.9623	3221.9246	0.0000	1.0000
K2PIG	14	-1597.4381	3222.8762	0.0000	1.0000
SYMI	17	-1595.9198	3225.8396	0.0000	1.0000
SYMG	17	-1596.3167	3226.6335	0.0000	1.0000
SYMIG	18	-1595.9723	3227.9447	0.0000	1.0000
GTR	19	-1599.8959	3237.7918	0.0000	1.0000
HKY	15	-1607.2172	3244.4344	0.0000	1.0000
F81I	15	-1612.3356	3254.6713	0.0000	1.0000
F81G	15	-1612.5366	3255.0731	0.0000	1.0000
F81IG	16	-1612.2619	3256.5237	0.0000	1.0000
K2P	12	-1622.9528	3269.9056	0.0000	1.0000
SYM	16	-1619.7032	3271.4063	0.0000	1.0000
JC69I	12	-1624.3937	3272.7874	0.0000	1.0000
JC69G	12	-1624.5837	3273.1673	0.0000	1.0000
JC69IG	13	-1624.3236	3274.6471	0.0000	1.0000
F81	14	-1635.4989	3298.9979	0.0000	1.0000
JC69	11	-1648.7166	3319.4333	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1584.1968	3267.5367	0.5553	0.5553
HKYG	16	-1584.5280	3268.1991	0.3988	0.9541
HKYIG	17	-1584.1231	3273.5858	0.0270	0.9811
K2PI	13	-1597.4547	3275.4631	0.0106	0.9916
K2PG	13	-1597.9623	3276.4784	0.0064	0.9980
GTRI	20	-1578.3270	3280.5828	0.0008	0.9988
GTRG	20	-1578.4782	3280.8852	0.0007	0.9995
K2PIG	14	-1597.4381	3281.6264	0.0005	1.0000
GTRIG	21	-1578.2831	3286.6915	0.0000	1.0000
SYMI	17	-1595.9198	3297.1791	0.0000	1.0000
SYMG	17	-1596.3167	3297.9730	0.0000	1.0000
SYMIG	18	-1595.9723	3303.4807	0.0000	1.0000
HKY	15	-1607.2172	3307.3811	0.0000	1.0000
GTR	19	-1599.8959	3317.5242	0.0000	1.0000
F81I	15	-1612.3356	3317.6179	0.0000	1.0000
F81G	15	-1612.5366	3318.0198	0.0000	1.0000
K2P	12	-1622.9528	3320.2629	0.0000	1.0000
JC69I	12	-1624.3937	3323.1448	0.0000	1.0000
JC69G	12	-1624.5837	3323.5246	0.0000	1.0000
F81IG	16	-1612.2619	3323.6668	0.0000	1.0000
JC69IG	13	-1624.3236	3329.2009	0.0000	1.0000
SYM	16	-1619.7032	3338.5494	0.0000	1.0000
F81	14	-1635.4989	3357.7481	0.0000	1.0000
JC69	11	-1648.7166	3365.5942	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11660.phylip = 1 - 491;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11660.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11660.phylip = 1 - 491;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11660.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11660.phylip = 1 - 491;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11660.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

