
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:59 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_21314.phylip" (19 taxa, 603 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	-3161.1730	6417.4428	0.9616	0.9616
GTRIG	45	-3163.2288	6423.8904	0.0383	0.9999
HKYG	40	-3175.2718	6436.3800	0.0001	1.0000
HKYIG	41	-3177.6370	6443.4130	0.0000	1.0000
K2PG	37	-3190.3739	6459.7248	0.0000	1.0000
K2PIG	38	-3190.0100	6461.2753	0.0000	1.0000
SYMG	41	-3186.6094	6461.3579	0.0000	1.0000
SYMIG	42	-3186.2004	6462.8508	0.0000	1.0000
GTRI	44	-3184.6885	6464.4738	0.0000	1.0000
HKYI	40	-3192.1632	6470.1627	0.0000	1.0000
K2PI	37	-3206.0532	6491.0834	0.0000	1.0000
SYMI	41	-3202.5981	6493.3353	0.0000	1.0000
F81G	39	-3245.5336	6574.6089	0.0000	1.0000
F81IG	40	-3245.1086	6576.0536	0.0000	1.0000
JC69G	36	-3257.4194	6591.5456	0.0000	1.0000
JC69IG	37	-3257.0484	6593.0738	0.0000	1.0000
GTR	43	-3250.7609	6594.2910	0.0000	1.0000
HKY	39	-3261.1105	6605.7627	0.0000	1.0000
F81I	39	-3261.7193	6606.9803	0.0000	1.0000
JC69I	36	-3272.6687	6622.0442	0.0000	1.0000
K2P	36	-3278.3158	6633.3384	0.0000	1.0000
SYM	40	-3274.5337	6634.9037	0.0000	1.0000
F81	38	-3329.4292	6740.1138	0.0000	1.0000
JC69	35	-3341.7673	6757.9791	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3161.1730	6410.3460	0.9550	0.9550
GTRIG	45	-3163.2288	6416.4577	0.0450	1.0000
HKYG	40	-3175.2718	6430.5437	0.0000	1.0000
HKYIG	41	-3177.6370	6437.2740	0.0000	1.0000
K2PG	37	-3190.3739	6454.7478	0.0000	1.0000
SYMG	41	-3186.6094	6455.2189	0.0000	1.0000
K2PIG	38	-3190.0100	6456.0200	0.0000	1.0000
SYMIG	42	-3186.2004	6456.4008	0.0000	1.0000
GTRI	44	-3184.6885	6457.3770	0.0000	1.0000
HKYI	40	-3192.1632	6464.3264	0.0000	1.0000
K2PI	37	-3206.0532	6486.1064	0.0000	1.0000
SYMI	41	-3202.5981	6487.1963	0.0000	1.0000
F81G	39	-3245.5336	6569.0672	0.0000	1.0000
F81IG	40	-3245.1086	6570.2173	0.0000	1.0000
JC69G	36	-3257.4194	6586.8388	0.0000	1.0000
GTR	43	-3250.7609	6587.5217	0.0000	1.0000
JC69IG	37	-3257.0484	6588.0968	0.0000	1.0000
HKY	39	-3261.1105	6600.2209	0.0000	1.0000
F81I	39	-3261.7193	6601.4385	0.0000	1.0000
JC69I	36	-3272.6687	6617.3375	0.0000	1.0000
K2P	36	-3278.3158	6628.6317	0.0000	1.0000
SYM	40	-3274.5337	6629.0674	0.0000	1.0000
F81	38	-3329.4292	6734.8585	0.0000	1.0000
JC69	35	-3341.7673	6753.5347	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3161.1730	6604.0303	0.7804	0.7804
HKYG	40	-3175.2718	6606.6204	0.2138	0.9942
GTRIG	45	-3163.2288	6614.5440	0.0041	0.9983
K2PG	37	-3190.3739	6617.6188	0.0009	0.9991
HKYIG	41	-3177.6370	6617.7526	0.0008	0.9999
K2PIG	38	-3190.0100	6623.2929	0.0001	1.0000
SYMG	41	-3186.6094	6635.6975	0.0000	1.0000
HKYI	40	-3192.1632	6640.4031	0.0000	1.0000
SYMIG	42	-3186.2004	6641.2813	0.0000	1.0000
K2PI	37	-3206.0532	6648.9773	0.0000	1.0000
GTRI	44	-3184.6885	6651.0614	0.0000	1.0000
SYMI	41	-3202.5981	6667.6749	0.0000	1.0000
F81G	39	-3245.5336	6740.7420	0.0000	1.0000
JC69G	36	-3257.4194	6745.3079	0.0000	1.0000
F81IG	40	-3245.1086	6746.2940	0.0000	1.0000
JC69IG	37	-3257.0484	6750.9677	0.0000	1.0000
HKY	39	-3261.1105	6771.8957	0.0000	1.0000
F81I	39	-3261.7193	6773.1133	0.0000	1.0000
JC69I	36	-3272.6687	6775.8065	0.0000	1.0000
GTR	43	-3250.7609	6776.8042	0.0000	1.0000
K2P	36	-3278.3158	6787.1007	0.0000	1.0000
SYM	40	-3274.5337	6805.1441	0.0000	1.0000
F81	38	-3329.4292	6902.1313	0.0000	1.0000
JC69	35	-3341.7673	6907.6018	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_21314.phylip = 1 - 603;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21314.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/chr2_21314.phylip = 1 - 603;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21314.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/chr2_21314.phylip = 1 - 603;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21314.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

