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

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

Minimum AIC  model: K2PI
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	13	-2404.3213	4835.1372	0.4677	0.4677
HKYI	16	-2402.2420	4837.2262	0.1646	0.6323
K2PIG	14	-2404.3445	4837.2604	0.1618	0.7940
K2PG	13	-2406.0640	4838.6225	0.0819	0.8759
HKYIG	17	-2402.2634	4839.3628	0.0565	0.9324
HKYG	16	-2403.8689	4840.4799	0.0323	0.9648
SYMI	17	-2403.3773	4841.5907	0.0186	0.9833
SYMIG	18	-2403.3786	4843.6930	0.0065	0.9898
SYMG	17	-2404.5016	4843.8393	0.0060	0.9959
GTRI	20	-2402.2586	4845.6695	0.0024	0.9983
GTRG	20	-2403.2664	4847.6851	0.0009	0.9992
GTRIG	21	-2402.2584	4847.7860	0.0008	1.0000
K2P	12	-2433.5483	4891.5199	0.0000	1.0000
HKY	15	-2430.7952	4892.2443	0.0000	1.0000
SYM	16	-2433.0324	4898.8070	0.0000	1.0000
JC69I	12	-2440.8005	4906.0243	0.0000	1.0000
F81I	15	-2438.2566	4907.1672	0.0000	1.0000
JC69IG	13	-2440.8327	4908.1600	0.0000	1.0000
GTR	19	-2435.0393	4909.1197	0.0000	1.0000
F81IG	16	-2438.2849	4909.3120	0.0000	1.0000
JC69G	12	-2442.5903	4909.6040	0.0000	1.0000
F81G	15	-2439.9815	4910.6169	0.0000	1.0000
JC69	11	-2467.8475	4958.0527	0.0000	1.0000
F81	14	-2464.8891	4958.3496	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	13	-2404.3213	4834.6426	0.4451	0.4451
HKYI	16	-2402.2420	4836.4840	0.1773	0.6224
K2PIG	14	-2404.3445	4836.6890	0.1600	0.7824
K2PG	13	-2406.0640	4838.1279	0.0779	0.8603
HKYIG	17	-2402.2634	4838.5268	0.0638	0.9242
HKYG	16	-2403.8689	4839.7377	0.0348	0.9590
SYMI	17	-2403.3773	4840.7547	0.0210	0.9800
SYMIG	18	-2403.3786	4842.7573	0.0077	0.9877
SYMG	17	-2404.5016	4843.0032	0.0068	0.9945
GTRI	20	-2402.2586	4844.5172	0.0032	0.9977
GTRIG	21	-2402.2584	4846.5168	0.0012	0.9988
GTRG	20	-2403.2664	4846.5329	0.0012	1.0000
K2P	12	-2433.5483	4891.0966	0.0000	1.0000
HKY	15	-2430.7952	4891.5903	0.0000	1.0000
SYM	16	-2433.0324	4898.0648	0.0000	1.0000
JC69I	12	-2440.8005	4905.6010	0.0000	1.0000
F81I	15	-2438.2566	4906.5133	0.0000	1.0000
JC69IG	13	-2440.8327	4907.6654	0.0000	1.0000
GTR	19	-2435.0393	4908.0786	0.0000	1.0000
F81IG	16	-2438.2849	4908.5698	0.0000	1.0000
JC69G	12	-2442.5903	4909.1807	0.0000	1.0000
F81G	15	-2439.9815	4909.9630	0.0000	1.0000
JC69	11	-2467.8475	4957.6950	0.0000	1.0000
F81	14	-2464.8891	4957.7781	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-2404.3213	4894.7036	0.8256	0.8256
K2PG	13	-2406.0640	4898.1889	0.1445	0.9701
K2PIG	14	-2404.3445	4901.3700	0.0295	0.9996
HKYI	16	-2402.2420	4910.4052	0.0003	0.9999
HKYG	16	-2403.8689	4913.6589	0.0001	1.0000
HKYIG	17	-2402.2634	4917.0680	0.0000	1.0000
SYMI	17	-2403.3773	4919.2959	0.0000	1.0000
SYMG	17	-2404.5016	4921.5445	0.0000	1.0000
SYMIG	18	-2403.3786	4925.9186	0.0000	1.0000
GTRI	20	-2402.2586	4936.9187	0.0000	1.0000
GTRG	20	-2403.2664	4938.9343	0.0000	1.0000
GTRIG	21	-2402.2584	4943.5383	0.0000	1.0000
K2P	12	-2433.5483	4946.5375	0.0000	1.0000
HKY	15	-2430.7952	4960.8914	0.0000	1.0000
JC69I	12	-2440.8005	4961.0419	0.0000	1.0000
JC69G	12	-2442.5903	4964.6215	0.0000	1.0000
JC69IG	13	-2440.8327	4967.7264	0.0000	1.0000
SYM	16	-2433.0324	4971.9860	0.0000	1.0000
F81I	15	-2438.2566	4975.8144	0.0000	1.0000
F81G	15	-2439.9815	4979.2641	0.0000	1.0000
F81IG	16	-2438.2849	4982.4910	0.0000	1.0000
GTR	19	-2435.0393	4995.8600	0.0000	1.0000
JC69	11	-2467.8475	5008.5158	0.0000	1.0000
F81	14	-2464.8891	5022.4592	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_6690.phylip = 1 - 750;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6690.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_6690.phylip = 1 - 750;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6690.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_6690.phylip = 1 - 750;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6690.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

