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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1626.8759	3295.2001	0.4808	0.4808
GTRIG	21	-1626.5209	3296.6376	0.2343	0.7151
GTRI	20	-1628.3948	3298.2378	0.1053	0.8204
HKYG	16	-1632.6666	3298.2647	0.1039	0.9243
HKYIG	17	-1632.2494	3299.5486	0.0547	0.9789
HKYI	16	-1634.2615	3301.4546	0.0211	1.0000
GTR	19	-1649.4678	3338.2437	0.0000	1.0000
K2PG	13	-1657.2552	3341.1306	0.0000	1.0000
HKY	15	-1655.2658	3341.3520	0.0000	1.0000
K2PIG	14	-1656.8256	3342.3680	0.0000	1.0000
SYMG	17	-1654.1073	3343.2642	0.0000	1.0000
K2PI	13	-1658.8438	3344.3077	0.0000	1.0000
SYMIG	18	-1653.7470	3344.6693	0.0000	1.0000
SYMI	17	-1655.5550	3346.1598	0.0000	1.0000
F81G	15	-1670.4001	3371.6208	0.0000	1.0000
F81IG	16	-1670.0318	3372.9951	0.0000	1.0000
F81I	15	-1671.8944	3374.6093	0.0000	1.0000
K2P	12	-1679.3695	3383.2697	0.0000	1.0000
SYM	16	-1675.7108	3384.3532	0.0000	1.0000
JC69G	12	-1690.1939	3404.9183	0.0000	1.0000
JC69IG	13	-1689.8014	3406.2230	0.0000	1.0000
JC69I	12	-1691.7382	3408.0069	0.0000	1.0000
F81	14	-1691.5705	3411.8577	0.0000	1.0000
JC69	11	-1711.4165	3445.2812	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1626.8759	3293.7518	0.4908	0.4908
GTRIG	21	-1626.5209	3295.0417	0.2575	0.7483
GTRI	20	-1628.3948	3296.7895	0.1075	0.8558
HKYG	16	-1632.6666	3297.3332	0.0819	0.9377
HKYIG	17	-1632.2494	3298.4989	0.0457	0.9834
HKYI	16	-1634.2615	3300.5231	0.0166	1.0000
GTR	19	-1649.4678	3336.9356	0.0000	1.0000
K2PG	13	-1657.2552	3340.5105	0.0000	1.0000
HKY	15	-1655.2658	3340.5315	0.0000	1.0000
K2PIG	14	-1656.8256	3341.6513	0.0000	1.0000
SYMG	17	-1654.1073	3342.2145	0.0000	1.0000
SYMIG	18	-1653.7470	3343.4940	0.0000	1.0000
K2PI	13	-1658.8438	3343.6876	0.0000	1.0000
SYMI	17	-1655.5550	3345.1100	0.0000	1.0000
F81G	15	-1670.4001	3370.8003	0.0000	1.0000
F81IG	16	-1670.0318	3372.0636	0.0000	1.0000
F81I	15	-1671.8944	3373.7888	0.0000	1.0000
K2P	12	-1679.3695	3382.7390	0.0000	1.0000
SYM	16	-1675.7108	3383.4217	0.0000	1.0000
JC69G	12	-1690.1939	3404.3877	0.0000	1.0000
JC69IG	13	-1689.8014	3405.6029	0.0000	1.0000
JC69I	12	-1691.7382	3407.4763	0.0000	1.0000
F81	14	-1691.5705	3411.1409	0.0000	1.0000
JC69	11	-1711.4165	3444.8330	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1632.6666	3367.7107	0.7899	0.7899
HKYI	16	-1634.2615	3370.9006	0.1603	0.9502
HKYIG	17	-1632.2494	3373.2750	0.0489	0.9991
GTRG	20	-1626.8759	3381.7237	0.0007	0.9998
GTRI	20	-1628.3948	3384.7614	0.0002	1.0000
GTRIG	21	-1626.5209	3387.4122	0.0000	1.0000
K2PG	13	-1657.2552	3397.6922	0.0000	1.0000
K2PI	13	-1658.8438	3400.8693	0.0000	1.0000
K2PIG	14	-1656.8256	3403.2316	0.0000	1.0000
HKY	15	-1655.2658	3406.5104	0.0000	1.0000
SYMG	17	-1654.1073	3416.9906	0.0000	1.0000
SYMI	17	-1655.5550	3419.8862	0.0000	1.0000
GTR	19	-1649.4678	3420.5089	0.0000	1.0000
SYMIG	18	-1653.7470	3422.6687	0.0000	1.0000
K2P	12	-1679.3695	3435.5222	0.0000	1.0000
F81G	15	-1670.4001	3436.7792	0.0000	1.0000
F81I	15	-1671.8944	3439.7677	0.0000	1.0000
F81IG	16	-1670.0318	3442.4411	0.0000	1.0000
SYM	16	-1675.7108	3453.7992	0.0000	1.0000
JC69G	12	-1690.1939	3457.1709	0.0000	1.0000
JC69I	12	-1691.7382	3460.2595	0.0000	1.0000
JC69IG	13	-1689.8014	3462.7846	0.0000	1.0000
F81	14	-1691.5705	3472.7213	0.0000	1.0000
JC69	11	-1711.4165	3493.2175	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_4403.phylip = 1 - 601;
 Partition Dummy = 1:591-tomodel-phylip/chr6_4403.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_4403.phylip = 1 - 601;
 Partition Dummy = 1:591-tomodel-phylip/chr6_4403.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

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

