
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:31 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/chrZ_1402.phylip" (19 taxa, 216 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: K2PIG
Minimum BIC  model: K2PIG


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	-630.6516	1354.0490	0.3305	0.3305
HKYIG	41	-626.7269	1355.2469	0.1816	0.5121
GTRIG	45	-620.9052	1356.1634	0.1148	0.6270
GTRG	44	-622.7578	1356.6736	0.0890	0.7160
K2PG	37	-633.5704	1356.9386	0.0779	0.7939
HKYG	40	-629.5064	1357.7557	0.0518	0.8457
SYMIG	42	-626.7256	1358.3299	0.0389	0.8846
K2PI	37	-634.3650	1358.5278	0.0352	0.9198
SYMG	41	-628.9606	1359.7143	0.0195	0.9392
GTRI	44	-624.3620	1359.8819	0.0179	0.9571
HKYI	40	-630.5854	1359.9137	0.0176	0.9747
JC69IG	37	-635.9583	1361.7144	0.0072	0.9819
SYMI	41	-630.0102	1361.8135	0.0068	0.9887
F81IG	40	-632.0402	1362.8232	0.0041	0.9928
JC69G	36	-638.2031	1363.2890	0.0033	0.9961
JC69I	36	-638.8530	1364.5887	0.0017	0.9978
F81G	39	-634.5441	1364.8155	0.0015	0.9993
F81I	39	-635.3041	1366.3355	0.0007	1.0000
K2P	36	-653.8600	1394.6027	0.0000	1.0000
GTR	43	-643.3266	1394.6532	0.0000	1.0000
HKY	39	-649.5213	1394.7698	0.0000	1.0000
SYM	40	-650.3902	1399.5232	0.0000	1.0000
JC69	35	-657.8311	1399.6623	0.0000	1.0000
F81	38	-653.7533	1400.2524	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-620.9052	1331.8105	0.5333	0.5333
GTRG	44	-622.7578	1333.5157	0.2273	0.7606
HKYIG	41	-626.7269	1335.4538	0.0863	0.8468
GTRI	44	-624.3620	1336.7240	0.0457	0.8926
K2PIG	38	-630.6516	1337.3032	0.0342	0.9268
SYMIG	42	-626.7256	1337.4513	0.0318	0.9585
HKYG	40	-629.5064	1339.0128	0.0146	0.9731
SYMG	41	-628.9606	1339.9212	0.0092	0.9823
K2PG	37	-633.5704	1341.1409	0.0050	0.9874
HKYI	40	-630.5854	1341.1708	0.0049	0.9923
SYMI	41	-630.0102	1342.0204	0.0032	0.9955
K2PI	37	-634.3650	1342.7300	0.0023	0.9978
F81IG	40	-632.0402	1344.0803	0.0012	0.9990
JC69IG	37	-635.9583	1345.9167	0.0005	0.9994
F81G	39	-634.5441	1347.0882	0.0003	0.9997
JC69G	36	-638.2031	1348.4063	0.0001	0.9998
F81I	39	-635.3041	1348.6082	0.0001	0.9999
JC69I	36	-638.8530	1349.7060	0.0001	1.0000
GTR	43	-643.3266	1372.6532	0.0000	1.0000
HKY	39	-649.5213	1377.0425	0.0000	1.0000
K2P	36	-653.8600	1379.7200	0.0000	1.0000
SYM	40	-650.3902	1380.7803	0.0000	1.0000
F81	38	-653.7533	1383.5067	0.0000	1.0000
JC69	35	-657.8311	1385.6623	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	38	-630.6516	1465.5638	0.4098	0.4098
K2PG	37	-633.5704	1466.0262	0.3252	0.7350
K2PI	37	-634.3650	1467.6153	0.1469	0.8820
JC69G	36	-638.2031	1469.9163	0.0465	0.9285
JC69IG	37	-635.9583	1470.8020	0.0299	0.9583
JC69I	36	-638.8530	1471.2160	0.0243	0.9826
HKYIG	41	-626.7269	1473.8402	0.0065	0.9891
HKYG	40	-629.5064	1474.0240	0.0060	0.9951
HKYI	40	-630.5854	1476.1820	0.0020	0.9971
SYMG	41	-628.9606	1478.3077	0.0007	0.9978
F81G	39	-634.5441	1478.7240	0.0006	0.9984
F81IG	40	-632.0402	1479.0915	0.0005	0.9989
SYMIG	42	-626.7256	1479.2130	0.0004	0.9993
F81I	39	-635.3041	1480.2440	0.0003	0.9996
SYMI	41	-630.0102	1480.4068	0.0002	0.9998
GTRG	44	-622.7578	1482.0279	0.0001	0.9999
GTRIG	45	-620.9052	1483.6980	0.0000	1.0000
GTRI	44	-624.3620	1485.2362	0.0000	1.0000
K2P	36	-653.8600	1501.2300	0.0000	1.0000
JC69	35	-657.8311	1503.7970	0.0000	1.0000
HKY	39	-649.5213	1508.6784	0.0000	1.0000
F81	38	-653.7533	1511.7673	0.0000	1.0000
SYM	40	-650.3902	1515.7915	0.0000	1.0000
GTR	43	-643.3266	1517.7902	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_1402.phylip = 1 - 216;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_1402.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) 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/chrZ_1402.phylip = 1 - 216;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_1402.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 (K2PIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_1402.phylip = 1 - 216;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_1402.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;


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

