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

Input data from file "683-tomodel-phylip/chr9_3515.phylip" (7 taxa, 532 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	16	-1708.9871	3451.0304	0.4141	0.4141
HKYI	16	-1709.3054	3451.6671	0.3012	0.7153
HKYIG	17	-1709.0115	3453.2136	0.1390	0.8543
GTRG	20	-1706.5805	3454.8048	0.0627	0.9170
GTRI	20	-1706.7962	3455.2362	0.0506	0.9675
GTRIG	21	-1706.5896	3456.9909	0.0210	0.9886
HKY	15	-1713.8126	3458.5554	0.0096	0.9982
GTR	19	-1711.2050	3461.8944	0.0018	1.0000
K2PG	13	-1728.0306	3482.7638	0.0000	1.0000
K2PI	13	-1728.5732	3483.8491	0.0000	1.0000
K2PIG	14	-1728.0306	3484.8735	0.0000	1.0000
SYMG	17	-1727.1537	3489.4981	0.0000	1.0000
SYMI	17	-1727.6768	3490.5443	0.0000	1.0000
SYMIG	18	-1727.1680	3491.6694	0.0000	1.0000
K2P	12	-1733.6978	3491.9967	0.0000	1.0000
SYM	16	-1732.4726	3498.0015	0.0000	1.0000
F81G	15	-1765.5077	3561.9456	0.0000	1.0000
F81I	15	-1765.6035	3562.1371	0.0000	1.0000
F81IG	16	-1765.5160	3564.0883	0.0000	1.0000
F81	14	-1769.8063	3568.4250	0.0000	1.0000
JC69G	12	-1782.5776	3589.7563	0.0000	1.0000
JC69I	12	-1782.8555	3590.3122	0.0000	1.0000
JC69IG	13	-1782.6013	3591.9052	0.0000	1.0000
JC69	11	-1787.3211	3597.1498	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1708.9871	3449.9741	0.3914	0.3914
HKYI	16	-1709.3054	3450.6108	0.2847	0.6762
HKYIG	17	-1709.0115	3452.0230	0.1405	0.8167
GTRG	20	-1706.5805	3453.1610	0.0796	0.8962
GTRI	20	-1706.7962	3453.5924	0.0641	0.9603
GTRIG	21	-1706.5896	3455.1792	0.0290	0.9893
HKY	15	-1713.8126	3457.6252	0.0085	0.9979
GTR	19	-1711.2050	3460.4100	0.0021	1.0000
K2PG	13	-1728.0306	3482.0611	0.0000	1.0000
K2PI	13	-1728.5732	3483.1464	0.0000	1.0000
K2PIG	14	-1728.0306	3484.0611	0.0000	1.0000
SYMG	17	-1727.1537	3488.3074	0.0000	1.0000
SYMI	17	-1727.6768	3489.3536	0.0000	1.0000
SYMIG	18	-1727.1680	3490.3361	0.0000	1.0000
K2P	12	-1733.6978	3491.3955	0.0000	1.0000
SYM	16	-1732.4726	3496.9451	0.0000	1.0000
F81G	15	-1765.5077	3561.0154	0.0000	1.0000
F81I	15	-1765.6035	3561.2069	0.0000	1.0000
F81IG	16	-1765.5160	3563.0320	0.0000	1.0000
F81	14	-1769.8063	3567.6126	0.0000	1.0000
JC69G	12	-1782.5776	3589.1552	0.0000	1.0000
JC69I	12	-1782.8555	3589.7111	0.0000	1.0000
JC69IG	13	-1782.6013	3591.2025	0.0000	1.0000
JC69	11	-1787.3211	3596.6421	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1708.9871	3518.4004	0.5115	0.5115
HKYI	16	-1709.3054	3519.0371	0.3721	0.8836
HKY	15	-1713.8126	3521.7748	0.0947	0.9783
HKYIG	17	-1709.0115	3524.7259	0.0216	0.9999
K2PG	13	-1728.0306	3537.6575	0.0000	0.9999
GTRG	20	-1706.5805	3538.6938	0.0000	1.0000
K2PI	13	-1728.5732	3538.7428	0.0000	1.0000
GTRI	20	-1706.7962	3539.1253	0.0000	1.0000
GTR	19	-1711.2050	3541.6662	0.0000	1.0000
K2P	12	-1733.6978	3542.7152	0.0000	1.0000
K2PIG	14	-1728.0306	3543.9341	0.0000	1.0000
GTRIG	21	-1706.5896	3544.9887	0.0000	1.0000
SYMG	17	-1727.1537	3561.0103	0.0000	1.0000
SYMI	17	-1727.6768	3562.0565	0.0000	1.0000
SYM	16	-1732.4726	3565.3714	0.0000	1.0000
SYMIG	18	-1727.1680	3567.3157	0.0000	1.0000
F81G	15	-1765.5077	3625.1650	0.0000	1.0000
F81I	15	-1765.6035	3625.3566	0.0000	1.0000
F81	14	-1769.8063	3627.4856	0.0000	1.0000
F81IG	16	-1765.5160	3631.4583	0.0000	1.0000
JC69G	12	-1782.5776	3640.4749	0.0000	1.0000
JC69I	12	-1782.8555	3641.0308	0.0000	1.0000
JC69	11	-1787.3211	3643.6852	0.0000	1.0000
JC69IG	13	-1782.6013	3646.7989	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr9_3515.phylip = 1 - 532;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3515.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr9_3515.phylip = 1 - 532;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3515.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr9_3515.phylip = 1 - 532;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3515.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

