
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:09 04/04/2011

Input data from file "683-tomodel-phylip/chr1_38361.phylip" (7 taxa, 1183 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-4390.2231	8820.4463	0.5722	0.5722
GTRIG	21	-4390.1004	8822.2008	0.2380	0.8102
GTRG	20	-4391.3269	8822.6537	0.1898	1.0000
SYMG	17	-4405.2368	8844.4737	0.0000	1.0000
SYMI	17	-4405.7357	8845.4714	0.0000	1.0000
SYMIG	18	-4405.0221	8846.0443	0.0000	1.0000
HKYG	16	-4408.5837	8849.1675	0.0000	1.0000
HKYI	16	-4409.1541	8850.3081	0.0000	1.0000
HKYIG	17	-4408.3688	8850.7377	0.0000	1.0000
K2PG	13	-4416.7857	8859.5715	0.0000	1.0000
K2PIG	14	-4416.6994	8861.3988	0.0000	1.0000
K2PI	13	-4417.7662	8861.5323	0.0000	1.0000
GTR	19	-4430.4541	8898.9083	0.0000	1.0000
SYM	16	-4445.0969	8922.1938	0.0000	1.0000
HKY	15	-4446.4923	8922.9845	0.0000	1.0000
K2P	12	-4455.6598	8935.3196	0.0000	1.0000
F81I	15	-4572.8928	9175.7856	0.0000	1.0000
F81IG	16	-4572.8876	9177.7752	0.0000	1.0000
F81G	15	-4574.2045	9178.4089	0.0000	1.0000
JC69I	12	-4584.7807	9193.5614	0.0000	1.0000
JC69G	12	-4585.7257	9195.4514	0.0000	1.0000
JC69IG	13	-4584.7369	9195.4739	0.0000	1.0000
F81	14	-4604.2142	9236.4285	0.0000	1.0000
JC69	11	-4615.9309	9253.8619	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	20	-4390.2231	8921.9624	0.6091	0.6091
GTRG	20	-4391.3269	8924.1699	0.2020	0.8111
K2PG	13	-4416.7857	8925.5570	0.1010	0.9121
K2PI	13	-4417.7662	8927.5178	0.0379	0.9500
GTRIG	21	-4390.1004	8928.7928	0.0200	0.9700
HKYG	16	-4408.5837	8930.3804	0.0091	0.9791
SYMG	17	-4405.2368	8930.7624	0.0075	0.9865
HKYI	16	-4409.1541	8931.5211	0.0051	0.9917
SYMI	17	-4405.7357	8931.7602	0.0045	0.9962
K2PIG	14	-4416.6994	8932.4601	0.0032	0.9994
HKYIG	17	-4408.3688	8937.0264	0.0003	0.9997
SYMIG	18	-4405.0221	8937.4089	0.0003	1.0000
GTR	19	-4430.4541	8995.3487	0.0000	1.0000
K2P	12	-4455.6598	8996.2293	0.0000	1.0000
HKY	15	-4446.4923	8999.1216	0.0000	1.0000
SYM	16	-4445.0969	9003.4067	0.0000	1.0000
F81I	15	-4572.8928	9251.9227	0.0000	1.0000
JC69I	12	-4584.7807	9254.4711	0.0000	1.0000
F81G	15	-4574.2045	9254.5461	0.0000	1.0000
JC69G	12	-4585.7257	9256.3611	0.0000	1.0000
F81IG	16	-4572.8876	9258.9881	0.0000	1.0000
JC69IG	13	-4584.7369	9261.4594	0.0000	1.0000
F81	14	-4604.2142	9307.4898	0.0000	1.0000
JC69	11	-4615.9309	9309.6958	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 683-tomodel-phylip/chr1_38361.phylip = 1 - 1183;
 Partition Dummy = 1:683-tomodel-phylip/chr1_38361.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 683-tomodel-phylip/chr1_38361.phylip = 1 - 1183;
 Partition Dummy = 1:683-tomodel-phylip/chr1_38361.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 (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_38361.phylip = 1 - 1183;
 Partition Dummy = 1:683-tomodel-phylip/chr1_38361.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;


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

