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

Input data from file "591-tomodel-phylip/chr2_9409.phylip" (7 taxa, 979 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	-3053.8066	6147.6131	0.2906	0.2906
GTRG	20	-3053.9046	6147.8093	0.2634	0.5540
HKYI	16	-3058.5712	6149.1424	0.1353	0.6893
GTRIG	21	-3053.6335	6149.2671	0.1271	0.8164
HKYG	16	-3058.6447	6149.2894	0.1257	0.9421
HKYIG	17	-3058.4195	6150.8390	0.0579	1.0000
GTR	19	-3072.2148	6182.4296	0.0000	1.0000
HKY	15	-3076.4125	6182.8250	0.0000	1.0000
F81I	15	-3111.7749	6253.5498	0.0000	1.0000
F81G	15	-3111.8922	6253.7844	0.0000	1.0000
F81IG	16	-3111.6758	6255.3517	0.0000	1.0000
F81	14	-3127.9467	6283.8935	0.0000	1.0000
K2PG	13	-3158.8266	6343.6531	0.0000	1.0000
K2PI	13	-3159.0514	6344.1028	0.0000	1.0000
K2PIG	14	-3158.7924	6345.5847	0.0000	1.0000
SYMG	17	-3156.1928	6346.3856	0.0000	1.0000
SYMI	17	-3156.5965	6347.1929	0.0000	1.0000
SYMIG	18	-3156.1664	6348.3328	0.0000	1.0000
K2P	12	-3177.2062	6378.4124	0.0000	1.0000
SYM	16	-3174.9015	6381.8030	0.0000	1.0000
JC69G	12	-3196.5173	6417.0346	0.0000	1.0000
JC69I	12	-3196.5556	6417.1112	0.0000	1.0000
JC69IG	13	-3196.3994	6418.7988	0.0000	1.0000
JC69	11	-3213.0456	6448.0912	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3058.5712	6227.3269	0.5085	0.5085
HKYG	16	-3058.6447	6227.4739	0.4725	0.9810
HKYIG	17	-3058.4195	6233.9101	0.0189	0.9999
GTRI	20	-3053.8066	6245.3438	0.0001	0.9999
GTRG	20	-3053.9046	6245.5399	0.0001	1.0000
GTRIG	21	-3053.6335	6251.8842	0.0000	1.0000
HKY	15	-3076.4125	6256.1230	0.0000	1.0000
GTR	19	-3072.2148	6275.2737	0.0000	1.0000
F81I	15	-3111.7749	6326.8477	0.0000	1.0000
F81G	15	-3111.8922	6327.0824	0.0000	1.0000
F81IG	16	-3111.6758	6333.5362	0.0000	1.0000
F81	14	-3127.9467	6352.3049	0.0000	1.0000
K2PG	13	-3158.8266	6407.1780	0.0000	1.0000
K2PI	13	-3159.0514	6407.6277	0.0000	1.0000
K2PIG	14	-3158.7924	6413.9962	0.0000	1.0000
SYMG	17	-3156.1928	6429.4566	0.0000	1.0000
SYMI	17	-3156.5965	6430.2640	0.0000	1.0000
SYMIG	18	-3156.1664	6436.2904	0.0000	1.0000
K2P	12	-3177.2062	6437.0508	0.0000	1.0000
SYM	16	-3174.9015	6459.9875	0.0000	1.0000
JC69G	12	-3196.5173	6475.6729	0.0000	1.0000
JC69I	12	-3196.5556	6475.7496	0.0000	1.0000
JC69IG	13	-3196.3994	6482.3237	0.0000	1.0000
JC69	11	-3213.0456	6501.8430	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/chr2_9409.phylip = 1 - 979;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9409.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/chr2_9409.phylip = 1 - 979;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9409.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/chr2_9409.phylip = 1 - 979;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9409.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

