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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3837.5650	7715.1300	0.7073	0.7073
GTRIG	21	-3837.5651	7717.1302	0.2602	0.9675
GTRI	20	-3840.9901	7721.9801	0.0230	0.9905
HKYG	16	-3846.2092	7724.4183	0.0068	0.9973
HKYIG	17	-3846.2094	7726.4187	0.0025	0.9998
HKYI	16	-3849.6440	7731.2880	0.0002	1.0000
GTR	19	-3855.2041	7748.4082	0.0000	1.0000
HKY	15	-3864.5959	7759.1918	0.0000	1.0000
K2PG	13	-3934.0716	7894.1431	0.0000	1.0000
K2PIG	14	-3934.0716	7896.1433	0.0000	1.0000
SYMG	17	-3933.3762	7900.7524	0.0000	1.0000
K2PI	13	-3937.7610	7901.5220	0.0000	1.0000
SYMIG	18	-3933.3763	7902.7526	0.0000	1.0000
SYMI	17	-3937.2012	7908.4023	0.0000	1.0000
K2P	12	-3958.8050	7941.6099	0.0000	1.0000
SYM	16	-3957.5743	7947.1486	0.0000	1.0000
F81G	15	-3984.4987	7998.9974	0.0000	1.0000
F81IG	16	-3984.4988	8000.9975	0.0000	1.0000
F81I	15	-3986.4507	8002.9014	0.0000	1.0000
F81	14	-4002.1687	8032.3375	0.0000	1.0000
JC69G	12	-4051.6154	8127.2307	0.0000	1.0000
JC69IG	13	-4051.6154	8129.2308	0.0000	1.0000
JC69I	12	-4053.6580	8131.3160	0.0000	1.0000
JC69	11	-4072.3294	8166.6588	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3846.2092	7804.3066	0.9367	0.9367
HKYI	16	-3849.6440	7811.1763	0.0302	0.9669
HKYIG	17	-3846.2094	7811.3000	0.0284	0.9952
GTRG	20	-3837.5650	7814.9903	0.0045	0.9997
GTRI	20	-3840.9901	7821.8404	0.0001	0.9999
GTRIG	21	-3837.5651	7821.9835	0.0001	1.0000
HKY	15	-3864.5959	7834.0870	0.0000	1.0000
GTR	19	-3855.2041	7843.2755	0.0000	1.0000
K2PG	13	-3934.0716	7959.0523	0.0000	1.0000
K2PIG	14	-3934.0716	7966.0455	0.0000	1.0000
K2PI	13	-3937.7610	7966.4312	0.0000	1.0000
SYMG	17	-3933.3762	7985.6337	0.0000	1.0000
SYMIG	18	-3933.3763	7992.6269	0.0000	1.0000
SYMI	17	-3937.2012	7993.2836	0.0000	1.0000
K2P	12	-3958.8050	8001.5261	0.0000	1.0000
SYM	16	-3957.5743	8027.0369	0.0000	1.0000
F81G	15	-3984.4987	8073.8926	0.0000	1.0000
F81I	15	-3986.4507	8077.7966	0.0000	1.0000
F81IG	16	-3984.4988	8080.8857	0.0000	1.0000
F81	14	-4002.1687	8102.2397	0.0000	1.0000
JC69G	12	-4051.6154	8187.1469	0.0000	1.0000
JC69I	12	-4053.6580	8191.2322	0.0000	1.0000
JC69IG	13	-4051.6154	8194.1400	0.0000	1.0000
JC69	11	-4072.3294	8221.5819	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 591-tomodel-phylip/chr17_3015.phylip = 1 - 1089;
 Partition Dummy = 1:591-tomodel-phylip/chr17_3015.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 591-tomodel-phylip/chr17_3015.phylip = 1 - 1089;
 Partition Dummy = 1:591-tomodel-phylip/chr17_3015.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 591-tomodel-phylip/chr17_3015.phylip = 1 - 1089;
 Partition Dummy = 1:591-tomodel-phylip/chr17_3015.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

