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

Minimum AIC  model: GTRG
Minimum AICc model: SYMG
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
SYMG	41	-2177.4353	4441.8047	0.1825	0.1825
K2PIG	38	-2180.8517	4441.9316	0.1713	0.3539
GTRG	44	-2174.1431	4441.9840	0.1669	0.5208
HKYIG	41	-2177.8129	4442.5598	0.1251	0.6459
SYMIG	42	-2177.3625	4443.9072	0.0638	0.7097
SYMI	41	-2178.5363	4444.0067	0.0607	0.7704
GTRIG	45	-2174.0465	4444.0583	0.0592	0.8296
K2PG	37	-2183.1060	4444.2178	0.0546	0.8842
GTRI	44	-2175.4052	4444.5082	0.0472	0.9314
HKYG	40	-2180.1138	4444.9200	0.0384	0.9699
K2PI	37	-2184.2062	4446.4181	0.0182	0.9881
HKYI	40	-2181.2829	4447.2582	0.0119	1.0000
SYM	40	-2199.4196	4483.5317	0.0000	1.0000
GTR	43	-2196.4244	4484.2856	0.0000	1.0000
K2P	36	-2205.7855	4487.3605	0.0000	1.0000
HKY	39	-2202.8574	4488.1720	0.0000	1.0000
JC69G	36	-2229.5812	4534.9518	0.0000	1.0000
JC69I	36	-2230.3821	4536.5536	0.0000	1.0000
F81G	39	-2227.0726	4536.6023	0.0000	1.0000
JC69IG	37	-2229.5590	4537.1236	0.0000	1.0000
F81I	39	-2227.9328	4538.3227	0.0000	1.0000
F81IG	40	-2227.0414	4538.7753	0.0000	1.0000
JC69	35	-2250.2268	4574.0332	0.0000	1.0000
F81	38	-2247.8995	4576.0273	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2174.1431	4436.2861	0.2298	0.2298
SYMG	41	-2177.4353	4436.8706	0.1715	0.4013
HKYIG	41	-2177.8129	4437.6257	0.1176	0.5189
K2PIG	38	-2180.8517	4437.7034	0.1131	0.6320
GTRIG	45	-2174.0465	4438.0929	0.0931	0.7251
SYMIG	42	-2177.3625	4438.7250	0.0679	0.7930
GTRI	44	-2175.4052	4438.8104	0.0650	0.8580
SYMI	41	-2178.5363	4439.0726	0.0570	0.9150
K2PG	37	-2183.1060	4440.2121	0.0323	0.9473
HKYG	40	-2180.1138	4440.2276	0.0320	0.9793
K2PI	37	-2184.2062	4442.4124	0.0107	0.9901
HKYI	40	-2181.2829	4442.5658	0.0099	1.0000
SYM	40	-2199.4196	4478.8393	0.0000	1.0000
GTR	43	-2196.4244	4478.8488	0.0000	1.0000
K2P	36	-2205.7855	4483.5710	0.0000	1.0000
HKY	39	-2202.8574	4483.7149	0.0000	1.0000
JC69G	36	-2229.5812	4531.1623	0.0000	1.0000
F81G	39	-2227.0726	4532.1452	0.0000	1.0000
JC69I	36	-2230.3821	4532.7641	0.0000	1.0000
JC69IG	37	-2229.5590	4533.1179	0.0000	1.0000
F81I	39	-2227.9328	4533.8656	0.0000	1.0000
F81IG	40	-2227.0414	4534.0828	0.0000	1.0000
JC69	35	-2250.2268	4570.4536	0.0000	1.0000
F81	38	-2247.8995	4571.7991	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-2183.1060	4610.6581	0.5933	0.5933
K2PIG	38	-2180.8517	4612.7561	0.2078	0.8011
K2PI	37	-2184.2062	4612.8584	0.1975	0.9986
HKYG	40	-2180.1138	4624.4936	0.0006	0.9992
SYMG	41	-2177.4353	4625.7433	0.0003	0.9995
HKYIG	41	-2177.8129	4626.4984	0.0002	0.9997
HKYI	40	-2181.2829	4626.8318	0.0002	0.9999
SYMI	41	-2178.5363	4627.9453	0.0001	1.0000
SYMIG	42	-2177.3625	4632.2043	0.0000	1.0000
GTRG	44	-2174.1431	4638.9787	0.0000	1.0000
GTRI	44	-2175.4052	4641.5030	0.0000	1.0000
GTRIG	45	-2174.0465	4645.3922	0.0000	1.0000
K2P	36	-2205.7855	4649.4104	0.0000	1.0000
SYM	40	-2199.4196	4663.1053	0.0000	1.0000
HKY	39	-2202.8574	4663.3743	0.0000	1.0000
GTR	43	-2196.4244	4676.9347	0.0000	1.0000
JC69G	36	-2229.5812	4697.0017	0.0000	1.0000
JC69I	36	-2230.3821	4698.6035	0.0000	1.0000
JC69IG	37	-2229.5590	4703.5640	0.0000	1.0000
F81G	39	-2227.0726	4711.8046	0.0000	1.0000
F81I	39	-2227.9328	4713.5250	0.0000	1.0000
F81IG	40	-2227.0414	4718.3488	0.0000	1.0000
JC69	35	-2250.2268	4731.6864	0.0000	1.0000
F81	38	-2247.8995	4746.8518	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7218.phylip = 1 - 740;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7218.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7218.phylip = 1 - 740;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7218.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal);
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/chr2_7218.phylip = 1 - 740;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7218.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

