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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	40	-1205.2258	2498.0972	0.4841	0.4841
HKYI	40	-1205.5945	2498.8347	0.3348	0.8188
HKYIG	41	-1205.1432	2500.3330	0.1583	0.9771
GTRG	44	-1204.1674	2505.6524	0.0111	0.9882
GTRI	44	-1204.5952	2506.5081	0.0072	0.9954
GTRIG	45	-1204.0504	2507.8651	0.0037	0.9991
K2PG	37	-1216.0638	2512.6368	0.0003	0.9994
K2PI	37	-1216.4051	2513.3194	0.0002	0.9996
HKY	39	-1214.0511	2513.3580	0.0002	0.9999
K2PIG	38	-1215.9965	2514.8701	0.0001	1.0000
SYMG	41	-1215.3124	2520.6716	0.0000	1.0000
GTR	43	-1213.1190	2521.1206	0.0000	1.0000
SYMI	41	-1215.6818	2521.4104	0.0000	1.0000
SYMIG	42	-1215.2402	2522.9395	0.0000	1.0000
K2P	36	-1224.4048	2526.9620	0.0000	1.0000
SYM	40	-1223.4840	2534.6137	0.0000	1.0000
F81G	39	-1230.8121	2546.8801	0.0000	1.0000
F81I	39	-1231.1671	2547.5901	0.0000	1.0000
F81IG	40	-1230.7455	2549.1367	0.0000	1.0000
JC69G	36	-1240.2646	2558.6817	0.0000	1.0000
JC69I	36	-1240.6057	2559.3638	0.0000	1.0000
F81	38	-1238.9426	2560.7623	0.0000	1.0000
JC69IG	37	-1240.2033	2560.9159	0.0000	1.0000
JC69	35	-1248.2231	2572.2526	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1205.2258	2490.4515	0.4542	0.4542
HKYI	40	-1205.5945	2491.1890	0.3142	0.7684
HKYIG	41	-1205.1432	2492.2863	0.1815	0.9499
GTRG	44	-1204.1674	2496.3348	0.0240	0.9739
GTRI	44	-1204.5952	2497.1904	0.0156	0.9895
GTRIG	45	-1204.0504	2498.1009	0.0099	0.9994
HKY	39	-1214.0511	2506.1021	0.0002	0.9996
K2PG	37	-1216.0638	2506.1275	0.0002	0.9998
K2PI	37	-1216.4051	2506.8101	0.0001	0.9999
K2PIG	38	-1215.9965	2507.9930	0.0001	1.0000
GTR	43	-1213.1190	2512.2379	0.0000	1.0000
SYMG	41	-1215.3124	2512.6249	0.0000	1.0000
SYMI	41	-1215.6818	2513.3637	0.0000	1.0000
SYMIG	42	-1215.2402	2514.4805	0.0000	1.0000
K2P	36	-1224.4048	2520.8096	0.0000	1.0000
SYM	40	-1223.4840	2526.9680	0.0000	1.0000
F81G	39	-1230.8121	2539.6243	0.0000	1.0000
F81I	39	-1231.1671	2540.3343	0.0000	1.0000
F81IG	40	-1230.7455	2541.4910	0.0000	1.0000
JC69G	36	-1240.2646	2552.5292	0.0000	1.0000
JC69I	36	-1240.6057	2553.2114	0.0000	1.0000
F81	38	-1238.9426	2553.8853	0.0000	1.0000
JC69IG	37	-1240.2033	2554.4066	0.0000	1.0000
JC69	35	-1248.2231	2566.4461	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1205.2258	2656.5608	0.4766	0.4766
HKYI	40	-1205.5945	2657.2983	0.3296	0.8062
K2PG	37	-1216.0638	2659.7786	0.0954	0.9016
K2PI	37	-1216.4051	2660.4612	0.0678	0.9694
HKYIG	41	-1205.1432	2662.5484	0.0239	0.9933
K2PIG	38	-1215.9965	2665.7969	0.0047	0.9980
HKY	39	-1214.0511	2668.0587	0.0015	0.9995
K2P	36	-1224.4048	2670.3080	0.0005	1.0000
GTRG	44	-1204.1674	2679.0550	0.0000	1.0000
GTRI	44	-1204.5952	2679.9107	0.0000	1.0000
SYMG	41	-1215.3124	2682.8869	0.0000	1.0000
SYMI	41	-1215.6818	2683.6257	0.0000	1.0000
GTRIG	45	-1204.0504	2684.9739	0.0000	1.0000
SYMIG	42	-1215.2402	2688.8953	0.0000	1.0000
GTR	43	-1213.1190	2690.8054	0.0000	1.0000
SYM	40	-1223.4840	2693.0773	0.0000	1.0000
F81G	39	-1230.8121	2701.5808	0.0000	1.0000
JC69G	36	-1240.2646	2702.0276	0.0000	1.0000
F81I	39	-1231.1671	2702.2909	0.0000	1.0000
JC69I	36	-1240.6057	2702.7098	0.0000	1.0000
F81IG	40	-1230.7455	2707.6003	0.0000	1.0000
JC69IG	37	-1240.2033	2708.0577	0.0000	1.0000
F81	38	-1238.9426	2711.6891	0.0000	1.0000
JC69	35	-1248.2231	2711.7918	0.0000	1.0000



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

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

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


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

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


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

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


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

