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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3825.4898	7682.9796	0.5425	0.5425
HKYIG	17	-3825.5101	7685.0202	0.1956	0.7380
HKYI	16	-3827.0440	7686.0880	0.1147	0.8527
GTRG	20	-3823.2356	7686.4711	0.0947	0.9473
GTRIG	21	-3823.2453	7688.4906	0.0345	0.9818
GTRI	20	-3824.8859	7689.7718	0.0182	1.0000
HKY	15	-3850.3643	7730.7286	0.0000	1.0000
GTR	19	-3848.8113	7735.6227	0.0000	1.0000
K2PG	13	-3888.0194	7802.0387	0.0000	1.0000
SYMG	17	-3884.0524	7802.1049	0.0000	1.0000
K2PIG	14	-3888.0356	7804.0713	0.0000	1.0000
SYMIG	18	-3884.0681	7804.1363	0.0000	1.0000
K2PI	13	-3889.5759	7805.1518	0.0000	1.0000
SYMI	17	-3885.6458	7805.2916	0.0000	1.0000
F81G	15	-3905.2767	7840.5535	0.0000	1.0000
F81IG	16	-3905.2870	7842.5741	0.0000	1.0000
F81I	15	-3906.3881	7842.7762	0.0000	1.0000
K2P	12	-3912.4603	7848.9207	0.0000	1.0000
SYM	16	-3909.3082	7850.6164	0.0000	1.0000
F81	14	-3928.6851	7885.3702	0.0000	1.0000
JC69G	12	-3955.3713	7934.7427	0.0000	1.0000
JC69IG	13	-3955.4020	7936.8040	0.0000	1.0000
JC69I	12	-3956.5941	7937.1882	0.0000	1.0000
JC69	11	-3978.9233	7979.8466	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3825.4898	7763.2023	0.8060	0.8060
HKYI	16	-3827.0440	7766.3107	0.1703	0.9763
HKYIG	17	-3825.5101	7770.2568	0.0237	1.0000
GTRG	20	-3823.2356	7786.7494	0.0000	1.0000
GTRI	20	-3824.8859	7790.0501	0.0000	1.0000
GTRIG	21	-3823.2453	7793.7828	0.0000	1.0000
HKY	15	-3850.3643	7805.9373	0.0000	1.0000
GTR	19	-3848.8113	7830.8871	0.0000	1.0000
K2PG	13	-3888.0194	7867.2196	0.0000	1.0000
K2PI	13	-3889.5759	7870.3327	0.0000	1.0000
K2PIG	14	-3888.0356	7874.2661	0.0000	1.0000
SYMG	17	-3884.0524	7887.3414	0.0000	1.0000
SYMI	17	-3885.6458	7890.5281	0.0000	1.0000
SYMIG	18	-3884.0681	7894.3867	0.0000	1.0000
K2P	12	-3912.4603	7909.0877	0.0000	1.0000
F81G	15	-3905.2767	7915.7622	0.0000	1.0000
F81I	15	-3906.3881	7917.9850	0.0000	1.0000
F81IG	16	-3905.2870	7922.7967	0.0000	1.0000
SYM	16	-3909.3082	7930.8391	0.0000	1.0000
F81	14	-3928.6851	7955.5650	0.0000	1.0000
JC69G	12	-3955.3713	7994.9096	0.0000	1.0000
JC69I	12	-3956.5941	7997.3552	0.0000	1.0000
JC69IG	13	-3955.4020	8001.9849	0.0000	1.0000
JC69	11	-3978.9233	8034.9997	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32389.phylip = 1 - 1112;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32389.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32389.phylip = 1 - 1112;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32389.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32389.phylip = 1 - 1112;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32389.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

