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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-1424.1528	2948.2334	0.4632	0.4632
GTRIG	45	-1423.9646	2950.4367	0.1539	0.6171
HKYG	40	-1430.3789	2950.5197	0.1477	0.7647
GTRI	44	-1425.3307	2950.5892	0.1426	0.9074
HKYIG	41	-1430.1245	2952.5295	0.0541	0.9614
HKYI	40	-1432.0628	2953.8875	0.0274	0.9888
K2PG	37	-1437.2921	2956.8793	0.0061	0.9950
K2PIG	38	-1437.1016	2958.9723	0.0022	0.9971
K2PI	37	-1438.7784	2959.8518	0.0014	0.9985
SYMG	41	-1434.2465	2960.7736	0.0009	0.9994
SYMI	41	-1435.2509	2962.7824	0.0003	0.9997
SYMIG	42	-1434.0980	2963.0103	0.0003	1.0000
HKY	39	-1451.0337	2989.3256	0.0000	1.0000
GTR	43	-1446.3906	2990.1446	0.0000	1.0000
K2P	36	-1462.1825	3004.2002	0.0000	1.0000
SYM	40	-1459.1827	3008.1272	0.0000	1.0000
F81G	39	-1469.8911	3027.0403	0.0000	1.0000
F81I	39	-1470.7766	3028.8113	0.0000	1.0000
F81IG	40	-1469.7660	3029.2938	0.0000	1.0000
JC69G	36	-1474.7838	3029.4029	0.0000	1.0000
JC69I	36	-1475.8506	3031.5365	0.0000	1.0000
JC69IG	37	-1474.6655	3031.6260	0.0000	1.0000
F81	38	-1489.0385	3062.8462	0.0000	1.0000
JC69	35	-1499.3226	3076.0352	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1424.1528	2936.3057	0.5164	0.5164
GTRIG	45	-1423.9646	2937.9291	0.2293	0.7458
GTRI	44	-1425.3307	2938.6615	0.1590	0.9048
HKYG	40	-1430.3789	2940.7578	0.0558	0.9606
HKYIG	41	-1430.1245	2942.2489	0.0265	0.9870
HKYI	40	-1432.0628	2944.1256	0.0104	0.9974
K2PG	37	-1437.2921	2948.5843	0.0011	0.9985
K2PIG	38	-1437.1016	2950.2031	0.0005	0.9990
SYMG	41	-1434.2465	2950.4930	0.0004	0.9994
K2PI	37	-1438.7784	2951.5569	0.0003	0.9997
SYMIG	42	-1434.0980	2952.1959	0.0002	0.9998
SYMI	41	-1435.2509	2952.5018	0.0002	1.0000
GTR	43	-1446.3906	2978.7812	0.0000	1.0000
HKY	39	-1451.0337	2980.0675	0.0000	1.0000
K2P	36	-1462.1825	2996.3649	0.0000	1.0000
SYM	40	-1459.1827	2998.3653	0.0000	1.0000
F81G	39	-1469.8911	3017.7821	0.0000	1.0000
F81IG	40	-1469.7660	3019.5319	0.0000	1.0000
F81I	39	-1470.7766	3019.5531	0.0000	1.0000
JC69G	36	-1474.7838	3021.5676	0.0000	1.0000
JC69IG	37	-1474.6655	3023.3310	0.0000	1.0000
JC69I	36	-1475.8506	3023.7012	0.0000	1.0000
F81	38	-1489.0385	3054.0769	0.0000	1.0000
JC69	35	-1499.3226	3068.6452	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1437.2921	3094.0774	0.6843	0.6843
K2PI	37	-1438.7784	3097.0499	0.1548	0.8391
HKYG	40	-1430.3789	3098.0476	0.0940	0.9331
K2PIG	38	-1437.1016	3099.6284	0.0426	0.9757
HKYI	40	-1432.0628	3101.4154	0.0174	0.9932
HKYIG	41	-1430.1245	3103.4710	0.0062	0.9994
GTRG	44	-1424.1528	3109.3245	0.0003	0.9997
GTRI	44	-1425.3307	3111.6802	0.0001	0.9998
SYMG	41	-1434.2465	3111.7150	0.0001	0.9999
SYMI	41	-1435.2509	3113.7239	0.0000	1.0000
GTRIG	45	-1423.9646	3114.8802	0.0000	1.0000
SYMIG	42	-1434.0980	3117.3502	0.0000	1.0000
HKY	39	-1451.0337	3133.4250	0.0000	1.0000
K2P	36	-1462.1825	3137.9257	0.0000	1.0000
GTR	43	-1446.3906	3147.8678	0.0000	1.0000
SYM	40	-1459.1827	3155.6551	0.0000	1.0000
JC69G	36	-1474.7838	3163.1284	0.0000	1.0000
JC69I	36	-1475.8506	3165.2620	0.0000	1.0000
JC69IG	37	-1474.6655	3168.8241	0.0000	1.0000
F81G	39	-1469.8911	3171.1397	0.0000	1.0000
F81I	39	-1470.7766	3172.9107	0.0000	1.0000
F81IG	40	-1469.7660	3176.8217	0.0000	1.0000
F81	38	-1489.0385	3203.5022	0.0000	1.0000
JC69	35	-1499.3226	3206.2738	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/chr4_13410.phylip = 1 - 377;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_13410.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/chr4_13410.phylip = 1 - 377;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_13410.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 (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_13410.phylip = 1 - 377;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_13410.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

