
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:17 04/04/2011

Input data from file "683-tomodel-phylip/chr7_6588.phylip" (7 taxa, 707 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	-1592.1567	3217.1018	0.4765	0.4765
HKYG	16	-1592.4082	3217.6048	0.3705	0.8471
HKYIG	17	-1592.4082	3219.7047	0.1297	0.9767
GTRI	20	-1591.6968	3224.6182	0.0111	0.9879
GTRG	20	-1591.9450	3225.1144	0.0087	0.9965
GTRIG	21	-1591.9450	3227.2389	0.0030	0.9995
HKY	15	-1600.1252	3230.9449	0.0005	1.0000
GTR	19	-1599.8771	3238.8605	0.0000	1.0000
F81I	15	-1615.0512	3260.7971	0.0000	1.0000
F81G	15	-1615.3474	3261.3895	0.0000	1.0000
F81IG	16	-1615.3313	3263.4510	0.0000	1.0000
F81	14	-1622.1105	3272.8279	0.0000	1.0000
K2PI	13	-1639.1982	3304.9217	0.0000	1.0000
K2PG	13	-1639.4637	3305.4527	0.0000	1.0000
K2PIG	14	-1639.4637	3307.5344	0.0000	1.0000
SYMI	17	-1637.8721	3310.6325	0.0000	1.0000
SYMG	17	-1638.0895	3311.0672	0.0000	1.0000
SYMIG	18	-1638.0894	3313.1731	0.0000	1.0000
K2P	12	-1646.6136	3317.6768	0.0000	1.0000
SYM	16	-1645.4335	3323.6553	0.0000	1.0000
JC69I	12	-1658.1710	3340.7915	0.0000	1.0000
JC69G	12	-1658.4560	3341.3616	0.0000	1.0000
JC69IG	13	-1658.4560	3343.4373	0.0000	1.0000
JC69	11	-1665.1329	3352.6456	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1592.1567	3216.3134	0.4707	0.4707
HKYG	16	-1592.4082	3216.8164	0.3660	0.8367
HKYIG	17	-1592.4082	3218.8164	0.1346	0.9713
GTRI	20	-1591.6968	3223.3937	0.0137	0.9850
GTRG	20	-1591.9450	3223.8899	0.0107	0.9956
GTRIG	21	-1591.9450	3225.8900	0.0039	0.9995
HKY	15	-1600.1252	3230.2503	0.0004	1.0000
GTR	19	-1599.8771	3237.7542	0.0000	1.0000
F81I	15	-1615.0512	3260.1024	0.0000	1.0000
F81G	15	-1615.3474	3260.6948	0.0000	1.0000
F81IG	16	-1615.3313	3262.6626	0.0000	1.0000
F81	14	-1622.1105	3272.2209	0.0000	1.0000
K2PI	13	-1639.1982	3304.3964	0.0000	1.0000
K2PG	13	-1639.4637	3304.9274	0.0000	1.0000
K2PIG	14	-1639.4637	3306.9275	0.0000	1.0000
SYMI	17	-1637.8721	3309.7442	0.0000	1.0000
SYMG	17	-1638.0895	3310.1790	0.0000	1.0000
SYMIG	18	-1638.0894	3312.1789	0.0000	1.0000
K2P	12	-1646.6136	3317.2273	0.0000	1.0000
SYM	16	-1645.4335	3322.8669	0.0000	1.0000
JC69I	12	-1658.1710	3340.3420	0.0000	1.0000
JC69G	12	-1658.4560	3340.9120	0.0000	1.0000
JC69IG	13	-1658.4560	3342.9121	0.0000	1.0000
JC69	11	-1665.1329	3352.2658	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1592.1567	3289.2899	0.5506	0.5506
HKYG	16	-1592.4082	3289.7929	0.4282	0.9788
HKYIG	17	-1592.4082	3296.3539	0.0161	0.9949
HKY	15	-1600.1252	3298.6658	0.0051	1.0000
GTRI	20	-1591.6968	3314.6143	0.0000	1.0000
GTRG	20	-1591.9450	3315.1106	0.0000	1.0000
GTRIG	21	-1591.9450	3321.6716	0.0000	1.0000
GTR	19	-1599.8771	3324.4138	0.0000	1.0000
F81I	15	-1615.0512	3328.5179	0.0000	1.0000
F81G	15	-1615.3474	3329.1103	0.0000	1.0000
F81IG	16	-1615.3313	3335.6391	0.0000	1.0000
F81	14	-1622.1105	3336.0754	0.0000	1.0000
K2PI	13	-1639.1982	3363.6898	0.0000	1.0000
K2PG	13	-1639.4637	3364.2208	0.0000	1.0000
K2PIG	14	-1639.4637	3370.7819	0.0000	1.0000
K2P	12	-1646.6136	3371.9596	0.0000	1.0000
SYMI	17	-1637.8721	3387.2817	0.0000	1.0000
SYMG	17	-1638.0895	3387.7165	0.0000	1.0000
SYMIG	18	-1638.0894	3394.2774	0.0000	1.0000
JC69I	12	-1658.1710	3395.0743	0.0000	1.0000
JC69G	12	-1658.4560	3395.6444	0.0000	1.0000
SYM	16	-1645.4335	3395.8434	0.0000	1.0000
JC69IG	13	-1658.4560	3402.2055	0.0000	1.0000
JC69	11	-1665.1329	3402.4371	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 683-tomodel-phylip/chr7_6588.phylip = 1 - 707;
 Partition Dummy = 1:683-tomodel-phylip/chr7_6588.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 683-tomodel-phylip/chr7_6588.phylip = 1 - 707;
 Partition Dummy = 1:683-tomodel-phylip/chr7_6588.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 683-tomodel-phylip/chr7_6588.phylip = 1 - 707;
 Partition Dummy = 1:683-tomodel-phylip/chr7_6588.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

