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

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	37	-1315.6174	2712.5386	0.3624	0.3624
K2PI	37	-1315.9584	2713.2208	0.2577	0.6201
HKYG	40	-1313.0093	2714.6050	0.1290	0.7490
K2PIG	38	-1315.6236	2714.9659	0.1077	0.8567
HKYI	40	-1313.3588	2715.3039	0.0909	0.9476
HKYIG	41	-1313.0142	2717.0679	0.0376	0.9853
SYMG	41	-1314.9899	2721.0192	0.0052	0.9905
SYMI	41	-1315.3508	2721.7411	0.0036	0.9941
GTRG	44	-1312.3869	2723.2499	0.0017	0.9958
SYMIG	42	-1314.9951	2723.4956	0.0015	0.9974
GTRI	44	-1312.7674	2724.0110	0.0012	0.9985
K2P	36	-1323.0128	2724.9271	0.0007	0.9993
GTRIG	45	-1312.3904	2725.7623	0.0005	0.9998
HKY	39	-1320.5962	2727.3387	0.0002	1.0000
SYM	40	-1322.1009	2732.7882	0.0000	1.0000
GTR	43	-1319.5686	2735.1214	0.0000	1.0000
JC69G	36	-1336.9559	2752.8134	0.0000	1.0000
JC69I	36	-1337.3911	2753.6838	0.0000	1.0000
JC69IG	37	-1336.9647	2755.2333	0.0000	1.0000
F81G	39	-1334.9174	2755.9810	0.0000	1.0000
F81I	39	-1335.3601	2756.8665	0.0000	1.0000
F81IG	40	-1334.9229	2758.4322	0.0000	1.0000
JC69	35	-1343.9776	2764.4667	0.0000	1.0000
F81	38	-1341.9306	2767.5799	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-1315.6174	2705.2347	0.2780	0.2780
K2PI	37	-1315.9584	2705.9169	0.1977	0.4757
HKYG	40	-1313.0093	2706.0186	0.1879	0.6636
HKYI	40	-1313.3588	2706.7175	0.1325	0.7960
K2PIG	38	-1315.6236	2707.2471	0.1016	0.8977
HKYIG	41	-1313.0142	2708.0285	0.0688	0.9664
SYMG	41	-1314.9899	2711.9799	0.0095	0.9760
SYMI	41	-1315.3508	2712.7017	0.0066	0.9826
GTRG	44	-1312.3869	2712.7737	0.0064	0.9890
GTRI	44	-1312.7674	2713.5348	0.0044	0.9934
SYMIG	42	-1314.9951	2713.9903	0.0035	0.9969
GTRIG	45	-1312.3904	2714.7809	0.0024	0.9992
K2P	36	-1323.0128	2718.0255	0.0005	0.9997
HKY	39	-1320.5962	2719.1925	0.0003	1.0000
SYM	40	-1322.1009	2724.2019	0.0000	1.0000
GTR	43	-1319.5686	2725.1373	0.0000	1.0000
JC69G	36	-1336.9559	2745.9119	0.0000	1.0000
JC69I	36	-1337.3911	2746.7823	0.0000	1.0000
F81G	39	-1334.9174	2747.8348	0.0000	1.0000
JC69IG	37	-1336.9647	2747.9294	0.0000	1.0000
F81I	39	-1335.3601	2748.7203	0.0000	1.0000
F81IG	40	-1334.9229	2749.8458	0.0000	1.0000
JC69	35	-1343.9776	2757.9551	0.0000	1.0000
F81	38	-1341.9306	2759.8611	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1315.6174	2854.9875	0.5635	0.5635
K2PI	37	-1315.9584	2855.6697	0.4006	0.9641
K2PIG	38	-1315.6236	2861.0472	0.0272	0.9913
K2P	36	-1323.0128	2863.7309	0.0071	0.9984
HKYG	40	-1313.0093	2867.9135	0.0009	0.9993
HKYI	40	-1313.3588	2868.6124	0.0006	0.9999
HKYIG	41	-1313.0142	2873.9707	0.0000	1.0000
HKY	39	-1320.5962	2877.0400	0.0000	1.0000
SYMG	41	-1314.9899	2877.9221	0.0000	1.0000
SYMI	41	-1315.3508	2878.6439	0.0000	1.0000
SYMIG	42	-1314.9951	2883.9799	0.0000	1.0000
SYM	40	-1322.1009	2886.0967	0.0000	1.0000
GTRG	44	-1312.3869	2890.8581	0.0000	1.0000
JC69G	36	-1336.9559	2891.6173	0.0000	1.0000
GTRI	44	-1312.7674	2891.6191	0.0000	1.0000
JC69I	36	-1337.3911	2892.4877	0.0000	1.0000
GTRIG	45	-1312.3904	2896.9126	0.0000	1.0000
JC69IG	37	-1336.9647	2897.6822	0.0000	1.0000
GTR	43	-1319.5686	2899.1743	0.0000	1.0000
JC69	35	-1343.9776	2899.6131	0.0000	1.0000
F81G	39	-1334.9174	2905.6823	0.0000	1.0000
F81I	39	-1335.3601	2906.5678	0.0000	1.0000
F81IG	40	-1334.9229	2911.7407	0.0000	1.0000
F81	38	-1341.9306	2913.6613	0.0000	1.0000



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

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

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


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

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


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

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


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

