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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	20	-2313.2264	4667.7104	0.5457	0.5457
GTRG	20	-2314.0546	4669.3668	0.2384	0.7841
GTRIG	21	-2313.0897	4669.5646	0.2159	1.0000
HKYI	16	-2335.2681	4703.3457	0.0000	1.0000
HKYIG	17	-2335.2668	4705.4456	0.0000	1.0000
HKYG	16	-2336.3592	4705.5279	0.0000	1.0000
GTR	19	-2341.3260	4721.7880	0.0000	1.0000
SYMI	17	-2345.8952	4726.7024	0.0000	1.0000
SYMG	17	-2346.9209	4728.7540	0.0000	1.0000
SYMIG	18	-2345.8789	4728.7786	0.0000	1.0000
K2PI	13	-2358.7914	4744.1221	0.0000	1.0000
K2PG	13	-2359.7421	4746.0234	0.0000	1.0000
K2PIG	14	-2358.7833	4746.1897	0.0000	1.0000
HKY	15	-2361.6281	4753.9694	0.0000	1.0000
SYM	16	-2371.4035	4775.6166	0.0000	1.0000
K2P	12	-2384.7593	4793.9802	0.0000	1.0000
F81I	15	-2382.5905	4795.8942	0.0000	1.0000
F81IG	16	-2382.5950	4797.9995	0.0000	1.0000
F81G	15	-2383.7094	4798.1321	0.0000	1.0000
JC69I	12	-2399.7721	4824.0057	0.0000	1.0000
JC69G	12	-2400.7863	4826.0342	0.0000	1.0000
JC69IG	13	-2399.7733	4826.0859	0.0000	1.0000
F81	14	-2406.4007	4841.4246	0.0000	1.0000
JC69	11	-2423.5987	4869.5873	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2313.2264	4666.4529	0.5380	0.5380
GTRG	20	-2314.0546	4668.1093	0.2350	0.7731
GTRIG	21	-2313.0897	4668.1793	0.2269	1.0000
HKYI	16	-2335.2681	4702.5361	0.0000	1.0000
HKYIG	17	-2335.2668	4704.5335	0.0000	1.0000
HKYG	16	-2336.3592	4704.7183	0.0000	1.0000
GTR	19	-2341.3260	4720.6520	0.0000	1.0000
SYMI	17	-2345.8952	4725.7903	0.0000	1.0000
SYMIG	18	-2345.8789	4727.7577	0.0000	1.0000
SYMG	17	-2346.9209	4727.8419	0.0000	1.0000
K2PI	13	-2358.7914	4743.5828	0.0000	1.0000
K2PG	13	-2359.7421	4745.4842	0.0000	1.0000
K2PIG	14	-2358.7833	4745.5666	0.0000	1.0000
HKY	15	-2361.6281	4753.2562	0.0000	1.0000
SYM	16	-2371.4035	4774.8071	0.0000	1.0000
K2P	12	-2384.7593	4793.5187	0.0000	1.0000
F81I	15	-2382.5905	4795.1810	0.0000	1.0000
F81IG	16	-2382.5950	4797.1899	0.0000	1.0000
F81G	15	-2383.7094	4797.4189	0.0000	1.0000
JC69I	12	-2399.7721	4823.5442	0.0000	1.0000
JC69IG	13	-2399.7733	4825.5466	0.0000	1.0000
JC69G	12	-2400.7863	4825.5726	0.0000	1.0000
F81	14	-2406.4007	4840.8014	0.0000	1.0000
JC69	11	-2423.5987	4869.1974	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	20	-2313.2264	4757.1577	0.6754	0.6754
GTRG	20	-2314.0546	4758.8141	0.2950	0.9704
GTRIG	21	-2313.0897	4763.4194	0.0295	0.9999
HKYI	16	-2335.2681	4775.1000	0.0001	1.0000
HKYG	16	-2336.3592	4777.2822	0.0000	1.0000
HKYIG	17	-2335.2668	4781.6326	0.0000	1.0000
K2PI	13	-2358.7914	4802.5410	0.0000	1.0000
SYMI	17	-2345.8952	4802.8894	0.0000	1.0000
K2PG	13	-2359.7421	4804.4423	0.0000	1.0000
SYMG	17	-2346.9209	4804.9410	0.0000	1.0000
GTR	19	-2341.3260	4806.8215	0.0000	1.0000
K2PIG	14	-2358.7833	4809.0599	0.0000	1.0000
SYMIG	18	-2345.8789	4809.3920	0.0000	1.0000
HKY	15	-2361.6281	4821.2848	0.0000	1.0000
SYM	16	-2371.4035	4847.3709	0.0000	1.0000
K2P	12	-2384.7593	4847.9416	0.0000	1.0000
F81I	15	-2382.5905	4863.2096	0.0000	1.0000
F81G	15	-2383.7094	4865.4475	0.0000	1.0000
F81IG	16	-2382.5950	4869.7538	0.0000	1.0000
JC69I	12	-2399.7721	4877.9671	0.0000	1.0000
JC69G	12	-2400.7863	4879.9955	0.0000	1.0000
JC69IG	13	-2399.7733	4884.5047	0.0000	1.0000
F81	14	-2406.4007	4904.2948	0.0000	1.0000
JC69	11	-2423.5987	4919.0850	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_13034.phylip = 1 - 689;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13034.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_13034.phylip = 1 - 689;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13034.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_13034.phylip = 1 - 689;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13034.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

