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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-3119.7331	6271.4663	0.6837	0.6837
GTRI	20	-3117.2345	6274.4690	0.1523	0.8360
HKYG	16	-3121.9400	6275.8800	0.0752	0.9113
GTRIG	21	-3117.2935	6276.5870	0.0528	0.9641
HKYIG	17	-3122.2634	6278.5269	0.0200	0.9841
GTRG	20	-3119.4952	6278.9904	0.0159	1.0000
HKY	15	-3149.2511	6328.5023	0.0000	1.0000
GTR	19	-3147.3227	6332.6454	0.0000	1.0000
F81I	15	-3177.4070	6384.8140	0.0000	1.0000
F81G	15	-3179.5509	6389.1018	0.0000	1.0000
F81IG	16	-3179.8423	6391.6845	0.0000	1.0000
K2PI	13	-3189.2089	6404.4178	0.0000	1.0000
SYMI	17	-3185.7486	6405.4972	0.0000	1.0000
K2PG	13	-3190.7668	6407.5336	0.0000	1.0000
SYMG	17	-3186.9432	6407.8864	0.0000	1.0000
K2PIG	14	-3190.7669	6409.5338	0.0000	1.0000
SYMIG	18	-3186.9433	6409.8865	0.0000	1.0000
F81	14	-3205.4071	6438.8143	0.0000	1.0000
K2P	12	-3218.3095	6460.6191	0.0000	1.0000
SYM	16	-3215.7325	6463.4651	0.0000	1.0000
JC69I	12	-3237.8459	6499.6917	0.0000	1.0000
JC69G	12	-3239.6058	6503.2116	0.0000	1.0000
JC69IG	13	-3239.8842	6505.7684	0.0000	1.0000
JC69	11	-3265.4937	6552.9873	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3119.7331	6351.9458	0.8989	0.8989
HKYG	16	-3121.9400	6356.3596	0.0989	0.9979
HKYIG	17	-3122.2634	6364.0364	0.0021	1.0000
GTRI	20	-3117.2345	6375.0684	0.0000	1.0000
GTRG	20	-3119.4952	6379.5898	0.0000	1.0000
GTRIG	21	-3117.2935	6382.2164	0.0000	1.0000
HKY	15	-3149.2511	6403.9519	0.0000	1.0000
GTR	19	-3147.3227	6428.2149	0.0000	1.0000
F81I	15	-3177.4070	6460.2636	0.0000	1.0000
F81G	15	-3179.5509	6464.5514	0.0000	1.0000
K2PI	13	-3189.2089	6469.8074	0.0000	1.0000
F81IG	16	-3179.8423	6472.1641	0.0000	1.0000
K2PG	13	-3190.7668	6472.9232	0.0000	1.0000
K2PIG	14	-3190.7669	6479.9535	0.0000	1.0000
SYMI	17	-3185.7486	6491.0068	0.0000	1.0000
SYMG	17	-3186.9432	6493.3959	0.0000	1.0000
SYMIG	18	-3186.9433	6500.4260	0.0000	1.0000
F81	14	-3205.4071	6509.2339	0.0000	1.0000
K2P	12	-3218.3095	6520.9788	0.0000	1.0000
SYM	16	-3215.7325	6543.9446	0.0000	1.0000
JC69I	12	-3237.8459	6560.0514	0.0000	1.0000
JC69G	12	-3239.6058	6563.5713	0.0000	1.0000
JC69IG	13	-3239.8842	6571.1581	0.0000	1.0000
JC69	11	-3265.4937	6608.3170	0.0000	1.0000



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

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

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


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

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


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

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

