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

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

Minimum AIC  model: GTR
Minimum AICc model: GTR
Minimum BIC  model: HKY


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
GTR	19	-1580.2162	3199.4953	0.3687	0.3687
GTRG	20	-1579.4523	3200.0810	0.2751	0.6438
GTRI	20	-1579.7262	3200.6289	0.2092	0.8530
GTRIG	21	-1579.3540	3202.0039	0.1052	0.9582
HKY	15	-1587.5314	3205.7304	0.0163	0.9745
HKYG	16	-1586.8042	3206.3660	0.0119	0.9864
HKYI	16	-1587.0267	3206.8110	0.0095	0.9959
HKYIG	17	-1586.8163	3208.4861	0.0041	1.0000
F81	14	-1604.6985	3237.9803	0.0000	1.0000
F81G	15	-1603.9751	3238.6178	0.0000	1.0000
F81I	15	-1604.1445	3238.9565	0.0000	1.0000
F81IG	16	-1603.9893	3240.7363	0.0000	1.0000
K2PG	13	-1621.0600	3268.6249	0.0000	1.0000
K2PI	13	-1621.4998	3269.5045	0.0000	1.0000
SYMG	17	-1617.6114	3270.0763	0.0000	1.0000
K2P	12	-1622.9144	3270.2609	0.0000	1.0000
K2PIG	14	-1621.0474	3270.6782	0.0000	1.0000
SYMI	17	-1618.1059	3271.0654	0.0000	1.0000
SYM	16	-1619.4539	3271.6655	0.0000	1.0000
SYMIG	18	-1617.6063	3272.1679	0.0000	1.0000
JC69G	12	-1636.6193	3297.6707	0.0000	1.0000
JC69	11	-1637.7729	3297.9110	0.0000	1.0000
JC69I	12	-1636.8011	3298.0342	0.0000	1.0000
JC69IG	13	-1636.6219	3299.7487	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	19	-1580.2162	3198.4323	0.3563	0.3563
GTRG	20	-1579.4523	3198.9045	0.2814	0.6376
GTRI	20	-1579.7262	3199.4525	0.2139	0.8515
GTRIG	21	-1579.3540	3200.7080	0.1142	0.9657
HKY	15	-1587.5314	3205.0628	0.0129	0.9787
HKYG	16	-1586.8042	3205.6084	0.0099	0.9885
HKYI	16	-1587.0267	3206.0533	0.0079	0.9964
HKYIG	17	-1586.8163	3207.6326	0.0036	1.0000
F81	14	-1604.6985	3237.3970	0.0000	1.0000
F81G	15	-1603.9751	3237.9502	0.0000	1.0000
F81I	15	-1604.1445	3238.2889	0.0000	1.0000
F81IG	16	-1603.9893	3239.9786	0.0000	1.0000
K2PG	13	-1621.0600	3268.1200	0.0000	1.0000
K2PI	13	-1621.4998	3268.9996	0.0000	1.0000
SYMG	17	-1617.6114	3269.2227	0.0000	1.0000
K2P	12	-1622.9144	3269.8288	0.0000	1.0000
K2PIG	14	-1621.0474	3270.0949	0.0000	1.0000
SYMI	17	-1618.1059	3270.2119	0.0000	1.0000
SYM	16	-1619.4539	3270.9079	0.0000	1.0000
SYMIG	18	-1617.6063	3271.2126	0.0000	1.0000
JC69G	12	-1636.6193	3297.2386	0.0000	1.0000
JC69	11	-1637.7729	3297.5458	0.0000	1.0000
JC69I	12	-1636.8011	3297.6021	0.0000	1.0000
JC69IG	13	-1636.6219	3299.2438	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1587.5314	3274.0609	0.8746	0.8746
HKYG	16	-1586.8042	3279.2063	0.0668	0.9413
HKYI	16	-1587.0267	3279.6512	0.0534	0.9948
GTR	19	-1580.2162	3285.8299	0.0024	0.9972
HKYIG	17	-1586.8163	3285.8304	0.0024	0.9997
GTRG	20	-1579.4523	3290.9019	0.0002	0.9998
GTRI	20	-1579.7262	3291.4499	0.0001	1.0000
GTRIG	21	-1579.3540	3297.3053	0.0000	1.0000
F81	14	-1604.6985	3301.7952	0.0000	1.0000
F81G	15	-1603.9751	3306.9483	0.0000	1.0000
F81I	15	-1604.1445	3307.2870	0.0000	1.0000
F81IG	16	-1603.9893	3313.5765	0.0000	1.0000
K2P	12	-1622.9144	3325.0272	0.0000	1.0000
K2PG	13	-1621.0600	3327.9183	0.0000	1.0000
K2PI	13	-1621.4998	3328.7979	0.0000	1.0000
K2PIG	14	-1621.0474	3334.4930	0.0000	1.0000
SYM	16	-1619.4539	3344.5058	0.0000	1.0000
SYMG	17	-1617.6114	3347.4205	0.0000	1.0000
JC69	11	-1637.7729	3348.1444	0.0000	1.0000
SYMI	17	-1618.1059	3348.4097	0.0000	1.0000
JC69G	12	-1636.6193	3352.4370	0.0000	1.0000
JC69I	12	-1636.8011	3352.8005	0.0000	1.0000
SYMIG	18	-1617.6063	3354.0103	0.0000	1.0000
JC69IG	13	-1636.6219	3359.0421	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_8064.phylip = 1 - 735;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8064.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_8064.phylip = 1 - 735;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8064.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_8064.phylip = 1 - 735;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8064.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

