
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:37 02/17/2011

Input data from file "183-phylip/chr1_31654.phylip" (29 taxa, 567 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1922.8810	3980.2285	0.5501	0.5501
HKYIG	61	-1922.0964	3981.1711	0.3434	0.8935
GTRG	64	-1920.0819	3984.7374	0.0577	0.9513
GTRIG	65	-1919.2583	3985.6423	0.0367	0.9880
HKYI	60	-1926.7722	3988.0109	0.0112	0.9992
GTRI	64	-1924.3898	3993.3534	0.0008	1.0000
K2PIG	58	-1938.3896	4006.2517	0.0000	1.0000
K2PG	57	-1939.9666	4006.9234	0.0000	1.0000
K2PI	57	-1942.6473	4012.2849	0.0000	1.0000
SYMIG	62	-1937.9585	4015.4170	0.0000	1.0000
SYMG	61	-1939.5768	4016.1317	0.0000	1.0000
SYMI	61	-1942.3953	4021.7688	0.0000	1.0000
HKY	59	-1968.6304	4069.2253	0.0000	1.0000
F81IG	60	-1967.9585	4070.3834	0.0000	1.0000
F81G	59	-1969.3580	4070.6805	0.0000	1.0000
GTR	63	-1965.7836	4073.5989	0.0000	1.0000
F81I	59	-1971.9169	4075.7982	0.0000	1.0000
JC69IG	57	-1981.5736	4090.1373	0.0000	1.0000
JC69G	56	-1982.9755	4090.4686	0.0000	1.0000
JC69I	56	-1985.5386	4095.5948	0.0000	1.0000
K2P	56	-1985.6297	4095.7771	0.0000	1.0000
SYM	60	-1984.3161	4103.0985	0.0000	1.0000
F81	58	-2013.3300	4156.1325	0.0000	1.0000
JC69	55	-2027.2644	4176.5836	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1922.8810	3965.7621	0.4195	0.4195
HKYIG	61	-1922.0964	3966.1929	0.3382	0.7577
GTRG	64	-1920.0819	3968.1637	0.1262	0.8839
GTRIG	65	-1919.2583	3968.5166	0.1058	0.9897
HKYI	60	-1926.7722	3973.5445	0.0086	0.9983
GTRI	64	-1924.3898	3976.7797	0.0017	1.0000
K2PIG	58	-1938.3896	3992.7792	0.0000	1.0000
K2PG	57	-1939.9666	3993.9332	0.0000	1.0000
K2PI	57	-1942.6473	3999.2947	0.0000	1.0000
SYMIG	62	-1937.9585	3999.9170	0.0000	1.0000
SYMG	61	-1939.5768	4001.1535	0.0000	1.0000
SYMI	61	-1942.3953	4006.7906	0.0000	1.0000
HKY	59	-1968.6304	4055.2608	0.0000	1.0000
F81IG	60	-1967.9585	4055.9170	0.0000	1.0000
F81G	59	-1969.3580	4056.7160	0.0000	1.0000
GTR	63	-1965.7836	4057.5671	0.0000	1.0000
F81I	59	-1971.9169	4061.8337	0.0000	1.0000
JC69IG	57	-1981.5736	4077.1472	0.0000	1.0000
JC69G	56	-1982.9755	4077.9510	0.0000	1.0000
JC69I	56	-1985.5386	4083.0772	0.0000	1.0000
K2P	56	-1985.6297	4083.2595	0.0000	1.0000
SYM	60	-1984.3161	4088.6321	0.0000	1.0000
F81	58	-2013.3300	4142.6601	0.0000	1.0000
JC69	55	-2027.2644	4164.5288	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1922.8810	4226.1836	0.8983	0.8983
HKYIG	61	-1922.0964	4230.9548	0.0827	0.9810
HKYI	60	-1926.7722	4233.9661	0.0183	0.9994
K2PG	57	-1939.9666	4241.3337	0.0005	0.9998
K2PIG	58	-1938.3896	4244.5201	0.0001	0.9999
GTRG	64	-1920.0819	4245.9467	0.0000	1.0000
K2PI	57	-1942.6473	4246.6952	0.0000	1.0000
GTRIG	65	-1919.2583	4250.6399	0.0000	1.0000
GTRI	64	-1924.3898	4254.5627	0.0000	1.0000
SYMG	61	-1939.5768	4265.9154	0.0000	1.0000
SYMIG	62	-1937.9585	4269.0193	0.0000	1.0000
SYMI	61	-1942.3953	4271.5525	0.0000	1.0000
HKY	59	-1968.6304	4311.3420	0.0000	1.0000
F81G	59	-1969.3580	4312.7972	0.0000	1.0000
F81IG	60	-1967.9585	4316.3385	0.0000	1.0000
F81I	59	-1971.9169	4317.9149	0.0000	1.0000
JC69G	56	-1982.9755	4321.0111	0.0000	1.0000
JC69IG	57	-1981.5736	4324.5476	0.0000	1.0000
JC69I	56	-1985.5386	4326.1373	0.0000	1.0000
K2P	56	-1985.6297	4326.3196	0.0000	1.0000
GTR	63	-1965.7836	4331.0097	0.0000	1.0000
SYM	60	-1984.3161	4349.0537	0.0000	1.0000
F81	58	-2013.3300	4394.4009	0.0000	1.0000
JC69	55	-2027.2644	4403.2486	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_31654.phylip = 1 - 567;
 Partition Dummy = 1:183-phylip/chr1_31654.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;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_31654.phylip = 1 - 567;
 Partition Dummy = 1:183-phylip/chr1_31654.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;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_31654.phylip = 1 - 567;
 Partition Dummy = 1:183-phylip/chr1_31654.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

