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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	64	-2051.1970	4249.2603	0.6500	0.6500
GTRIG	65	-2050.5005	4250.5011	0.3495	0.9995
GTRI	64	-2058.3839	4263.6340	0.0005	1.0000
HKYG	60	-2072.2047	4280.8588	0.0000	1.0000
HKYIG	61	-2074.3992	4287.8343	0.0000	1.0000
HKYI	60	-2076.8519	4290.1533	0.0000	1.0000
SYMG	61	-2092.5740	4324.1840	0.0000	1.0000
K2PG	57	-2098.5661	4325.8912	0.0000	1.0000
SYMIG	62	-2094.9328	4331.4998	0.0000	1.0000
K2PIG	58	-2100.6855	4332.6819	0.0000	1.0000
SYMI	61	-2100.4422	4339.9204	0.0000	1.0000
K2PI	57	-2105.9183	4340.5954	0.0000	1.0000
GTR	63	-2108.4940	4361.2324	0.0000	1.0000
F81G	59	-2123.5978	4381.0700	0.0000	1.0000
F81I	59	-2126.9843	4387.8431	0.0000	1.0000
F81IG	60	-2126.3774	4389.2042	0.0000	1.0000
HKY	59	-2128.9010	4391.6765	0.0000	1.0000
JC69G	56	-2143.5861	4413.3905	0.0000	1.0000
JC69IG	57	-2146.6092	4421.9772	0.0000	1.0000
JC69I	56	-2150.6755	4427.5693	0.0000	1.0000
SYM	60	-2150.8022	4438.0538	0.0000	1.0000
K2P	56	-2159.1344	4444.4870	0.0000	1.0000
F81	58	-2180.0262	4491.3633	0.0000	1.0000
JC69	55	-2202.4052	4528.4992	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2051.1970	4230.3940	0.5751	0.5751
GTRIG	65	-2050.5005	4231.0011	0.4245	0.9996
GTRI	64	-2058.3839	4244.7677	0.0004	1.0000
HKYG	60	-2072.2047	4264.4094	0.0000	1.0000
HKYIG	61	-2074.3992	4270.7983	0.0000	1.0000
HKYI	60	-2076.8519	4273.7039	0.0000	1.0000
SYMG	61	-2092.5740	4307.1480	0.0000	1.0000
K2PG	57	-2098.5661	4311.1323	0.0000	1.0000
SYMIG	62	-2094.9328	4313.8655	0.0000	1.0000
K2PIG	58	-2100.6855	4317.3709	0.0000	1.0000
SYMI	61	-2100.4422	4322.8843	0.0000	1.0000
K2PI	57	-2105.9183	4325.8365	0.0000	1.0000
GTR	63	-2108.4940	4342.9881	0.0000	1.0000
F81G	59	-2123.5978	4365.1956	0.0000	1.0000
F81I	59	-2126.9843	4371.9687	0.0000	1.0000
F81IG	60	-2126.3774	4372.7548	0.0000	1.0000
HKY	59	-2128.9010	4375.8021	0.0000	1.0000
JC69G	56	-2143.5861	4399.1722	0.0000	1.0000
JC69IG	57	-2146.6092	4407.2183	0.0000	1.0000
JC69I	56	-2150.6755	4413.3511	0.0000	1.0000
SYM	60	-2150.8022	4421.6043	0.0000	1.0000
K2P	56	-2159.1344	4430.2688	0.0000	1.0000
F81	58	-2180.0262	4476.0523	0.0000	1.0000
JC69	55	-2202.4052	4514.8103	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2051.1970	4500.8924	0.9173	0.9173
GTRIG	65	-2050.5005	4505.7260	0.0818	0.9991
GTRI	64	-2058.3839	4515.2661	0.0007	0.9998
HKYG	60	-2072.2047	4518.0016	0.0002	1.0000
HKYI	60	-2076.8519	4527.2961	0.0000	1.0000
HKYIG	61	-2074.3992	4528.6170	0.0000	1.0000
K2PG	57	-2098.5661	4552.0449	0.0000	1.0000
K2PIG	58	-2100.6855	4562.5101	0.0000	1.0000
SYMG	61	-2092.5740	4564.9667	0.0000	1.0000
K2PI	57	-2105.9183	4566.7491	0.0000	1.0000
SYMIG	62	-2094.9328	4575.9108	0.0000	1.0000
SYMI	61	-2100.4422	4580.7031	0.0000	1.0000
GTR	63	-2108.4940	4609.2599	0.0000	1.0000
F81G	59	-2123.5978	4614.5612	0.0000	1.0000
F81I	59	-2126.9843	4621.3343	0.0000	1.0000
HKY	59	-2128.9010	4625.1677	0.0000	1.0000
F81IG	60	-2126.3774	4626.3470	0.0000	1.0000
JC69G	56	-2143.5861	4635.8583	0.0000	1.0000
JC69IG	57	-2146.6092	4648.1309	0.0000	1.0000
JC69I	56	-2150.6755	4650.0371	0.0000	1.0000
K2P	56	-2159.1344	4666.9548	0.0000	1.0000
SYM	60	-2150.8022	4675.1965	0.0000	1.0000
F81	58	-2180.0262	4721.1915	0.0000	1.0000
JC69	55	-2202.4052	4747.2698	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 183-phylip/chr3_454.phylip = 1 - 506;
 Partition Dummy = 1:183-phylip/chr3_454.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 183-phylip/chr3_454.phylip = 1 - 506;
 Partition Dummy = 1:183-phylip/chr3_454.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 (GTRG)]

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


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

