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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2028.1910	4096.3820	0.4115	0.4115
GTRG	20	-2028.4821	4096.9642	0.3076	0.7191
GTRIG	21	-2028.1237	4098.2475	0.1619	0.8811
HKYI	16	-2034.1870	4100.3740	0.0559	0.9370
HKYG	16	-2034.4702	4100.9404	0.0421	0.9791
HKYIG	17	-2034.1872	4102.3745	0.0206	0.9997
GTR	19	-2036.5263	4111.0526	0.0003	0.9999
HKY	15	-2042.4300	4114.8599	0.0000	1.0000
F81I	15	-2044.3706	4118.7412	0.0000	1.0000
F81G	15	-2044.5184	4119.0367	0.0000	1.0000
F81IG	16	-2044.3710	4120.7419	0.0000	1.0000
F81	14	-2052.0244	4132.0488	0.0000	1.0000
SYMI	17	-2073.3132	4180.6263	0.0000	1.0000
SYMG	17	-2073.3960	4180.7920	0.0000	1.0000
K2PI	13	-2078.1902	4182.3804	0.0000	1.0000
SYMIG	18	-2073.1998	4182.3997	0.0000	1.0000
K2PG	13	-2078.3450	4182.6899	0.0000	1.0000
K2PIG	14	-2078.1904	4184.3808	0.0000	1.0000
SYM	16	-2080.6572	4193.3144	0.0000	1.0000
K2P	12	-2085.8635	4195.7270	0.0000	1.0000
JC69I	12	-2086.6536	4197.3072	0.0000	1.0000
JC69G	12	-2086.7582	4197.5164	0.0000	1.0000
JC69IG	13	-2086.6532	4199.3064	0.0000	1.0000
JC69	11	-2094.0639	4210.1278	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2034.1870	4180.9242	0.5571	0.5571
HKYG	16	-2034.4702	4181.4906	0.4197	0.9769
HKYIG	17	-2034.1872	4187.9591	0.0165	0.9934
HKY	15	-2042.4300	4190.3758	0.0049	0.9983
F81I	15	-2044.3706	4194.2570	0.0007	0.9990
F81G	15	-2044.5184	4194.5526	0.0006	0.9997
GTRI	20	-2028.1910	4197.0697	0.0002	0.9998
GTRG	20	-2028.4821	4197.6519	0.0001	1.0000
F81IG	16	-2044.3710	4201.2921	0.0000	1.0000
F81	14	-2052.0244	4202.5303	0.0000	1.0000
GTRIG	21	-2028.1237	4203.9696	0.0000	1.0000
GTR	19	-2036.5263	4206.7060	0.0000	1.0000
K2PI	13	-2078.1902	4247.8275	0.0000	1.0000
K2PG	13	-2078.3450	4248.1369	0.0000	1.0000
K2PIG	14	-2078.1904	4254.8622	0.0000	1.0000
K2P	12	-2085.8635	4256.1396	0.0000	1.0000
JC69I	12	-2086.6536	4257.7199	0.0000	1.0000
JC69G	12	-2086.7582	4257.9290	0.0000	1.0000
JC69IG	13	-2086.6532	4264.7534	0.0000	1.0000
JC69	11	-2094.0639	4265.5061	0.0000	1.0000
SYMI	17	-2073.3132	4266.2109	0.0000	1.0000
SYMG	17	-2073.3960	4266.3766	0.0000	1.0000
SYMIG	18	-2073.1998	4273.0187	0.0000	1.0000
SYM	16	-2080.6572	4273.8646	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/chr11_3210.phylip = 1 - 1135;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3210.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/chr11_3210.phylip = 1 - 1135;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3210.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3210.phylip = 1 - 1135;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3210.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

