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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: SYMG


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
GTRIG	45	-4077.1672	8248.6157	0.7171	0.7171
GTRG	44	-4079.5473	8251.1856	0.1984	0.9155
SYMIG	42	-4082.8495	8253.4228	0.0648	0.9803
SYMG	41	-4085.1295	8255.8059	0.0197	1.0000
GTRI	44	-4094.1183	8280.3275	0.0000	1.0000
SYMI	41	-4098.3089	8282.1646	0.0000	1.0000
K2PIG	38	-4105.0109	8289.0649	0.0000	1.0000
K2PG	37	-4107.2940	8291.4720	0.0000	1.0000
HKYIG	41	-4104.0435	8293.6338	0.0000	1.0000
HKYG	40	-4106.3413	8296.0571	0.0000	1.0000
K2PI	37	-4119.9514	8316.7869	0.0000	1.0000
HKYI	40	-4119.1702	8321.7148	0.0000	1.0000
JC69IG	37	-4218.6280	8514.1402	0.0000	1.0000
JC69G	36	-4220.5855	8515.9004	0.0000	1.0000
F81IG	40	-4218.0134	8519.4013	0.0000	1.0000
F81G	39	-4219.9725	8521.1515	0.0000	1.0000
SYM	40	-4220.6545	8524.6836	0.0000	1.0000
GTR	43	-4218.2844	8526.4738	0.0000	1.0000
JC69I	36	-4232.7545	8540.2384	0.0000	1.0000
F81I	39	-4232.1873	8545.5812	0.0000	1.0000
K2P	36	-4239.6029	8553.9353	0.0000	1.0000
HKY	39	-4238.8632	8558.9330	0.0000	1.0000
JC69	35	-4348.8229	8770.2252	0.0000	1.0000
F81	38	-4348.1490	8775.3412	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-4077.1672	8244.3344	0.7468	0.7468
GTRG	44	-4079.5473	8247.0947	0.1879	0.9347
SYMIG	42	-4082.8495	8249.6991	0.0511	0.9858
SYMG	41	-4085.1295	8252.2591	0.0142	1.0000
GTRI	44	-4094.1183	8276.2366	0.0000	1.0000
SYMI	41	-4098.3089	8278.6177	0.0000	1.0000
K2PIG	38	-4105.0109	8286.0217	0.0000	1.0000
K2PG	37	-4107.2940	8288.5879	0.0000	1.0000
HKYIG	41	-4104.0435	8290.0869	0.0000	1.0000
HKYG	40	-4106.3413	8292.6826	0.0000	1.0000
K2PI	37	-4119.9514	8313.9028	0.0000	1.0000
HKYI	40	-4119.1702	8318.3403	0.0000	1.0000
JC69IG	37	-4218.6280	8511.2561	0.0000	1.0000
JC69G	36	-4220.5855	8513.1709	0.0000	1.0000
F81IG	40	-4218.0134	8516.0268	0.0000	1.0000
F81G	39	-4219.9725	8517.9449	0.0000	1.0000
SYM	40	-4220.6545	8521.3091	0.0000	1.0000
GTR	43	-4218.2844	8522.5688	0.0000	1.0000
JC69I	36	-4232.7545	8537.5089	0.0000	1.0000
F81I	39	-4232.1873	8542.3746	0.0000	1.0000
K2P	36	-4239.6029	8551.2058	0.0000	1.0000
HKY	39	-4238.8632	8555.7264	0.0000	1.0000
JC69	35	-4348.8229	8767.6459	0.0000	1.0000
F81	38	-4348.1490	8772.2981	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	41	-4085.1295	8454.0066	0.7584	0.7584
SYMIG	42	-4082.8495	8456.3673	0.2330	0.9914
GTRG	44	-4079.5473	8463.6042	0.0062	0.9976
GTRIG	45	-4077.1672	8465.7646	0.0021	0.9998
K2PG	37	-4107.2940	8470.6528	0.0002	0.9999
K2PIG	38	-4105.0109	8473.0073	0.0001	1.0000
SYMI	41	-4098.3089	8480.3652	0.0000	1.0000
HKYG	40	-4106.3413	8489.5095	0.0000	1.0000
HKYIG	41	-4104.0435	8491.8344	0.0000	1.0000
GTRI	44	-4094.1183	8492.7461	0.0000	1.0000
K2PI	37	-4119.9514	8495.9677	0.0000	1.0000
HKYI	40	-4119.1702	8515.1672	0.0000	1.0000
JC69G	36	-4220.5855	8690.3151	0.0000	1.0000
JC69IG	37	-4218.6280	8693.3209	0.0000	1.0000
F81G	39	-4219.9725	8709.8511	0.0000	1.0000
F81IG	40	-4218.0134	8712.8537	0.0000	1.0000
JC69I	36	-4232.7545	8714.6531	0.0000	1.0000
SYM	40	-4220.6545	8718.1359	0.0000	1.0000
K2P	36	-4239.6029	8728.3499	0.0000	1.0000
GTR	43	-4218.2844	8734.1576	0.0000	1.0000
F81I	39	-4232.1873	8734.2808	0.0000	1.0000
HKY	39	-4238.8632	8747.6326	0.0000	1.0000
JC69	35	-4348.8229	8939.8694	0.0000	1.0000
F81	38	-4348.1490	8959.2836	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/chr2_17644.phylip = 1 - 1013;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_17644.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 (GTRIG)]

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

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


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

