
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:25 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/chr1_37896.phylip" (19 taxa, 416 characters)

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


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.8699	3120.4864	0.5751	0.5751
HKYG	40	-1516.9953	3122.7373	0.1866	0.7618
HKYIG	41	-1515.8585	3122.9256	0.1699	0.9316
K2PI	37	-1522.4461	3126.3313	0.0309	0.9626
GTRI	44	-1514.4506	3127.5751	0.0166	0.9792
K2PG	37	-1523.4678	3128.3747	0.0111	0.9903
GTRG	44	-1515.5340	3129.7418	0.0056	0.9959
K2PIG	38	-1524.0792	3132.0204	0.0018	0.9977
GTRIG	45	-1515.5342	3132.2576	0.0016	0.9993
SYMI	41	-1521.7663	3134.7412	0.0005	0.9998
SYMG	41	-1522.8799	3136.9684	0.0002	1.0000
SYMIG	42	-1522.8801	3139.4439	0.0000	1.0000
F81I	39	-1542.8495	3171.9970	0.0000	1.0000
F81G	39	-1543.9501	3174.1980	0.0000	1.0000
JC69I	36	-1548.1767	3175.3823	0.0000	1.0000
F81IG	40	-1543.9503	3176.6472	0.0000	1.0000
JC69G	36	-1549.2290	3177.4870	0.0000	1.0000
JC69IG	37	-1549.2292	3179.8975	0.0000	1.0000
HKY	39	-1567.9502	3222.1982	0.0000	1.0000
K2P	36	-1573.6635	3226.3560	0.0000	1.0000
GTR	43	-1566.9170	3230.0060	0.0000	1.0000
SYM	40	-1573.0893	3234.9253	0.0000	1.0000
F81	38	-1592.4733	3268.8086	0.0000	1.0000
JC69	35	-1597.3324	3271.2964	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1515.8699	3111.7397	0.5409	0.5409
HKYIG	41	-1515.8585	3113.7170	0.2012	0.7421
HKYG	40	-1516.9953	3113.9907	0.1755	0.9176
GTRI	44	-1514.4506	3116.9013	0.0410	0.9586
K2PI	37	-1522.4461	3118.8921	0.0151	0.9737
GTRG	44	-1515.5340	3119.0680	0.0139	0.9876
K2PG	37	-1523.4678	3120.9356	0.0054	0.9930
GTRIG	45	-1515.5342	3121.0684	0.0051	0.9981
K2PIG	38	-1524.0792	3124.1583	0.0011	0.9992
SYMI	41	-1521.7663	3125.5326	0.0005	0.9998
SYMG	41	-1522.8799	3127.7598	0.0002	0.9999
SYMIG	42	-1522.8801	3129.7603	0.0001	1.0000
F81I	39	-1542.8495	3163.6991	0.0000	1.0000
F81G	39	-1543.9501	3165.9001	0.0000	1.0000
F81IG	40	-1543.9503	3167.9005	0.0000	1.0000
JC69I	36	-1548.1767	3168.3533	0.0000	1.0000
JC69G	36	-1549.2290	3170.4580	0.0000	1.0000
JC69IG	37	-1549.2292	3172.4583	0.0000	1.0000
HKY	39	-1567.9502	3213.9003	0.0000	1.0000
K2P	36	-1573.6635	3219.3270	0.0000	1.0000
GTR	43	-1566.9170	3219.8340	0.0000	1.0000
SYM	40	-1573.0893	3226.1786	0.0000	1.0000
F81	38	-1592.4733	3260.9465	0.0000	1.0000
JC69	35	-1597.3324	3264.6648	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-1522.4461	3268.0275	0.6730	0.6730
K2PG	37	-1523.4678	3270.0709	0.2423	0.9153
HKYI	40	-1515.8699	3272.9671	0.0569	0.9722
HKYG	40	-1516.9953	3275.2181	0.0185	0.9907
K2PIG	38	-1524.0792	3277.3244	0.0064	0.9972
HKYIG	41	-1515.8585	3278.9751	0.0028	1.0000
SYMI	41	-1521.7663	3290.7907	0.0000	1.0000
SYMG	41	-1522.8799	3293.0179	0.0000	1.0000
GTRI	44	-1514.4506	3294.2514	0.0000	1.0000
GTRG	44	-1515.5340	3296.4181	0.0000	1.0000
SYMIG	42	-1522.8801	3299.0491	0.0000	1.0000
GTRIG	45	-1515.5342	3302.4493	0.0000	1.0000
JC69I	36	-1548.1767	3313.4580	0.0000	1.0000
JC69G	36	-1549.2290	3315.5626	0.0000	1.0000
F81I	39	-1542.8495	3320.8958	0.0000	1.0000
JC69IG	37	-1549.2292	3321.5937	0.0000	1.0000
F81G	39	-1543.9501	3323.0968	0.0000	1.0000
F81IG	40	-1543.9503	3329.1279	0.0000	1.0000
K2P	36	-1573.6635	3364.4317	0.0000	1.0000
HKY	39	-1567.9502	3371.0970	0.0000	1.0000
SYM	40	-1573.0893	3387.4060	0.0000	1.0000
GTR	43	-1566.9170	3393.1535	0.0000	1.0000
JC69	35	-1597.3324	3405.7388	0.0000	1.0000
F81	38	-1592.4733	3414.1126	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/chr1_37896.phylip = 1 - 416;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_37896.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/chr1_37896.phylip = 1 - 416;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_37896.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 (K2PI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_37896.phylip = 1 - 416;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_37896.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

