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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1504.7546	3042.6765	0.4028	0.4028
HKYI	16	-1504.9614	3043.0901	0.3275	0.7303
HKYIG	17	-1504.6148	3044.5457	0.1582	0.8884
GTRG	20	-1502.3792	3046.5767	0.0573	0.9457
GTRI	20	-1502.8919	3047.6020	0.0343	0.9800
GTRIG	21	-1502.3409	3048.6861	0.0200	1.0000
HKY	15	-1517.5150	3066.0578	0.0000	1.0000
GTR	19	-1515.4481	3070.5377	0.0000	1.0000
F81G	15	-1524.6257	3080.2792	0.0000	1.0000
F81I	15	-1524.7931	3080.6141	0.0000	1.0000
F81IG	16	-1524.4788	3082.1249	0.0000	1.0000
F81	14	-1536.6369	3102.1712	0.0000	1.0000
SYMI	17	-1553.6076	3142.5313	0.0000	1.0000
SYMG	17	-1554.1616	3143.6394	0.0000	1.0000
SYMIG	18	-1553.6093	3144.6928	0.0000	1.0000
K2PI	13	-1559.7775	3146.3312	0.0000	1.0000
K2PG	13	-1560.2958	3147.3677	0.0000	1.0000
K2PIG	14	-1559.7741	3148.4456	0.0000	1.0000
SYM	16	-1567.3969	3167.9611	0.0000	1.0000
K2P	12	-1572.2104	3169.0846	0.0000	1.0000
JC69I	12	-1572.8274	3170.3185	0.0000	1.0000
JC69G	12	-1573.3410	3171.3459	0.0000	1.0000
JC69IG	13	-1572.8244	3172.4248	0.0000	1.0000
JC69	11	-1584.4708	3191.5022	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1504.7546	3041.5091	0.3807	0.3807
HKYI	16	-1504.9614	3041.9227	0.3096	0.6903
HKYIG	17	-1504.6148	3043.2296	0.1611	0.8514
GTRG	20	-1502.3792	3044.7585	0.0750	0.9264
GTRI	20	-1502.8919	3045.7839	0.0449	0.9713
GTRIG	21	-1502.3409	3046.6818	0.0287	1.0000
HKY	15	-1517.5150	3065.0300	0.0000	1.0000
GTR	19	-1515.4481	3068.8962	0.0000	1.0000
F81G	15	-1524.6257	3079.2514	0.0000	1.0000
F81I	15	-1524.7931	3079.5863	0.0000	1.0000
F81IG	16	-1524.4788	3080.9575	0.0000	1.0000
F81	14	-1536.6369	3101.2738	0.0000	1.0000
SYMI	17	-1553.6076	3141.2151	0.0000	1.0000
SYMG	17	-1554.1616	3142.3233	0.0000	1.0000
SYMIG	18	-1553.6093	3143.2187	0.0000	1.0000
K2PI	13	-1559.7775	3145.5550	0.0000	1.0000
K2PG	13	-1560.2958	3146.5916	0.0000	1.0000
K2PIG	14	-1559.7741	3147.5482	0.0000	1.0000
SYM	16	-1567.3969	3166.7937	0.0000	1.0000
K2P	12	-1572.2104	3168.4207	0.0000	1.0000
JC69I	12	-1572.8274	3169.6547	0.0000	1.0000
JC69G	12	-1573.3410	3170.6821	0.0000	1.0000
JC69IG	13	-1572.8244	3171.6487	0.0000	1.0000
JC69	11	-1584.4708	3190.9417	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1504.7546	3108.3894	0.5360	0.5360
HKYI	16	-1504.9614	3108.8030	0.4359	0.9719
HKYIG	17	-1504.6148	3114.2899	0.0280	0.9999
HKY	15	-1517.5150	3127.7302	0.0000	1.0000
GTRG	20	-1502.3792	3128.3588	0.0000	1.0000
GTRI	20	-1502.8919	3129.3842	0.0000	1.0000
GTRIG	21	-1502.3409	3134.4621	0.0000	1.0000
F81G	15	-1524.6257	3141.9516	0.0000	1.0000
F81I	15	-1524.7931	3142.2865	0.0000	1.0000
F81IG	16	-1524.4788	3147.8378	0.0000	1.0000
GTR	19	-1515.4481	3148.3165	0.0000	1.0000
F81	14	-1536.6369	3159.7940	0.0000	1.0000
K2PI	13	-1559.7775	3199.8953	0.0000	1.0000
K2PG	13	-1560.2958	3200.9318	0.0000	1.0000
K2PIG	14	-1559.7741	3206.0684	0.0000	1.0000
SYMI	17	-1553.6076	3212.2754	0.0000	1.0000
SYMG	17	-1554.1616	3213.3836	0.0000	1.0000
SYMIG	18	-1553.6093	3218.4590	0.0000	1.0000
K2P	12	-1572.2104	3218.5809	0.0000	1.0000
JC69I	12	-1572.8274	3219.8149	0.0000	1.0000
JC69G	12	-1573.3410	3220.8423	0.0000	1.0000
JC69IG	13	-1572.8244	3225.9889	0.0000	1.0000
SYM	16	-1567.3969	3233.6740	0.0000	1.0000
JC69	11	-1584.4708	3236.9219	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/chr11_3450.phylip = 1 - 483;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3450.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/chr11_3450.phylip = 1 - 483;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3450.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/chr11_3450.phylip = 1 - 483;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3450.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

