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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1684.2872	3409.9447	0.2774	0.2774
GTRG	20	-1684.5050	3410.3804	0.2231	0.5004
HKYI	16	-1688.8607	3410.6031	0.1996	0.7000
HKYG	16	-1689.0462	3410.9740	0.1658	0.8658
GTRIG	21	-1684.5051	3412.5200	0.0765	0.9423
HKYIG	17	-1689.0462	3413.0859	0.0577	1.0000
K2PI	13	-1702.2282	3431.0434	0.0000	1.0000
K2PG	13	-1702.2984	3431.1839	0.0000	1.0000
SYMI	17	-1698.8870	3432.7674	0.0000	1.0000
SYMG	17	-1698.9718	3432.9370	0.0000	1.0000
K2PIG	14	-1702.2984	3433.2753	0.0000	1.0000
SYMIG	18	-1698.9717	3435.0557	0.0000	1.0000
GTR	19	-1709.6178	3458.4734	0.0000	1.0000
HKY	15	-1715.2954	3461.3675	0.0000	1.0000
SYM	16	-1723.5074	3479.8965	0.0000	1.0000
K2P	12	-1728.2797	3481.0618	0.0000	1.0000
F81G	15	-1737.0020	3504.7807	0.0000	1.0000
F81I	15	-1737.2150	3505.2067	0.0000	1.0000
F81IG	16	-1736.9998	3506.8814	0.0000	1.0000
JC69G	12	-1747.1490	3518.8005	0.0000	1.0000
JC69I	12	-1747.3865	3519.2754	0.0000	1.0000
JC69IG	13	-1747.1336	3520.8544	0.0000	1.0000
F81	14	-1760.1508	3548.9801	0.0000	1.0000
JC69	11	-1770.3607	3563.1458	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1684.2872	3408.5744	0.3027	0.3027
GTRG	20	-1684.5050	3409.0101	0.2434	0.5461
HKYI	16	-1688.8607	3409.7214	0.1706	0.7166
HKYG	16	-1689.0462	3410.0923	0.1417	0.8583
GTRIG	21	-1684.5051	3411.0102	0.0895	0.9479
HKYIG	17	-1689.0462	3412.0924	0.0521	1.0000
K2PI	13	-1702.2282	3430.4563	0.0000	1.0000
K2PG	13	-1702.2984	3430.5968	0.0000	1.0000
SYMI	17	-1698.8870	3431.7739	0.0000	1.0000
SYMG	17	-1698.9718	3431.9435	0.0000	1.0000
K2PIG	14	-1702.2984	3432.5968	0.0000	1.0000
SYMIG	18	-1698.9717	3433.9435	0.0000	1.0000
GTR	19	-1709.6178	3457.2356	0.0000	1.0000
HKY	15	-1715.2954	3460.5908	0.0000	1.0000
SYM	16	-1723.5074	3479.0148	0.0000	1.0000
K2P	12	-1728.2797	3480.5594	0.0000	1.0000
F81G	15	-1737.0020	3504.0040	0.0000	1.0000
F81I	15	-1737.2150	3504.4300	0.0000	1.0000
F81IG	16	-1736.9998	3505.9997	0.0000	1.0000
JC69G	12	-1747.1490	3518.2981	0.0000	1.0000
JC69I	12	-1747.3865	3518.7729	0.0000	1.0000
JC69IG	13	-1747.1336	3520.2673	0.0000	1.0000
F81	14	-1760.1508	3548.3016	0.0000	1.0000
JC69	11	-1770.3607	3562.7214	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1688.8607	3480.9542	0.5226	0.5226
HKYG	16	-1689.0462	3481.3251	0.4342	0.9568
HKYIG	17	-1689.0462	3487.7772	0.0172	0.9741
K2PI	13	-1702.2282	3488.3330	0.0131	0.9871
K2PG	13	-1702.2984	3488.4734	0.0122	0.9993
K2PIG	14	-1702.2984	3494.9255	0.0005	0.9998
GTRI	20	-1684.2872	3497.6154	0.0001	0.9999
GTRG	20	-1684.5050	3498.0510	0.0001	1.0000
GTRIG	21	-1684.5051	3504.5033	0.0000	1.0000
SYMI	17	-1698.8870	3507.4587	0.0000	1.0000
SYMG	17	-1698.9718	3507.6284	0.0000	1.0000
SYMIG	18	-1698.9717	3514.0804	0.0000	1.0000
HKY	15	-1715.2954	3527.3716	0.0000	1.0000
K2P	12	-1728.2797	3533.9840	0.0000	1.0000
GTR	19	-1709.6178	3541.8245	0.0000	1.0000
SYM	16	-1723.5074	3550.2476	0.0000	1.0000
F81G	15	-1737.0020	3570.7848	0.0000	1.0000
F81I	15	-1737.2150	3571.2108	0.0000	1.0000
JC69G	12	-1747.1490	3571.7227	0.0000	1.0000
JC69I	12	-1747.3865	3572.1975	0.0000	1.0000
F81IG	16	-1736.9998	3577.2325	0.0000	1.0000
JC69IG	13	-1747.1336	3578.1439	0.0000	1.0000
F81	14	-1760.1508	3610.6303	0.0000	1.0000
JC69	11	-1770.3607	3611.6939	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10216.phylip = 1 - 634;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10216.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10216.phylip = 1 - 634;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10216.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

