
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:08 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/chr13_3036.phylip" (19 taxa, 347 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	-1044.0564	2178.8318	0.5820	0.5820
HKYIG	41	-1043.7781	2180.8480	0.2124	0.7944
HKYI	40	-1045.7788	2182.2765	0.1040	0.8984
GTRG	44	-1041.0302	2183.1729	0.0664	0.9648
GTRIG	45	-1040.7976	2185.3494	0.0224	0.9872
GTRI	44	-1042.6736	2186.4598	0.0128	1.0000
HKY	39	-1058.9852	2206.1333	0.0000	1.0000
GTR	43	-1055.9641	2210.4167	0.0000	1.0000
K2PG	37	-1068.7871	2220.6746	0.0000	1.0000
K2PI	37	-1070.5807	2224.2617	0.0000	1.0000
K2PIG	38	-1069.3862	2224.3958	0.0000	1.0000
SYMG	41	-1068.1290	2229.5498	0.0000	1.0000
SYMI	41	-1069.8646	2233.0210	0.0000	1.0000
SYMIG	42	-1068.6854	2233.2524	0.0000	1.0000
K2P	36	-1085.2221	2251.0377	0.0000	1.0000
SYM	40	-1086.4598	2263.6385	0.0000	1.0000
F81G	39	-1088.7360	2265.6348	0.0000	1.0000
F81IG	40	-1089.1358	2268.9906	0.0000	1.0000
F81I	39	-1090.4394	2269.0417	0.0000	1.0000
F81	38	-1103.8119	2293.2472	0.0000	1.0000
JC69G	36	-1107.6549	2295.9034	0.0000	1.0000
JC69I	36	-1109.7721	2300.1378	0.0000	1.0000
JC69IG	37	-1108.5832	2300.2667	0.0000	1.0000
JC69	35	-1126.1035	2330.3099	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1044.0564	2168.1128	0.4366	0.4366
HKYIG	41	-1043.7781	2169.5562	0.2121	0.6487
GTRG	44	-1041.0302	2170.0604	0.1649	0.8136
HKYI	40	-1045.7788	2171.5575	0.0780	0.8916
GTRIG	45	-1040.7976	2171.5952	0.0765	0.9681
GTRI	44	-1042.6736	2173.3472	0.0319	1.0000
HKY	39	-1058.9852	2195.9704	0.0000	1.0000
GTR	43	-1055.9641	2197.9283	0.0000	1.0000
K2PG	37	-1068.7871	2211.5743	0.0000	1.0000
K2PIG	38	-1069.3862	2214.7725	0.0000	1.0000
K2PI	37	-1070.5807	2215.1614	0.0000	1.0000
SYMG	41	-1068.1290	2218.2580	0.0000	1.0000
SYMIG	42	-1068.6854	2221.3708	0.0000	1.0000
SYMI	41	-1069.8646	2221.7292	0.0000	1.0000
K2P	36	-1085.2221	2242.4442	0.0000	1.0000
SYM	40	-1086.4598	2252.9195	0.0000	1.0000
F81G	39	-1088.7360	2255.4719	0.0000	1.0000
F81IG	40	-1089.1358	2258.2717	0.0000	1.0000
F81I	39	-1090.4394	2258.8788	0.0000	1.0000
F81	38	-1103.8119	2283.6238	0.0000	1.0000
JC69G	36	-1107.6549	2287.3098	0.0000	1.0000
JC69IG	37	-1108.5832	2291.1664	0.0000	1.0000
JC69I	36	-1109.7721	2291.5443	0.0000	1.0000
JC69	35	-1126.1035	2322.2070	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1044.0564	2322.0858	0.8001	0.8001
HKYI	40	-1045.7788	2325.5305	0.1429	0.9431
HKYIG	41	-1043.7781	2327.3785	0.0567	0.9998
GTRG	44	-1041.0302	2339.4307	0.0001	1.0000
GTRI	44	-1042.6736	2342.7175	0.0000	1.0000
GTRIG	45	-1040.7976	2344.8148	0.0000	1.0000
HKY	39	-1058.9852	2346.0941	0.0000	1.0000
K2PG	37	-1068.7871	2353.9993	0.0000	1.0000
K2PI	37	-1070.5807	2357.5864	0.0000	1.0000
K2PIG	38	-1069.3862	2361.0468	0.0000	1.0000
GTR	43	-1055.9641	2363.4492	0.0000	1.0000
SYMG	41	-1068.1290	2376.0804	0.0000	1.0000
SYMI	41	-1069.8646	2379.5515	0.0000	1.0000
K2P	36	-1085.2221	2381.0199	0.0000	1.0000
SYMIG	42	-1068.6854	2383.0425	0.0000	1.0000
F81G	39	-1088.7360	2405.5956	0.0000	1.0000
SYM	40	-1086.4598	2406.8925	0.0000	1.0000
F81I	39	-1090.4394	2409.0025	0.0000	1.0000
F81IG	40	-1089.1358	2412.2447	0.0000	1.0000
JC69G	36	-1107.6549	2425.8855	0.0000	1.0000
F81	38	-1103.8119	2429.8982	0.0000	1.0000
JC69I	36	-1109.7721	2430.1200	0.0000	1.0000
JC69IG	37	-1108.5832	2433.5914	0.0000	1.0000
JC69	35	-1126.1035	2456.9333	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/chr13_3036.phylip = 1 - 347;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_3036.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/chr13_3036.phylip = 1 - 347;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_3036.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/chr13_3036.phylip = 1 - 347;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_3036.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

