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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-2041.4664	4231.2271	0.4841	0.4841
HKYIG	61	-2046.9771	4231.9457	0.3380	0.8221
GTRG	64	-2044.1796	4234.0614	0.1174	0.9394
HKYG	60	-2049.9708	4235.3847	0.0606	1.0000
HKYI	60	-2061.2550	4257.9531	0.0000	1.0000
GTRI	64	-2056.4545	4258.6111	0.0000	1.0000
K2PIG	58	-2065.3746	4261.1273	0.0000	1.0000
K2PG	57	-2068.0795	4264.0207	0.0000	1.0000
SYMIG	62	-2062.2748	4265.1005	0.0000	1.0000
SYMG	61	-2064.2637	4266.5190	0.0000	1.0000
K2PI	57	-2075.9328	4279.7272	0.0000	1.0000
SYMI	61	-2073.2370	4284.4655	0.0000	1.0000
F81G	59	-2101.4824	4335.8701	0.0000	1.0000
GTR	63	-2096.4224	4335.9658	0.0000	1.0000
HKY	59	-2101.6929	4336.2911	0.0000	1.0000
F81IG	60	-2102.2028	4339.8486	0.0000	1.0000
JC69IG	57	-2112.7803	4353.4222	0.0000	1.0000
JC69G	56	-2115.1297	4355.6151	0.0000	1.0000
F81I	59	-2113.1370	4359.1794	0.0000	1.0000
K2P	56	-2119.8817	4365.1191	0.0000	1.0000
SYM	60	-2117.9421	4371.3272	0.0000	1.0000
JC69I	56	-2123.3052	4371.9660	0.0000	1.0000
F81	58	-2150.9627	4432.3036	0.0000	1.0000
JC69	55	-2165.2057	4453.2715	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2041.4664	4212.9329	0.6987	0.6987
HKYIG	61	-2046.9771	4215.9541	0.1543	0.8530
GTRG	64	-2044.1796	4216.3592	0.1260	0.9790
HKYG	60	-2049.9708	4219.9417	0.0210	1.0000
GTRI	64	-2056.4545	4240.9090	0.0000	1.0000
HKYI	60	-2061.2550	4242.5101	0.0000	1.0000
K2PIG	58	-2065.3746	4246.7491	0.0000	1.0000
SYMIG	62	-2062.2748	4248.5497	0.0000	1.0000
K2PG	57	-2068.0795	4250.1590	0.0000	1.0000
SYMG	61	-2064.2637	4250.5275	0.0000	1.0000
K2PI	57	-2075.9328	4265.8655	0.0000	1.0000
SYMI	61	-2073.2370	4268.4739	0.0000	1.0000
GTR	63	-2096.4224	4318.8448	0.0000	1.0000
F81G	59	-2101.4824	4320.9648	0.0000	1.0000
HKY	59	-2101.6929	4321.3859	0.0000	1.0000
F81IG	60	-2102.2028	4324.4056	0.0000	1.0000
JC69IG	57	-2112.7803	4339.5606	0.0000	1.0000
JC69G	56	-2115.1297	4342.2595	0.0000	1.0000
F81I	59	-2113.1370	4344.2741	0.0000	1.0000
K2P	56	-2119.8817	4351.7634	0.0000	1.0000
SYM	60	-2117.9421	4355.8842	0.0000	1.0000
JC69I	56	-2123.3052	4358.6104	0.0000	1.0000
F81	58	-2150.9627	4417.9255	0.0000	1.0000
JC69	55	-2165.2057	4440.4114	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2049.9708	4476.8777	0.5361	0.5361
HKYIG	61	-2046.9771	4477.1724	0.4627	0.9988
GTRG	64	-2044.1796	4490.4243	0.0006	0.9994
GTRIG	65	-2041.4664	4491.2802	0.0004	0.9998
K2PG	57	-2068.0795	4494.2482	0.0001	0.9999
K2PIG	58	-2065.3746	4495.1206	0.0001	1.0000
HKYI	60	-2061.2550	4499.4461	0.0000	1.0000
K2PI	57	-2075.9328	4509.9547	0.0000	1.0000
SYMG	61	-2064.2637	4511.7458	0.0000	1.0000
SYMIG	62	-2062.2748	4514.0502	0.0000	1.0000
GTRI	64	-2056.4545	4514.9740	0.0000	1.0000
SYMI	61	-2073.2370	4529.6922	0.0000	1.0000
F81G	59	-2101.4824	4573.6186	0.0000	1.0000
HKY	59	-2101.6929	4574.0396	0.0000	1.0000
F81IG	60	-2102.2028	4581.3416	0.0000	1.0000
JC69G	56	-2115.1297	4582.0664	0.0000	1.0000
JC69IG	57	-2112.7803	4583.6498	0.0000	1.0000
GTR	63	-2096.4224	4588.6276	0.0000	1.0000
K2P	56	-2119.8817	4591.5704	0.0000	1.0000
F81I	59	-2113.1370	4596.9278	0.0000	1.0000
JC69I	56	-2123.3052	4598.4173	0.0000	1.0000
SYM	60	-2117.9421	4612.8202	0.0000	1.0000
F81	58	-2150.9627	4666.2969	0.0000	1.0000
JC69	55	-2165.2057	4675.9361	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_5441.phylip = 1 - 535;
 Partition Dummy = 1:183-phylip/chr3_5441.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_5441.phylip = 1 - 535;
 Partition Dummy = 1:183-phylip/chr3_5441.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

