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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	44	-3634.7733	7362.8549	0.6946	0.6946
GTRIG	45	-3634.4712	7364.4995	0.3052	0.9998
GTRI	44	-3643.1591	7379.6265	0.0002	1.0000
HKYG	40	-3650.0134	7384.4002	0.0000	1.0000
HKYIG	41	-3649.7652	7386.1286	0.0000	1.0000
HKYI	40	-3655.1424	7394.6581	0.0000	1.0000
SYMG	41	-3680.1153	7446.8288	0.0000	1.0000
K2PG	37	-3684.7762	7447.2868	0.0000	1.0000
SYMIG	42	-3679.8986	7448.6261	0.0000	1.0000
K2PIG	38	-3684.5679	7449.0774	0.0000	1.0000
SYMI	41	-3686.6463	7459.8908	0.0000	1.0000
K2PI	37	-3691.1723	7460.0790	0.0000	1.0000
GTR	43	-3709.9670	7510.9996	0.0000	1.0000
HKY	39	-3720.0981	7522.3507	0.0000	1.0000
F81IG	40	-3737.1028	7558.5790	0.0000	1.0000
F81G	39	-3741.8826	7565.9196	0.0000	1.0000
F81I	39	-3747.5721	7577.2986	0.0000	1.0000
SYM	40	-3746.5904	7577.5541	0.0000	1.0000
K2P	36	-3751.5554	7578.6439	0.0000	1.0000
JC69G	36	-3761.4896	7598.5124	0.0000	1.0000
JC69IG	37	-3761.1190	7599.9723	0.0000	1.0000
JC69I	36	-3765.7494	7607.0320	0.0000	1.0000
F81	38	-3804.2431	7688.4276	0.0000	1.0000
JC69	35	-3821.0816	7715.5010	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3634.7733	7357.5466	0.6676	0.6676
GTRIG	45	-3634.4712	7358.9424	0.3322	0.9998
GTRI	44	-3643.1591	7374.3182	0.0002	1.0000
HKYG	40	-3650.0134	7380.0268	0.0000	1.0000
HKYIG	41	-3649.7652	7381.5304	0.0000	1.0000
HKYI	40	-3655.1424	7390.2848	0.0000	1.0000
SYMG	41	-3680.1153	7442.2307	0.0000	1.0000
K2PG	37	-3684.7762	7443.5524	0.0000	1.0000
SYMIG	42	-3679.8986	7443.7972	0.0000	1.0000
K2PIG	38	-3684.5679	7445.1359	0.0000	1.0000
SYMI	41	-3686.6463	7455.2927	0.0000	1.0000
K2PI	37	-3691.1723	7456.3446	0.0000	1.0000
GTR	43	-3709.9670	7505.9340	0.0000	1.0000
HKY	39	-3720.0981	7518.1962	0.0000	1.0000
F81IG	40	-3737.1028	7554.2057	0.0000	1.0000
F81G	39	-3741.8826	7561.7652	0.0000	1.0000
F81I	39	-3747.5721	7573.1441	0.0000	1.0000
SYM	40	-3746.5904	7573.1808	0.0000	1.0000
K2P	36	-3751.5554	7575.1107	0.0000	1.0000
JC69G	36	-3761.4896	7594.9793	0.0000	1.0000
JC69IG	37	-3761.1190	7596.2379	0.0000	1.0000
JC69I	36	-3765.7494	7603.4988	0.0000	1.0000
F81	38	-3804.2431	7684.4861	0.0000	1.0000
JC69	35	-3821.0816	7712.1632	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3634.7733	7563.1717	0.8288	0.8288
HKYG	40	-3650.0134	7566.9587	0.1248	0.9535
GTRIG	45	-3634.4712	7569.2408	0.0399	0.9934
HKYIG	41	-3649.7652	7573.1356	0.0057	0.9991
HKYI	40	-3655.1424	7577.2167	0.0007	0.9998
GTRI	44	-3643.1591	7579.9433	0.0002	1.0000
K2PG	37	-3684.7762	7616.4644	0.0000	1.0000
K2PIG	38	-3684.5679	7622.7212	0.0000	1.0000
K2PI	37	-3691.1723	7629.2566	0.0000	1.0000
SYMG	41	-3680.1153	7633.8359	0.0000	1.0000
SYMIG	42	-3679.8986	7640.0757	0.0000	1.0000
SYMI	41	-3686.6463	7646.8979	0.0000	1.0000
HKY	39	-3720.0981	7700.4549	0.0000	1.0000
GTR	43	-3709.9670	7706.8858	0.0000	1.0000
F81IG	40	-3737.1028	7741.1376	0.0000	1.0000
K2P	36	-3751.5554	7743.3495	0.0000	1.0000
F81G	39	-3741.8826	7744.0238	0.0000	1.0000
F81I	39	-3747.5721	7755.4028	0.0000	1.0000
SYM	40	-3746.5904	7760.1127	0.0000	1.0000
JC69G	36	-3761.4896	7763.2180	0.0000	1.0000
JC69IG	37	-3761.1190	7769.1499	0.0000	1.0000
JC69I	36	-3765.7494	7771.7375	0.0000	1.0000
F81	38	-3804.2431	7862.0714	0.0000	1.0000
JC69	35	-3821.0816	7875.7286	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/chrZ_9608.phylip = 1 - 791;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9608.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 (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9608.phylip = 1 - 791;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9608.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 BIC model (GTRG)]

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


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

