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

Input data from file "591-tomodel-phylip/chr8_310.phylip" (7 taxa, 561 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	-1633.9557	3300.9114	0.4804	0.4804
HKYI	16	-1634.8364	3302.6727	0.1991	0.6796
HKYIG	17	-1633.8893	3302.9057	0.1772	0.8568
GTRG	20	-1631.4825	3304.5205	0.0791	0.9358
GTRI	20	-1632.2847	3306.1249	0.0354	0.9713
GTRIG	21	-1631.4162	3306.5466	0.0287	1.0000
HKY	15	-1646.6054	3324.0915	0.0000	1.0000
GTR	19	-1643.6586	3326.7220	0.0000	1.0000
SYMG	17	-1653.0947	3341.3166	0.0000	1.0000
SYMI	17	-1653.6567	3342.4404	0.0000	1.0000
SYMIG	18	-1652.9826	3343.2271	0.0000	1.0000
K2PG	13	-1658.3889	3343.4431	0.0000	1.0000
K2PI	13	-1658.9427	3344.5509	0.0000	1.0000
K2PIG	14	-1658.3048	3345.3789	0.0000	1.0000
SYM	16	-1665.8674	3364.7347	0.0000	1.0000
K2P	12	-1671.0790	3366.7274	0.0000	1.0000
F81G	15	-1670.9323	3372.7453	0.0000	1.0000
F81I	15	-1671.6703	3374.2214	0.0000	1.0000
F81IG	16	-1670.8962	3374.7923	0.0000	1.0000
F81	14	-1682.1189	3393.0071	0.0000	1.0000
JC69G	12	-1691.6683	3407.9060	0.0000	1.0000
JC69I	12	-1692.3532	3409.2758	0.0000	1.0000
JC69IG	13	-1691.6236	3409.9127	0.0000	1.0000
JC69	11	-1703.0839	3428.6488	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1633.9557	3299.9114	0.4530	0.4530
HKYI	16	-1634.8364	3301.6727	0.1878	0.6407
HKYIG	17	-1633.8893	3301.7786	0.1781	0.8188
GTRG	20	-1631.4825	3302.9649	0.0984	0.9172
GTRI	20	-1632.2847	3304.5693	0.0441	0.9613
GTRIG	21	-1631.4162	3304.8323	0.0387	1.0000
HKY	15	-1646.6054	3323.2108	0.0000	1.0000
GTR	19	-1643.6586	3325.3172	0.0000	1.0000
SYMG	17	-1653.0947	3340.1895	0.0000	1.0000
SYMI	17	-1653.6567	3341.3134	0.0000	1.0000
SYMIG	18	-1652.9826	3341.9651	0.0000	1.0000
K2PG	13	-1658.3889	3342.7777	0.0000	1.0000
K2PI	13	-1658.9427	3343.8855	0.0000	1.0000
K2PIG	14	-1658.3048	3344.6096	0.0000	1.0000
SYM	16	-1665.8674	3363.7347	0.0000	1.0000
K2P	12	-1671.0790	3366.1581	0.0000	1.0000
F81G	15	-1670.9323	3371.8646	0.0000	1.0000
F81I	15	-1671.6703	3373.3407	0.0000	1.0000
F81IG	16	-1670.8962	3373.7923	0.0000	1.0000
F81	14	-1682.1189	3392.2378	0.0000	1.0000
JC69G	12	-1691.6683	3407.3367	0.0000	1.0000
JC69I	12	-1692.3532	3408.7064	0.0000	1.0000
JC69IG	13	-1691.6236	3409.2473	0.0000	1.0000
JC69	11	-1703.0839	3428.1679	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1633.9557	3369.1870	0.6850	0.6850
HKYI	16	-1634.8364	3370.9482	0.2840	0.9690
HKYIG	17	-1633.8893	3375.3839	0.0309	0.9999
HKY	15	-1646.6054	3388.1566	0.0001	1.0000
GTRG	20	-1631.4825	3389.5594	0.0000	1.0000
GTRI	20	-1632.2847	3391.1638	0.0000	1.0000
GTRIG	21	-1631.4162	3395.7565	0.0000	1.0000
K2PG	13	-1658.3889	3399.0641	0.0000	1.0000
K2PI	13	-1658.9427	3400.1718	0.0000	1.0000
K2PIG	14	-1658.3048	3405.2257	0.0000	1.0000
GTR	19	-1643.6586	3407.5819	0.0000	1.0000
SYMG	17	-1653.0947	3413.7947	0.0000	1.0000
SYMI	17	-1653.6567	3414.9186	0.0000	1.0000
K2P	12	-1671.0790	3418.1147	0.0000	1.0000
SYMIG	18	-1652.9826	3419.9001	0.0000	1.0000
SYM	16	-1665.8674	3433.0103	0.0000	1.0000
F81G	15	-1670.9323	3436.8104	0.0000	1.0000
F81I	15	-1671.6703	3438.2865	0.0000	1.0000
F81IG	16	-1670.8962	3443.0679	0.0000	1.0000
F81	14	-1682.1189	3452.8539	0.0000	1.0000
JC69G	12	-1691.6683	3459.2933	0.0000	1.0000
JC69I	12	-1692.3532	3460.6631	0.0000	1.0000
JC69IG	13	-1691.6236	3465.5337	0.0000	1.0000
JC69	11	-1703.0839	3475.7948	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 591-tomodel-phylip/chr8_310.phylip = 1 - 561;
 Partition Dummy = 1:591-tomodel-phylip/chr8_310.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 591-tomodel-phylip/chr8_310.phylip = 1 - 561;
 Partition Dummy = 1:591-tomodel-phylip/chr8_310.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 591-tomodel-phylip/chr8_310.phylip = 1 - 561;
 Partition Dummy = 1:591-tomodel-phylip/chr8_310.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

