
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:44 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_22616.phylip" (19 taxa, 423 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	-1368.1298	2824.8460	0.5562	0.5562
HKYIG	41	-1367.6318	2826.3031	0.2684	0.8246
HKYI	40	-1369.3994	2827.3852	0.1563	0.9809
GTRG	44	-1367.0980	2832.6721	0.0111	0.9920
GTRIG	45	-1366.6315	2834.2444	0.0051	0.9971
GTRI	44	-1368.4572	2835.3906	0.0029	0.9999
K2PG	37	-1381.4645	2844.2328	0.0000	1.0000
K2PIG	38	-1380.9260	2845.5707	0.0000	1.0000
K2PI	37	-1382.4509	2846.2058	0.0000	1.0000
SYMG	41	-1380.6742	2852.3879	0.0000	1.0000
SYMIG	42	-1380.1501	2853.8056	0.0000	1.0000
SYMI	41	-1381.5989	2854.2372	0.0000	1.0000
HKY	39	-1389.3412	2864.8286	0.0000	1.0000
GTR	43	-1388.0705	2872.1252	0.0000	1.0000
K2P	36	-1403.2235	2885.3485	0.0000	1.0000
SYM	40	-1402.3335	2893.2533	0.0000	1.0000
F81G	39	-1413.9091	2913.9644	0.0000	1.0000
F81IG	40	-1413.3891	2915.3647	0.0000	1.0000
F81I	39	-1415.3103	2916.7668	0.0000	1.0000
JC69G	36	-1425.0248	2928.9512	0.0000	1.0000
JC69IG	37	-1424.4903	2930.2844	0.0000	1.0000
JC69I	36	-1426.3128	2931.5272	0.0000	1.0000
F81	38	-1434.7643	2953.2474	0.0000	1.0000
JC69	35	-1446.2637	2969.0390	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1368.1298	2816.2597	0.5048	0.5048
HKYIG	41	-1367.6318	2817.2637	0.3055	0.8103
HKYI	40	-1369.3994	2818.7988	0.1418	0.9521
GTRG	44	-1367.0980	2822.1959	0.0259	0.9781
GTRIG	45	-1366.6315	2823.2630	0.0152	0.9933
GTRI	44	-1368.4572	2824.9144	0.0067	1.0000
K2PG	37	-1381.4645	2836.9289	0.0000	1.0000
K2PIG	38	-1380.9260	2837.8519	0.0000	1.0000
K2PI	37	-1382.4509	2838.9019	0.0000	1.0000
SYMG	41	-1380.6742	2843.3485	0.0000	1.0000
SYMIG	42	-1380.1501	2844.3003	0.0000	1.0000
SYMI	41	-1381.5989	2845.1978	0.0000	1.0000
HKY	39	-1389.3412	2856.6823	0.0000	1.0000
GTR	43	-1388.0705	2862.1410	0.0000	1.0000
K2P	36	-1403.2235	2878.4469	0.0000	1.0000
SYM	40	-1402.3335	2884.6669	0.0000	1.0000
F81G	39	-1413.9091	2905.8182	0.0000	1.0000
F81IG	40	-1413.3891	2906.7783	0.0000	1.0000
F81I	39	-1415.3103	2908.6206	0.0000	1.0000
JC69G	36	-1425.0248	2922.0496	0.0000	1.0000
JC69IG	37	-1424.4903	2922.9805	0.0000	1.0000
JC69I	36	-1426.3128	2924.6257	0.0000	1.0000
F81	38	-1434.7643	2945.5286	0.0000	1.0000
JC69	35	-1446.2637	2962.5274	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1368.1298	2978.1545	0.7239	0.7239
HKYI	40	-1369.3994	2980.6937	0.2034	0.9272
HKYIG	41	-1367.6318	2983.2059	0.0579	0.9852
K2PG	37	-1381.4645	2986.6817	0.0102	0.9953
K2PI	37	-1382.4509	2988.6546	0.0038	0.9991
K2PIG	38	-1380.9260	2991.6521	0.0008	1.0000
GTRG	44	-1367.0980	3000.2803	0.0000	1.0000
GTRI	44	-1368.4572	3002.9987	0.0000	1.0000
GTRIG	45	-1366.6315	3005.3947	0.0000	1.0000
SYMG	41	-1380.6742	3009.2907	0.0000	1.0000
SYMI	41	-1381.5989	3011.1401	0.0000	1.0000
SYMIG	42	-1380.1501	3014.2899	0.0000	1.0000
HKY	39	-1389.3412	3014.5299	0.0000	1.0000
K2P	36	-1403.2235	3024.1523	0.0000	1.0000
GTR	43	-1388.0705	3036.1780	0.0000	1.0000
SYM	40	-1402.3335	3046.5618	0.0000	1.0000
F81G	39	-1413.9091	3063.6657	0.0000	1.0000
F81I	39	-1415.3103	3066.4681	0.0000	1.0000
JC69G	36	-1425.0248	3067.7550	0.0000	1.0000
F81IG	40	-1413.3891	3068.6732	0.0000	1.0000
JC69I	36	-1426.3128	3070.3311	0.0000	1.0000
JC69IG	37	-1424.4903	3072.7333	0.0000	1.0000
F81	38	-1434.7643	3099.3287	0.0000	1.0000
JC69	35	-1446.2637	3104.1854	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_22616.phylip = 1 - 423;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22616.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_22616.phylip = 1 - 423;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22616.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_22616.phylip = 1 - 423;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22616.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

