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

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

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


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
HKYI	16	-1642.9083	3318.8429	0.4084	0.4084
HKYG	16	-1643.1373	3319.3011	0.3248	0.7332
HKYIG	17	-1642.8435	3320.8438	0.1502	0.8834
GTRI	20	-1640.6287	3322.8544	0.0550	0.9384
GTRG	20	-1640.9019	3323.4007	0.0418	0.9802
GTRIG	21	-1640.5683	3324.8965	0.0198	1.0000
HKY	15	-1657.5318	3345.9675	0.0000	1.0000
GTR	19	-1655.5070	3350.4562	0.0000	1.0000
K2PI	13	-1668.7040	3364.0908	0.0000	1.0000
K2PG	13	-1668.7867	3364.2562	0.0000	1.0000
K2PIG	14	-1668.2986	3365.3867	0.0000	1.0000
SYMI	17	-1667.1164	3369.3897	0.0000	1.0000
SYMG	17	-1667.1474	3369.4518	0.0000	1.0000
SYMIG	18	-1667.0027	3371.3008	0.0000	1.0000
F81I	15	-1675.4070	3381.7180	0.0000	1.0000
F81G	15	-1675.5729	3382.0498	0.0000	1.0000
F81IG	16	-1675.3767	3383.7798	0.0000	1.0000
K2P	12	-1684.2487	3393.0816	0.0000	1.0000
SYM	16	-1683.0213	3399.0691	0.0000	1.0000
F81	14	-1688.9770	3406.7435	0.0000	1.0000
JC69G	12	-1698.3469	3421.2780	0.0000	1.0000
JC69I	12	-1698.3535	3421.2912	0.0000	1.0000
JC69IG	13	-1698.2584	3423.1997	0.0000	1.0000
JC69	11	-1712.3856	3447.2646	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1642.9083	3317.8165	0.3887	0.3887
HKYG	16	-1643.1373	3318.2747	0.3091	0.6978
HKYIG	17	-1642.8435	3319.6869	0.1526	0.8503
GTRI	20	-1640.6287	3321.2575	0.0696	0.9199
GTRG	20	-1640.9019	3321.8038	0.0529	0.9728
GTRIG	21	-1640.5683	3323.1365	0.0272	1.0000
HKY	15	-1657.5318	3345.0636	0.0000	1.0000
GTR	19	-1655.5070	3349.0140	0.0000	1.0000
K2PI	13	-1668.7040	3363.4079	0.0000	1.0000
K2PG	13	-1668.7867	3363.5733	0.0000	1.0000
K2PIG	14	-1668.2986	3364.5972	0.0000	1.0000
SYMI	17	-1667.1164	3368.2328	0.0000	1.0000
SYMG	17	-1667.1474	3368.2949	0.0000	1.0000
SYMIG	18	-1667.0027	3370.0053	0.0000	1.0000
F81I	15	-1675.4070	3380.8140	0.0000	1.0000
F81G	15	-1675.5729	3381.1458	0.0000	1.0000
F81IG	16	-1675.3767	3382.7534	0.0000	1.0000
K2P	12	-1684.2487	3392.4973	0.0000	1.0000
SYM	16	-1683.0213	3398.0427	0.0000	1.0000
F81	14	-1688.9770	3405.9540	0.0000	1.0000
JC69G	12	-1698.3469	3420.6937	0.0000	1.0000
JC69I	12	-1698.3535	3420.7070	0.0000	1.0000
JC69IG	13	-1698.2584	3422.5168	0.0000	1.0000
JC69	11	-1712.3856	3446.7712	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1642.9083	3386.6877	0.5432	0.5432
HKYG	16	-1643.1373	3387.1459	0.4320	0.9752
HKYIG	17	-1642.8435	3392.8625	0.0248	1.0000
GTRI	20	-1640.6287	3407.3464	0.0000	1.0000
GTRG	20	-1640.9019	3407.8928	0.0000	1.0000
HKY	15	-1657.5318	3409.6303	0.0000	1.0000
GTRIG	21	-1640.5683	3413.5300	0.0000	1.0000
K2PI	13	-1668.7040	3419.3658	0.0000	1.0000
K2PG	13	-1668.7867	3419.5312	0.0000	1.0000
K2PIG	14	-1668.2986	3424.8595	0.0000	1.0000
GTR	19	-1655.5070	3430.7986	0.0000	1.0000
SYMI	17	-1667.1164	3441.4084	0.0000	1.0000
SYMG	17	-1667.1474	3441.4705	0.0000	1.0000
K2P	12	-1684.2487	3444.1507	0.0000	1.0000
F81I	15	-1675.4070	3445.3808	0.0000	1.0000
F81G	15	-1675.5729	3445.7126	0.0000	1.0000
SYMIG	18	-1667.0027	3447.4854	0.0000	1.0000
F81IG	16	-1675.3767	3451.6246	0.0000	1.0000
F81	14	-1688.9770	3466.2163	0.0000	1.0000
SYM	16	-1683.0213	3466.9139	0.0000	1.0000
JC69G	12	-1698.3469	3472.3471	0.0000	1.0000
JC69I	12	-1698.3535	3472.3604	0.0000	1.0000
JC69IG	13	-1698.2584	3478.4746	0.0000	1.0000
JC69	11	-1712.3856	3494.1201	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 591-tomodel-phylip/chr2_9336.phylip = 1 - 547;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9336.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 591-tomodel-phylip/chr2_9336.phylip = 1 - 547;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9336.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 591-tomodel-phylip/chr2_9336.phylip = 1 - 547;
 Partition Dummy = 1:591-tomodel-phylip/chr2_9336.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

