
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:10 02/24/2011

Input data from file "/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_19948.phylip" (19 taxa, 371 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	40	-1515.0949	3120.1293	0.7033	0.7033
HKYIG	41	-1515.0915	3122.6512	0.1993	0.9026
HKYG	40	-1517.1588	3124.2570	0.0893	0.9919
GTRI	44	-1514.7652	3129.6777	0.0059	0.9978
GTRIG	45	-1514.7702	3132.2789	0.0016	0.9994
GTRG	44	-1517.0697	3134.2866	0.0006	1.0000
F81I	39	-1537.2481	3161.9222	0.0000	1.0000
F81IG	40	-1537.2283	3164.3960	0.0000	1.0000
F81G	39	-1538.6258	3164.6777	0.0000	1.0000
K2PI	37	-1543.9564	3170.3572	0.0000	1.0000
K2PG	37	-1544.0838	3170.6120	0.0000	1.0000
K2PIG	38	-1543.6341	3172.1959	0.0000	1.0000
SYMG	41	-1540.2034	3172.8749	0.0000	1.0000
SYMI	41	-1541.2312	3174.9305	0.0000	1.0000
SYMIG	42	-1540.2033	3175.4188	0.0000	1.0000
HKY	39	-1554.5032	3196.4324	0.0000	1.0000
JC69I	36	-1560.3546	3200.6851	0.0000	1.0000
JC69G	36	-1560.7889	3201.5539	0.0000	1.0000
JC69IG	37	-1560.1702	3202.7849	0.0000	1.0000
GTR	43	-1553.9164	3205.4046	0.0000	1.0000
F81	38	-1573.4961	3231.9200	0.0000	1.0000
K2P	36	-1579.4043	3238.7846	0.0000	1.0000
SYM	40	-1575.0858	3240.1111	0.0000	1.0000
JC69	35	-1594.9515	3267.4254	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1515.0949	3110.1899	0.6521	0.6521
HKYIG	41	-1515.0915	3112.1831	0.2407	0.8929
HKYG	40	-1517.1588	3114.3176	0.0828	0.9757
GTRI	44	-1514.7652	3117.5304	0.0166	0.9923
GTRIG	45	-1514.7702	3119.5404	0.0061	0.9983
GTRG	44	-1517.0697	3122.1393	0.0017	1.0000
F81I	39	-1537.2481	3152.4963	0.0000	1.0000
F81IG	40	-1537.2283	3154.4566	0.0000	1.0000
F81G	39	-1538.6258	3155.2517	0.0000	1.0000
K2PI	37	-1543.9564	3161.9127	0.0000	1.0000
K2PG	37	-1544.0838	3162.1676	0.0000	1.0000
SYMG	41	-1540.2034	3162.4068	0.0000	1.0000
K2PIG	38	-1543.6341	3163.2682	0.0000	1.0000
SYMIG	42	-1540.2033	3164.4066	0.0000	1.0000
SYMI	41	-1541.2312	3164.4624	0.0000	1.0000
HKY	39	-1554.5032	3187.0064	0.0000	1.0000
JC69I	36	-1560.3546	3192.7091	0.0000	1.0000
JC69G	36	-1560.7889	3193.5778	0.0000	1.0000
GTR	43	-1553.9164	3193.8327	0.0000	1.0000
JC69IG	37	-1560.1702	3194.3405	0.0000	1.0000
F81	38	-1573.4961	3222.9923	0.0000	1.0000
SYM	40	-1575.0858	3230.1717	0.0000	1.0000
K2P	36	-1579.4043	3230.8086	0.0000	1.0000
JC69	35	-1594.9515	3259.9031	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1515.0949	3266.8379	0.8481	0.8481
HKYG	40	-1517.1588	3270.9656	0.1077	0.9558
HKYIG	41	-1515.0915	3272.7474	0.0442	1.0000
GTRI	44	-1514.7652	3289.8433	0.0000	1.0000
GTRG	44	-1517.0697	3294.4522	0.0000	1.0000
GTRIG	45	-1514.7702	3295.7695	0.0000	1.0000
F81I	39	-1537.2481	3305.2281	0.0000	1.0000
K2PI	37	-1543.9564	3306.8122	0.0000	1.0000
K2PG	37	-1544.0838	3307.0670	0.0000	1.0000
F81G	39	-1538.6258	3307.9836	0.0000	1.0000
F81IG	40	-1537.2283	3311.1046	0.0000	1.0000
K2PIG	38	-1543.6341	3312.0839	0.0000	1.0000
SYMG	41	-1540.2034	3322.9711	0.0000	1.0000
SYMI	41	-1541.2312	3325.0267	0.0000	1.0000
SYMIG	42	-1540.2033	3328.8870	0.0000	1.0000
JC69I	36	-1560.3546	3333.6924	0.0000	1.0000
JC69G	36	-1560.7889	3334.5611	0.0000	1.0000
JC69IG	37	-1560.1702	3339.2400	0.0000	1.0000
HKY	39	-1554.5032	3339.7383	0.0000	1.0000
GTR	43	-1553.9164	3362.2294	0.0000	1.0000
K2P	36	-1579.4043	3371.7919	0.0000	1.0000
F81	38	-1573.4961	3371.8079	0.0000	1.0000
SYM	40	-1575.0858	3386.8198	0.0000	1.0000
JC69	35	-1594.9515	3396.9701	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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_19948.phylip = 1 - 371;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_19948.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_19948.phylip = 1 - 371;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_19948.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_19948.phylip = 1 - 371;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_19948.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

