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

Minimum AIC  model: HKYG
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	-1226.8378	2536.3280	0.3701	0.3701
HKYG	40	-1223.5765	2537.3394	0.2232	0.5932
K2PIG	38	-1226.6361	2538.4203	0.1300	0.7232
K2PI	37	-1228.1902	2539.0327	0.0957	0.8190
HKYIG	41	-1223.3739	2539.4769	0.0766	0.8956
HKYI	40	-1224.8702	2539.9267	0.0612	0.9568
GTRG	44	-1221.0826	2542.6180	0.0159	0.9728
SYMG	41	-1225.3294	2543.3877	0.0108	0.9836
GTRIG	45	-1220.9109	2544.8818	0.0051	0.9887
GTRI	44	-1222.2934	2545.0396	0.0047	0.9935
SYMIG	42	-1225.1465	2545.5804	0.0036	0.9971
SYMI	41	-1226.6535	2546.0359	0.0029	1.0000
K2P	36	-1242.1056	2564.3829	0.0000	1.0000
HKY	39	-1239.0340	2565.7275	0.0000	1.0000
GTR	43	-1236.5645	2570.9911	0.0000	1.0000
SYM	40	-1240.4171	2571.0205	0.0000	1.0000
JC69G	36	-1246.9072	2573.9863	0.0000	1.0000
JC69I	36	-1247.9652	2576.1023	0.0000	1.0000
JC69IG	37	-1246.7578	2576.1678	0.0000	1.0000
F81G	39	-1244.4161	2576.4916	0.0000	1.0000
F81I	39	-1245.4449	2578.5492	0.0000	1.0000
F81IG	40	-1244.2711	2578.7285	0.0000	1.0000
JC69	35	-1260.7880	2599.2824	0.0000	1.0000
F81	38	-1258.3216	2601.7913	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1223.5765	2527.1531	0.2830	0.2830
K2PG	37	-1226.8378	2527.6757	0.2179	0.5008
HKYIG	41	-1223.3739	2528.7479	0.1275	0.6283
K2PIG	38	-1226.6361	2529.2722	0.0981	0.7264
HKYI	40	-1224.8702	2529.7403	0.0776	0.8040
GTRG	44	-1221.0826	2530.1652	0.0628	0.8667
K2PI	37	-1228.1902	2530.3804	0.0564	0.9231
GTRIG	45	-1220.9109	2531.8219	0.0274	0.9505
GTRI	44	-1222.2934	2532.5868	0.0187	0.9692
SYMG	41	-1225.3294	2532.6587	0.0180	0.9872
SYMIG	42	-1225.1465	2534.2929	0.0080	0.9952
SYMI	41	-1226.6535	2535.3069	0.0048	1.0000
HKY	39	-1239.0340	2556.0681	0.0000	1.0000
K2P	36	-1242.1056	2556.2112	0.0000	1.0000
GTR	43	-1236.5645	2559.1290	0.0000	1.0000
SYM	40	-1240.4171	2560.8341	0.0000	1.0000
JC69G	36	-1246.9072	2565.8145	0.0000	1.0000
F81G	39	-1244.4161	2566.8322	0.0000	1.0000
JC69IG	37	-1246.7578	2567.5155	0.0000	1.0000
JC69I	36	-1247.9652	2567.9305	0.0000	1.0000
F81IG	40	-1244.2711	2568.5421	0.0000	1.0000
F81I	39	-1245.4449	2568.8897	0.0000	1.0000
JC69	35	-1260.7880	2591.5760	0.0000	1.0000
F81	38	-1258.3216	2592.6431	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1226.8378	2671.7686	0.7530	0.7530
K2PI	37	-1228.1902	2674.4733	0.1948	0.9478
K2PIG	38	-1226.6361	2677.2595	0.0484	0.9962
HKYG	40	-1223.5765	2682.9292	0.0028	0.9990
HKYI	40	-1224.8702	2685.5164	0.0008	0.9998
HKYIG	41	-1223.3739	2688.4184	0.0002	1.0000
SYMG	41	-1225.3294	2692.3292	0.0000	1.0000
SYMI	41	-1226.6535	2694.9774	0.0000	1.0000
K2P	36	-1242.1056	2696.4097	0.0000	1.0000
SYMIG	42	-1225.1465	2697.8578	0.0000	1.0000
GTRG	44	-1221.0826	2701.5189	0.0000	1.0000
GTRI	44	-1222.2934	2703.9405	0.0000	1.0000
JC69G	36	-1246.9072	2706.0130	0.0000	1.0000
GTRIG	45	-1220.9109	2707.0700	0.0000	1.0000
HKY	39	-1239.0340	2707.9498	0.0000	1.0000
JC69I	36	-1247.9652	2708.1290	0.0000	1.0000
JC69IG	37	-1246.7578	2711.6084	0.0000	1.0000
SYM	40	-1240.4171	2716.6103	0.0000	1.0000
F81G	39	-1244.4161	2718.7139	0.0000	1.0000
F81I	39	-1245.4449	2720.7714	0.0000	1.0000
F81IG	40	-1244.2711	2724.3183	0.0000	1.0000
GTR	43	-1236.5645	2726.5883	0.0000	1.0000
JC69	35	-1260.7880	2727.8801	0.0000	1.0000
F81	38	-1258.3216	2740.6304	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/chr5_1597.phylip = 1 - 363;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_1597.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 (K2PG)]

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

