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

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

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	13	-1683.6166	3394.0700	0.3667	0.3667
HKYG	16	-1680.8679	3394.9950	0.2309	0.5976
K2PIG	14	-1683.6002	3396.1682	0.1284	0.7260
HKYIG	17	-1680.8516	3397.1231	0.0797	0.8057
K2PI	13	-1685.4632	3397.7632	0.0578	0.8635
SYMG	17	-1681.6383	3398.6965	0.0363	0.8998
GTRG	20	-1678.4229	3398.8083	0.0343	0.9341
HKYI	16	-1682.8717	3399.0026	0.0311	0.9652
SYMIG	18	-1681.6147	3400.8201	0.0125	0.9778
GTRIG	21	-1678.3975	3400.9589	0.0117	0.9895
SYMI	17	-1683.5257	3402.4714	0.0055	0.9950
GTRI	20	-1680.3418	3402.6462	0.0050	1.0000
K2P	12	-1700.4905	3425.6966	0.0000	1.0000
HKY	15	-1698.3079	3427.7244	0.0000	1.0000
SYM	16	-1698.3268	3429.9130	0.0000	1.0000
GTR	19	-1695.4783	3430.7282	0.0000	1.0000
JC69G	12	-1727.7411	3480.1978	0.0000	1.0000
F81G	15	-1725.0987	3481.3059	0.0000	1.0000
JC69IG	13	-1727.7400	3482.3167	0.0000	1.0000
F81IG	16	-1725.0992	3483.4577	0.0000	1.0000
JC69I	12	-1729.7713	3484.2581	0.0000	1.0000
F81I	15	-1727.0772	3485.2628	0.0000	1.0000
JC69	11	-1744.3679	3511.3398	0.0000	1.0000
F81	14	-1741.6650	3512.2977	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	13	-1683.6166	3393.2332	0.3167	0.3167
HKYG	16	-1680.8679	3393.7358	0.2463	0.5630
K2PIG	14	-1683.6002	3395.2005	0.1184	0.6814
HKYIG	17	-1680.8516	3395.7031	0.0921	0.7735
GTRG	20	-1678.4229	3396.8457	0.0520	0.8255
K2PI	13	-1685.4632	3396.9264	0.0500	0.8754
SYMG	17	-1681.6383	3397.2765	0.0419	0.9174
HKYI	16	-1682.8717	3397.7433	0.0332	0.9506
GTRIG	21	-1678.3975	3398.7949	0.0196	0.9702
SYMIG	18	-1681.6147	3399.2294	0.0158	0.9860
GTRI	20	-1680.3418	3400.6836	0.0076	0.9936
SYMI	17	-1683.5257	3401.0515	0.0064	1.0000
K2P	12	-1700.4905	3424.9810	0.0000	1.0000
HKY	15	-1698.3079	3426.6158	0.0000	1.0000
SYM	16	-1698.3268	3428.6537	0.0000	1.0000
GTR	19	-1695.4783	3428.9567	0.0000	1.0000
JC69G	12	-1727.7411	3479.4822	0.0000	1.0000
F81G	15	-1725.0987	3480.1973	0.0000	1.0000
JC69IG	13	-1727.7400	3481.4799	0.0000	1.0000
F81IG	16	-1725.0992	3482.1984	0.0000	1.0000
JC69I	12	-1729.7713	3483.5425	0.0000	1.0000
F81I	15	-1727.0772	3484.1543	0.0000	1.0000
JC69	11	-1744.3679	3510.7357	0.0000	1.0000
F81	14	-1741.6650	3511.3299	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1683.6166	3446.6245	0.8280	0.8280
K2PI	13	-1685.4632	3450.3177	0.1306	0.9586
K2PIG	14	-1683.6002	3452.6988	0.0397	0.9984
HKYG	16	-1680.8679	3459.4481	0.0014	0.9997
HKYI	16	-1682.8717	3463.4557	0.0002	0.9999
HKYIG	17	-1680.8516	3465.5225	0.0001	1.0000
SYMG	17	-1681.6383	3467.0959	0.0000	1.0000
SYMI	17	-1683.5257	3470.8708	0.0000	1.0000
SYMIG	18	-1681.6147	3473.1558	0.0000	1.0000
K2P	12	-1700.4905	3474.2653	0.0000	1.0000
GTRG	20	-1678.4229	3478.9862	0.0000	1.0000
GTRI	20	-1680.3418	3482.8240	0.0000	1.0000
GTRIG	21	-1678.3975	3485.0424	0.0000	1.0000
HKY	15	-1698.3079	3488.2212	0.0000	1.0000
SYM	16	-1698.3268	3494.3661	0.0000	1.0000
GTR	19	-1695.4783	3506.9901	0.0000	1.0000
JC69G	12	-1727.7411	3528.7665	0.0000	1.0000
JC69I	12	-1729.7713	3532.8268	0.0000	1.0000
JC69IG	13	-1727.7400	3534.8712	0.0000	1.0000
F81G	15	-1725.0987	3541.8027	0.0000	1.0000
F81I	15	-1727.0772	3545.7596	0.0000	1.0000
F81IG	16	-1725.0992	3547.9108	0.0000	1.0000
JC69	11	-1744.3679	3555.9130	0.0000	1.0000
F81	14	-1741.6650	3568.8282	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_3262.phylip = 1 - 449;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3262.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_3262.phylip = 1 - 449;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3262.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_3262.phylip = 1 - 449;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3262.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

