
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:42 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_25261.phylip" (19 taxa, 348 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	-1264.1275	2618.9390	0.7434	0.7434
HKYIG	41	-1264.5537	2622.3622	0.1342	0.8777
HKYI	40	-1266.3944	2623.4728	0.0770	0.9547
GTRG	44	-1262.1833	2625.4359	0.0289	0.9836
GTRIG	45	-1261.9327	2627.5740	0.0099	0.9935
GTRI	44	-1263.6778	2628.4248	0.0065	1.0000
K2PG	37	-1279.6488	2642.3685	0.0000	1.0000
K2PIG	38	-1279.3095	2644.2113	0.0000	1.0000
K2PI	37	-1280.7159	2644.5028	0.0000	1.0000
SYMG	41	-1279.2578	2651.7705	0.0000	1.0000
SYMI	41	-1280.2124	2653.6796	0.0000	1.0000
SYMIG	42	-1278.9519	2653.7463	0.0000	1.0000
HKY	39	-1284.1772	2656.4842	0.0000	1.0000
GTR	43	-1282.8925	2664.2324	0.0000	1.0000
K2P	36	-1298.3388	2677.2435	0.0000	1.0000
SYM	40	-1297.7639	2686.2118	0.0000	1.0000
F81G	39	-1308.3062	2704.7423	0.0000	1.0000
F81IG	40	-1309.1932	2709.0704	0.0000	1.0000
F81I	39	-1310.8104	2709.7507	0.0000	1.0000
JC69G	36	-1319.0927	2718.7512	0.0000	1.0000
JC69I	36	-1319.9722	2720.5103	0.0000	1.0000
JC69IG	37	-1318.8463	2720.7636	0.0000	1.0000
F81	38	-1326.1868	2737.9658	0.0000	1.0000
JC69	35	-1336.3859	2750.8487	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1264.1275	2608.2550	0.6406	0.6406
HKYIG	41	-1264.5537	2611.1073	0.1539	0.7945
GTRG	44	-1262.1833	2612.3666	0.0820	0.8765
HKYI	40	-1266.3944	2612.7888	0.0664	0.9428
GTRIG	45	-1261.9327	2613.8654	0.0388	0.9816
GTRI	44	-1263.6778	2615.3555	0.0184	1.0000
K2PG	37	-1279.6488	2633.2976	0.0000	1.0000
K2PIG	38	-1279.3095	2634.6191	0.0000	1.0000
K2PI	37	-1280.7159	2635.4318	0.0000	1.0000
SYMG	41	-1279.2578	2640.5156	0.0000	1.0000
SYMIG	42	-1278.9519	2641.9037	0.0000	1.0000
SYMI	41	-1280.2124	2642.4247	0.0000	1.0000
HKY	39	-1284.1772	2646.3544	0.0000	1.0000
GTR	43	-1282.8925	2651.7850	0.0000	1.0000
K2P	36	-1298.3388	2668.6776	0.0000	1.0000
SYM	40	-1297.7639	2675.5277	0.0000	1.0000
F81G	39	-1308.3062	2694.6124	0.0000	1.0000
F81IG	40	-1309.1932	2698.3864	0.0000	1.0000
F81I	39	-1310.8104	2699.6208	0.0000	1.0000
JC69G	36	-1319.0927	2710.1853	0.0000	1.0000
JC69IG	37	-1318.8463	2711.6926	0.0000	1.0000
JC69I	36	-1319.9722	2711.9444	0.0000	1.0000
F81	38	-1326.1868	2728.3736	0.0000	1.0000
JC69	35	-1336.3859	2742.7718	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1264.1275	2762.3431	0.8769	0.8769
HKYI	40	-1266.3944	2766.8769	0.0909	0.9678
HKYIG	41	-1264.5537	2769.0476	0.0307	0.9985
K2PG	37	-1279.6488	2775.8291	0.0010	0.9995
K2PI	37	-1280.7159	2777.9633	0.0004	0.9999
K2PIG	38	-1279.3095	2781.0028	0.0001	0.9999
GTRG	44	-1262.1833	2781.8635	0.0001	1.0000
GTRI	44	-1263.6778	2784.8524	0.0000	1.0000
GTRIG	45	-1261.9327	2787.2145	0.0000	1.0000
HKY	39	-1284.1772	2796.5903	0.0000	1.0000
SYMG	41	-1279.2578	2798.4559	0.0000	1.0000
SYMI	41	-1280.2124	2800.3650	0.0000	1.0000
SYMIG	42	-1278.9519	2803.6962	0.0000	1.0000
K2P	36	-1298.3388	2807.3569	0.0000	1.0000
GTR	43	-1282.8925	2817.4297	0.0000	1.0000
SYM	40	-1297.7639	2829.6158	0.0000	1.0000
F81G	39	-1308.3062	2844.8483	0.0000	1.0000
JC69G	36	-1319.0927	2848.8646	0.0000	1.0000
F81I	39	-1310.8104	2849.8567	0.0000	1.0000
JC69I	36	-1319.9722	2850.6237	0.0000	1.0000
F81IG	40	-1309.1932	2852.4745	0.0000	1.0000
JC69IG	37	-1318.8463	2854.2241	0.0000	1.0000
F81	38	-1326.1868	2874.7573	0.0000	1.0000
JC69	35	-1336.3859	2877.5988	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/chr2_25261.phylip = 1 - 348;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_25261.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/chr2_25261.phylip = 1 - 348;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_25261.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/chr2_25261.phylip = 1 - 348;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_25261.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

