
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:55 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/chr5_10413.phylip" (19 taxa, 476 characters)

Minimum AIC  model: K2PIG
Minimum AICc model: K2PIG
Minimum BIC  model: K2PG


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
K2PIG	38	-1605.4566	3293.6957	0.2826	0.2826
K2PG	37	-1606.7707	3293.9615	0.2475	0.5301
SYMIG	42	-1601.5082	3295.3582	0.1231	0.6532
SYMG	41	-1602.8203	3295.5761	0.1104	0.7636
K2PI	37	-1608.1667	3296.7535	0.0613	0.8249
HKYIG	41	-1603.6140	3297.1635	0.0499	0.8748
HKYG	40	-1605.0275	3297.5952	0.0402	0.9150
GTRIG	45	-1599.3704	3298.3688	0.0273	0.9423
SYMI	41	-1604.3744	3298.6843	0.0233	0.9657
GTRG	44	-1600.8243	3298.8366	0.0216	0.9873
HKYI	40	-1606.5413	3300.6229	0.0089	0.9961
GTRI	44	-1602.5469	3302.2818	0.0039	1.0000
JC69IG	37	-1657.4146	3395.2492	0.0000	1.0000
JC69G	36	-1658.6167	3395.3017	0.0000	1.0000
JC69I	36	-1659.9740	3398.0162	0.0000	1.0000
F81IG	40	-1655.9649	3399.4700	0.0000	1.0000
F81G	39	-1657.2140	3399.5839	0.0000	1.0000
F81I	39	-1658.5820	3402.3199	0.0000	1.0000
SYM	40	-1662.6102	3412.7606	0.0000	1.0000
K2P	36	-1668.2722	3414.6127	0.0000	1.0000
GTR	43	-1660.6986	3416.1564	0.0000	1.0000
HKY	39	-1666.4408	3418.0376	0.0000	1.0000
JC69	35	-1718.1274	3511.9821	0.0000	1.0000
F81	38	-1716.8321	3516.4469	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PIG	38	-1605.4566	3286.9131	0.2011	0.2011
SYMIG	42	-1601.5082	3287.0164	0.1910	0.3921
K2PG	37	-1606.7707	3287.5414	0.1469	0.5390
SYMG	41	-1602.8203	3287.6406	0.1398	0.6788
GTRIG	45	-1599.3704	3288.7409	0.0806	0.7595
HKYIG	41	-1603.6140	3289.2281	0.0632	0.8227
GTRG	44	-1600.8243	3289.6486	0.0512	0.8739
HKYG	40	-1605.0275	3290.0550	0.0418	0.9157
K2PI	37	-1608.1667	3290.3334	0.0364	0.9521
SYMI	41	-1604.3744	3290.7488	0.0296	0.9817
HKYI	40	-1606.5413	3293.0826	0.0092	0.9909
GTRI	44	-1602.5469	3293.0939	0.0091	1.0000
JC69IG	37	-1657.4146	3388.8291	0.0000	1.0000
JC69G	36	-1658.6167	3389.2333	0.0000	1.0000
F81IG	40	-1655.9649	3391.9297	0.0000	1.0000
JC69I	36	-1659.9740	3391.9479	0.0000	1.0000
F81G	39	-1657.2140	3392.4279	0.0000	1.0000
F81I	39	-1658.5820	3395.1639	0.0000	1.0000
SYM	40	-1662.6102	3405.2203	0.0000	1.0000
GTR	43	-1660.6986	3407.3971	0.0000	1.0000
K2P	36	-1668.2722	3408.5444	0.0000	1.0000
HKY	39	-1666.4408	3410.8816	0.0000	1.0000
JC69	35	-1718.1274	3506.2548	0.0000	1.0000
F81	38	-1716.8321	3509.6643	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1606.7707	3441.6618	0.7046	0.7046
K2PI	37	-1608.1667	3444.4538	0.1744	0.8790
K2PIG	38	-1605.4566	3445.1990	0.1202	0.9992
HKYG	40	-1605.0275	3456.6717	0.0004	0.9996
SYMG	41	-1602.8203	3458.4228	0.0002	0.9998
HKYI	40	-1606.5413	3459.6994	0.0001	0.9999
HKYIG	41	-1603.6140	3460.0102	0.0001	0.9999
SYMI	41	-1604.3744	3461.5309	0.0000	1.0000
SYMIG	42	-1601.5082	3461.9639	0.0000	1.0000
GTRG	44	-1600.8243	3472.9270	0.0000	1.0000
GTRIG	45	-1599.3704	3476.1847	0.0000	1.0000
GTRI	44	-1602.5469	3476.3722	0.0000	1.0000
JC69G	36	-1658.6167	3539.1884	0.0000	1.0000
JC69I	36	-1659.9740	3541.9029	0.0000	1.0000
JC69IG	37	-1657.4146	3542.9496	0.0000	1.0000
F81G	39	-1657.2140	3554.8792	0.0000	1.0000
F81I	39	-1658.5820	3557.6152	0.0000	1.0000
K2P	36	-1668.2722	3558.4994	0.0000	1.0000
F81IG	40	-1655.9649	3558.5465	0.0000	1.0000
SYM	40	-1662.6102	3571.8371	0.0000	1.0000
HKY	39	-1666.4408	3573.3329	0.0000	1.0000
GTR	43	-1660.6986	3586.5101	0.0000	1.0000
JC69	35	-1718.1274	3652.0444	0.0000	1.0000
F81	38	-1716.8321	3667.9502	0.0000	1.0000



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

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

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


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

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


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

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


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

