
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:57 04/30/2011

Input data from file "591-tomodel-temp-phylip/chr10_4414.phylip" (7 taxa, 677 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	-2057.2323	4155.7451	0.3206	0.3206
HKYG	16	-2061.4925	4155.8093	0.3105	0.6311
GTRIG	21	-2057.2166	4157.8438	0.1123	0.7434
HKYIG	17	-2061.4767	4157.8820	0.1101	0.8535
HKYI	16	-2062.7402	4158.3047	0.0892	0.9427
GTRI	20	-2058.9544	4159.1893	0.0573	1.0000
GTR	19	-2089.7147	4218.5862	0.0000	1.0000
HKY	15	-2094.4536	4219.6334	0.0000	1.0000
F81G	15	-2102.2885	4235.3032	0.0000	1.0000
F81I	15	-2103.0180	4236.7622	0.0000	1.0000
F81IG	16	-2102.2674	4237.3590	0.0000	1.0000
K2PG	13	-2111.6049	4249.7588	0.0000	1.0000
K2PI	13	-2112.3867	4251.3223	0.0000	1.0000
K2PIG	14	-2111.5503	4251.7350	0.0000	1.0000
SYMG	17	-2110.6294	4256.1874	0.0000	1.0000
SYMIG	18	-2110.5229	4258.0854	0.0000	1.0000
SYMI	17	-2111.6460	4258.2206	0.0000	1.0000
F81	14	-2132.2137	4293.0618	0.0000	1.0000
K2P	12	-2143.4060	4311.2819	0.0000	1.0000
JC69G	12	-2144.2097	4312.8893	0.0000	1.0000
JC69I	12	-2144.7060	4313.8818	0.0000	1.0000
JC69IG	13	-2144.1481	4314.8451	0.0000	1.0000
SYM	16	-2142.4800	4317.7842	0.0000	1.0000
JC69	11	-2174.0010	4370.3989	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2057.2323	4154.4646	0.3530	0.3530
HKYG	16	-2061.4925	4154.9850	0.2721	0.6251
GTRIG	21	-2057.2166	4156.4331	0.1319	0.7571
HKYIG	17	-2061.4767	4156.9533	0.1017	0.8588
HKYI	16	-2062.7402	4157.4805	0.0781	0.9369
GTRI	20	-2058.9544	4157.9089	0.0631	1.0000
GTR	19	-2089.7147	4217.4294	0.0000	1.0000
HKY	15	-2094.4536	4218.9072	0.0000	1.0000
F81G	15	-2102.2885	4234.5770	0.0000	1.0000
F81I	15	-2103.0180	4236.0360	0.0000	1.0000
F81IG	16	-2102.2674	4236.5347	0.0000	1.0000
K2PG	13	-2111.6049	4249.2098	0.0000	1.0000
K2PI	13	-2112.3867	4250.7733	0.0000	1.0000
K2PIG	14	-2111.5503	4251.1006	0.0000	1.0000
SYMG	17	-2110.6294	4255.2587	0.0000	1.0000
SYMIG	18	-2110.5229	4257.0458	0.0000	1.0000
SYMI	17	-2111.6460	4257.2919	0.0000	1.0000
F81	14	-2132.2137	4292.4274	0.0000	1.0000
K2P	12	-2143.4060	4310.8120	0.0000	1.0000
JC69G	12	-2144.2097	4312.4194	0.0000	1.0000
JC69I	12	-2144.7060	4313.4119	0.0000	1.0000
JC69IG	13	-2144.1481	4314.2961	0.0000	1.0000
SYM	16	-2142.4800	4316.9600	0.0000	1.0000
JC69	11	-2174.0010	4370.0019	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2061.4925	4227.2678	0.7539	0.7539
HKYI	16	-2062.7402	4229.7632	0.2165	0.9704
HKYIG	17	-2061.4767	4233.7537	0.0294	0.9999
GTRG	20	-2057.2323	4244.8180	0.0001	1.0000
GTRI	20	-2058.9544	4248.2623	0.0000	1.0000
GTRIG	21	-2057.2166	4251.3042	0.0000	1.0000
HKY	15	-2094.4536	4286.6722	0.0000	1.0000
F81G	15	-2102.2885	4302.3421	0.0000	1.0000
GTR	19	-2089.7147	4303.2651	0.0000	1.0000
F81I	15	-2103.0180	4303.8011	0.0000	1.0000
K2PG	13	-2111.6049	4307.9395	0.0000	1.0000
F81IG	16	-2102.2674	4308.8175	0.0000	1.0000
K2PI	13	-2112.3867	4309.5030	0.0000	1.0000
K2PIG	14	-2111.5503	4314.3480	0.0000	1.0000
SYMG	17	-2110.6294	4332.0591	0.0000	1.0000
SYMI	17	-2111.6460	4334.0923	0.0000	1.0000
SYMIG	18	-2110.5229	4338.3639	0.0000	1.0000
F81	14	-2132.2137	4355.6748	0.0000	1.0000
K2P	12	-2143.4060	4365.0241	0.0000	1.0000
JC69G	12	-2144.2097	4366.6315	0.0000	1.0000
JC69I	12	-2144.7060	4367.6240	0.0000	1.0000
JC69IG	13	-2144.1481	4373.0258	0.0000	1.0000
SYM	16	-2142.4800	4389.2427	0.0000	1.0000
JC69	11	-2174.0010	4419.6963	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-temp-phylip/chr10_4414.phylip = 1 - 677;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4414.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-temp-phylip/chr10_4414.phylip = 1 - 677;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4414.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-temp-phylip/chr10_4414.phylip = 1 - 677;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4414.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

