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

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

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


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	-2782.1581	5700.8847	0.3568	0.3568
GTRI	64	-2782.2593	5701.0870	0.3225	0.6793
GTRIG	65	-2782.1598	5703.1649	0.1141	0.7934
HKYG	60	-2788.0267	5703.5611	0.0936	0.8869
HKYI	60	-2788.1483	5703.8043	0.0829	0.9698
HKYIG	61	-2788.0295	5705.8249	0.0302	1.0000
K2PG	57	-2817.7181	5756.1970	0.0000	1.0000
K2PI	57	-2817.8261	5756.4130	0.0000	1.0000
K2PIG	58	-2817.7180	5758.4410	0.0000	1.0000
SYMG	61	-2814.4423	5758.6505	0.0000	1.0000
SYMI	61	-2814.5096	5758.7851	0.0000	1.0000
SYMIG	62	-2814.4429	5760.9145	0.0000	1.0000
GTR	63	-2815.7434	5765.7831	0.0000	1.0000
HKY	59	-2821.0098	5767.2738	0.0000	1.0000
F81I	59	-2838.5832	5802.4204	0.0000	1.0000
F81G	59	-2838.7531	5802.7603	0.0000	1.0000
F81IG	60	-2838.7528	5805.0133	0.0000	1.0000
K2P	56	-2848.3638	5815.2485	0.0000	1.0000
SYM	60	-2844.6917	5816.8911	0.0000	1.0000
JC69I	56	-2862.8467	5844.2143	0.0000	1.0000
JC69G	56	-2863.0230	5844.5669	0.0000	1.0000
JC69IG	57	-2863.0228	5846.8063	0.0000	1.0000
F81	58	-2868.8193	5860.6437	0.0000	1.0000
JC69	55	-2892.3338	5900.9534	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2782.1581	5692.3163	0.3819	0.3819
GTRI	64	-2782.2593	5692.5185	0.3451	0.7270
GTRIG	65	-2782.1598	5694.3196	0.1402	0.8672
HKYG	60	-2788.0267	5696.0534	0.0589	0.9262
HKYI	60	-2788.1483	5696.2966	0.0522	0.9784
HKYIG	61	-2788.0295	5698.0590	0.0216	1.0000
K2PG	57	-2817.7181	5749.4363	0.0000	1.0000
K2PI	57	-2817.8261	5749.6523	0.0000	1.0000
SYMG	61	-2814.4423	5750.8845	0.0000	1.0000
SYMI	61	-2814.5096	5751.0192	0.0000	1.0000
K2PIG	58	-2817.7180	5751.4359	0.0000	1.0000
SYMIG	62	-2814.4429	5752.8858	0.0000	1.0000
GTR	63	-2815.7434	5757.4868	0.0000	1.0000
HKY	59	-2821.0098	5760.0197	0.0000	1.0000
F81I	59	-2838.5832	5795.1663	0.0000	1.0000
F81G	59	-2838.7531	5795.5062	0.0000	1.0000
F81IG	60	-2838.7528	5797.5056	0.0000	1.0000
K2P	56	-2848.3638	5808.7276	0.0000	1.0000
SYM	60	-2844.6917	5809.3834	0.0000	1.0000
JC69I	56	-2862.8467	5837.6933	0.0000	1.0000
JC69G	56	-2863.0230	5838.0460	0.0000	1.0000
JC69IG	57	-2863.0228	5840.0456	0.0000	1.0000
F81	58	-2868.8193	5853.6386	0.0000	1.0000
JC69	55	-2892.3338	5894.6677	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2788.0267	5992.6407	0.5216	0.5216
HKYI	60	-2788.1483	5992.8839	0.4619	0.9835
HKYIG	61	-2788.0295	5999.5895	0.0162	0.9997
GTRG	64	-2782.1581	6008.6761	0.0002	0.9998
GTRI	64	-2782.2593	6008.8783	0.0002	1.0000
GTRIG	65	-2782.1598	6015.6225	0.0000	1.0000
K2PG	57	-2817.7181	6031.1942	0.0000	1.0000
K2PI	57	-2817.8261	6031.4102	0.0000	1.0000
K2PIG	58	-2817.7180	6038.1370	0.0000	1.0000
HKY	59	-2821.0098	6051.6639	0.0000	1.0000
SYMG	61	-2814.4423	6052.4150	0.0000	1.0000
SYMI	61	-2814.5096	6052.5497	0.0000	1.0000
SYMIG	62	-2814.4429	6059.3594	0.0000	1.0000
GTR	63	-2815.7434	6068.9035	0.0000	1.0000
K2P	56	-2848.3638	6085.5424	0.0000	1.0000
F81I	59	-2838.5832	6086.8105	0.0000	1.0000
F81G	59	-2838.7531	6087.1505	0.0000	1.0000
F81IG	60	-2838.7528	6094.0930	0.0000	1.0000
SYM	60	-2844.6917	6105.9707	0.0000	1.0000
JC69I	56	-2862.8467	6114.5082	0.0000	1.0000
JC69G	56	-2863.0230	6114.8609	0.0000	1.0000
JC69IG	57	-2863.0228	6121.8036	0.0000	1.0000
F81	58	-2868.8193	6140.3397	0.0000	1.0000
JC69	55	-2892.3338	6166.5394	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/chr11_3210.phylip = 1 - 1036;
 Partition Dummy = 1:183-phylip/chr11_3210.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/chr11_3210.phylip = 1 - 1036;
 Partition Dummy = 1:183-phylip/chr11_3210.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr11_3210.phylip = 1 - 1036;
 Partition Dummy = 1:183-phylip/chr11_3210.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

