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

Input data from file "591-tomodel-phylip/chr2_24807.phylip" (7 taxa, 476 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	-1580.7161	3194.6173	0.4267	0.4267
HKYG	16	-1581.1471	3195.4793	0.2773	0.7040
HKYIG	17	-1580.7155	3196.7673	0.1456	0.8497
GTRI	20	-1578.1479	3198.1419	0.0732	0.9229
GTRG	20	-1578.4839	3198.8140	0.0523	0.9753
GTRIG	21	-1578.1396	3200.3145	0.0247	1.0000
HKY	15	-1592.9648	3216.9730	0.0000	1.0000
GTR	19	-1590.4954	3220.6574	0.0000	1.0000
K2PI	13	-1598.5059	3223.7997	0.0000	1.0000
K2PG	13	-1599.4959	3225.7796	0.0000	1.0000
K2PIG	14	-1598.5205	3225.9520	0.0000	1.0000
SYMI	17	-1596.0127	3227.3616	0.0000	1.0000
SYMG	17	-1597.0331	3229.4025	0.0000	1.0000
SYMIG	18	-1596.0281	3229.5529	0.0000	1.0000
F81I	15	-1601.0468	3233.1371	0.0000	1.0000
F81G	15	-1601.7591	3234.5618	0.0000	1.0000
F81IG	16	-1601.0538	3235.2928	0.0000	1.0000
K2P	12	-1610.9532	3246.5803	0.0000	1.0000
SYM	16	-1608.4784	3250.1419	0.0000	1.0000
F81	14	-1612.8128	3254.5366	0.0000	1.0000
JC69I	12	-1616.5322	3257.7382	0.0000	1.0000
JC69G	12	-1617.4889	3259.6516	0.0000	1.0000
JC69IG	13	-1616.5460	3259.8800	0.0000	1.0000
JC69	11	-1628.2952	3279.1595	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1580.7161	3193.4321	0.3974	0.3974
HKYG	16	-1581.1471	3194.2941	0.2583	0.6557
HKYIG	17	-1580.7155	3195.4311	0.1463	0.8020
GTRI	20	-1578.1479	3196.2958	0.0949	0.8969
GTRG	20	-1578.4839	3196.9678	0.0678	0.9648
GTRIG	21	-1578.1396	3198.2792	0.0352	1.0000
HKY	15	-1592.9648	3215.9296	0.0000	1.0000
GTR	19	-1590.4954	3218.9908	0.0000	1.0000
K2PI	13	-1598.5059	3223.0118	0.0000	1.0000
K2PG	13	-1599.4959	3224.9917	0.0000	1.0000
K2PIG	14	-1598.5205	3225.0409	0.0000	1.0000
SYMI	17	-1596.0127	3226.0253	0.0000	1.0000
SYMIG	18	-1596.0281	3228.0562	0.0000	1.0000
SYMG	17	-1597.0331	3228.0663	0.0000	1.0000
F81I	15	-1601.0468	3232.0936	0.0000	1.0000
F81G	15	-1601.7591	3233.5183	0.0000	1.0000
F81IG	16	-1601.0538	3234.1077	0.0000	1.0000
K2P	12	-1610.9532	3245.9065	0.0000	1.0000
SYM	16	-1608.4784	3248.9568	0.0000	1.0000
F81	14	-1612.8128	3253.6255	0.0000	1.0000
JC69I	12	-1616.5322	3257.0643	0.0000	1.0000
JC69G	12	-1617.4889	3258.9777	0.0000	1.0000
JC69IG	13	-1616.5460	3259.0921	0.0000	1.0000
JC69	11	-1628.2952	3278.5905	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1580.7161	3260.0788	0.5896	0.5896
HKYG	16	-1581.1471	3260.9408	0.3831	0.9727
HKYIG	17	-1580.7155	3266.2432	0.0270	0.9997
K2PI	13	-1598.5059	3277.1622	0.0001	0.9998
HKY	15	-1592.9648	3278.4108	0.0001	0.9999
K2PG	13	-1599.4959	3279.1421	0.0000	0.9999
GTRI	20	-1578.1479	3279.6041	0.0000	1.0000
GTRG	20	-1578.4839	3280.2762	0.0000	1.0000
K2PIG	14	-1598.5205	3283.3568	0.0000	1.0000
GTRIG	21	-1578.1396	3285.7530	0.0000	1.0000
F81I	15	-1601.0468	3294.5749	0.0000	1.0000
K2P	12	-1610.9532	3295.8915	0.0000	1.0000
F81G	15	-1601.7591	3295.9995	0.0000	1.0000
SYMI	17	-1596.0127	3296.8374	0.0000	1.0000
GTR	19	-1590.4954	3298.1337	0.0000	1.0000
SYMG	17	-1597.0331	3298.8784	0.0000	1.0000
F81IG	16	-1601.0538	3300.7543	0.0000	1.0000
SYMIG	18	-1596.0281	3303.0337	0.0000	1.0000
JC69I	12	-1616.5322	3307.0493	0.0000	1.0000
JC69G	12	-1617.4889	3308.9628	0.0000	1.0000
F81	14	-1612.8128	3311.9414	0.0000	1.0000
JC69IG	13	-1616.5460	3313.2425	0.0000	1.0000
SYM	16	-1608.4784	3315.6034	0.0000	1.0000
JC69	11	-1628.2952	3324.4101	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/chr2_24807.phylip = 1 - 476;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24807.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/chr2_24807.phylip = 1 - 476;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24807.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/chr2_24807.phylip = 1 - 476;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24807.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

