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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-2821.8484	5684.9124	0.2159	0.2159
SYMG	17	-2825.3052	5685.4922	0.1616	0.3774
HKYG	16	-2826.6899	5686.1625	0.1155	0.4930
GTRIG	21	-2821.6373	5686.6138	0.0922	0.5852
GTRI	20	-2822.9706	5687.1567	0.0703	0.6555
SYMIG	18	-2825.0867	5687.1604	0.0702	0.7256
K2PG	13	-2830.3678	5687.2572	0.0668	0.7925
SYMI	17	-2826.4123	5687.7065	0.0534	0.8459
HKYI	16	-2827.5338	5687.8503	0.0497	0.8956
HKYIG	17	-2826.4904	5687.8626	0.0494	0.9450
K2PIG	14	-2830.2057	5689.0139	0.0278	0.9727
K2PI	13	-2831.2641	5689.0497	0.0273	1.0000
F81I	15	-2879.4582	5789.6061	0.0000	1.0000
F81G	15	-2879.7348	5790.1592	0.0000	1.0000
F81IG	16	-2879.1297	5791.0421	0.0000	1.0000
JC69I	12	-2885.4373	5795.3209	0.0000	1.0000
JC69G	12	-2885.7017	5795.8498	0.0000	1.0000
JC69IG	13	-2885.1014	5796.7243	0.0000	1.0000
GTR	19	-2879.4977	5798.0936	0.0000	1.0000
SYM	16	-2882.9668	5798.7164	0.0000	1.0000
HKY	15	-2884.4035	5799.4967	0.0000	1.0000
K2P	12	-2888.3251	5801.0966	0.0000	1.0000
F81	14	-2933.1290	5894.8605	0.0000	1.0000
JC69	11	-2938.3512	5899.0796	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2821.8484	5683.6968	0.2424	0.2424
SYMG	17	-2825.3052	5684.6104	0.1535	0.3959
GTRIG	21	-2821.6373	5685.2746	0.1101	0.5061
HKYG	16	-2826.6899	5685.3798	0.1045	0.6106
GTRI	20	-2822.9706	5685.9411	0.0789	0.6895
SYMIG	18	-2825.0867	5686.1733	0.0703	0.7597
K2PG	13	-2830.3678	5686.7357	0.0530	0.8128
SYMI	17	-2826.4123	5686.8246	0.0507	0.8635
HKYIG	17	-2826.4904	5686.9808	0.0469	0.9105
HKYI	16	-2827.5338	5687.0676	0.0449	0.9554
K2PIG	14	-2830.2057	5688.4114	0.0230	0.9784
K2PI	13	-2831.2641	5688.5282	0.0216	1.0000
F81I	15	-2879.4582	5788.9165	0.0000	1.0000
F81G	15	-2879.7348	5789.4695	0.0000	1.0000
F81IG	16	-2879.1297	5790.2593	0.0000	1.0000
JC69I	12	-2885.4373	5794.8746	0.0000	1.0000
JC69G	12	-2885.7017	5795.4034	0.0000	1.0000
JC69IG	13	-2885.1014	5796.2029	0.0000	1.0000
GTR	19	-2879.4977	5796.9953	0.0000	1.0000
SYM	16	-2882.9668	5797.9337	0.0000	1.0000
HKY	15	-2884.4035	5798.8070	0.0000	1.0000
K2P	12	-2888.3251	5800.6502	0.0000	1.0000
F81	14	-2933.1290	5894.2580	0.0000	1.0000
JC69	11	-2938.3512	5898.7024	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2830.3678	5746.1207	0.6870	0.6870
K2PI	13	-2831.2641	5747.9132	0.2803	0.9673
K2PIG	14	-2830.2057	5752.3645	0.0303	0.9976
HKYG	16	-2826.6899	5758.4690	0.0014	0.9990
HKYI	16	-2827.5338	5760.1568	0.0006	0.9996
SYMG	17	-2825.3052	5762.2677	0.0002	0.9999
SYMI	17	-2826.4123	5764.4819	0.0001	0.9999
HKYIG	17	-2826.4904	5764.6381	0.0001	1.0000
SYMIG	18	-2825.0867	5768.3987	0.0000	1.0000
GTRG	20	-2821.8484	5775.0584	0.0000	1.0000
GTRI	20	-2822.9706	5777.3027	0.0000	1.0000
GTRIG	21	-2821.6373	5781.2043	0.0000	1.0000
JC69I	12	-2885.4373	5849.6915	0.0000	1.0000
JC69G	12	-2885.7017	5850.2204	0.0000	1.0000
K2P	12	-2888.3251	5855.4672	0.0000	1.0000
JC69IG	13	-2885.1014	5855.5879	0.0000	1.0000
F81I	15	-2879.4582	5857.4376	0.0000	1.0000
F81G	15	-2879.7348	5857.9907	0.0000	1.0000
F81IG	16	-2879.1297	5863.3486	0.0000	1.0000
HKY	15	-2884.4035	5867.3282	0.0000	1.0000
SYM	16	-2882.9668	5871.0229	0.0000	1.0000
GTR	19	-2879.4977	5883.7888	0.0000	1.0000
JC69	11	-2938.3512	5948.9513	0.0000	1.0000
F81	14	-2933.1290	5958.2111	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/chr7_4728.phylip = 1 - 712;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4728.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/chr7_4728.phylip = 1 - 712;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4728.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 (K2PG)]

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

