
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:23 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/chr1_34724.phylip" (19 taxa, 463 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1455.7698	3009.0132	0.5401	0.5401
GTRI	44	-1456.3659	3010.2055	0.2975	0.8376
GTRIG	45	-1455.7442	3011.4165	0.1624	1.0000
HKYG	40	-1473.5946	3034.9618	0.0000	1.0000
HKYI	40	-1474.4486	3036.6698	0.0000	1.0000
HKYIG	41	-1473.4589	3037.0982	0.0000	1.0000
GTR	43	-1475.9733	3046.9775	0.0000	1.0000
SYMG	41	-1484.0526	3058.2857	0.0000	1.0000
SYMI	41	-1484.7965	3059.7736	0.0000	1.0000
SYMIG	42	-1484.0022	3060.6044	0.0000	1.0000
HKY	39	-1491.9686	3069.3131	0.0000	1.0000
K2PG	37	-1495.6198	3071.8561	0.0000	1.0000
K2PIG	38	-1495.5114	3074.0133	0.0000	1.0000
K2PI	37	-1496.7760	3074.1685	0.0000	1.0000
SYM	40	-1500.4282	3088.6289	0.0000	1.0000
K2P	36	-1512.2536	3102.7606	0.0000	1.0000
F81G	39	-1521.3441	3128.0640	0.0000	1.0000
F81IG	40	-1521.1399	3130.0524	0.0000	1.0000
F81I	39	-1522.3845	3130.1448	0.0000	1.0000
JC69G	36	-1536.6876	3151.6287	0.0000	1.0000
JC69IG	37	-1536.4995	3153.6155	0.0000	1.0000
JC69I	36	-1537.8326	3153.9188	0.0000	1.0000
F81	38	-1538.2021	3159.3948	0.0000	1.0000
JC69	35	-1553.0578	3182.0172	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1455.7698	2999.5395	0.5186	0.5186
GTRI	44	-1456.3659	3000.7319	0.2857	0.8043
GTRIG	45	-1455.7442	3001.4884	0.1957	1.0000
HKYG	40	-1473.5946	3027.1893	0.0000	1.0000
HKYI	40	-1474.4486	3028.8973	0.0000	1.0000
HKYIG	41	-1473.4589	3028.9177	0.0000	1.0000
GTR	43	-1475.9733	3037.9465	0.0000	1.0000
SYMG	41	-1484.0526	3050.1051	0.0000	1.0000
SYMI	41	-1484.7965	3051.5931	0.0000	1.0000
SYMIG	42	-1484.0022	3052.0044	0.0000	1.0000
HKY	39	-1491.9686	3061.9372	0.0000	1.0000
K2PG	37	-1495.6198	3065.2396	0.0000	1.0000
K2PIG	38	-1495.5114	3067.0227	0.0000	1.0000
K2PI	37	-1496.7760	3067.5521	0.0000	1.0000
SYM	40	-1500.4282	3080.8563	0.0000	1.0000
K2P	36	-1512.2536	3096.5071	0.0000	1.0000
F81G	39	-1521.3441	3120.6882	0.0000	1.0000
F81IG	40	-1521.1399	3122.2799	0.0000	1.0000
F81I	39	-1522.3845	3122.7689	0.0000	1.0000
JC69G	36	-1536.6876	3145.3752	0.0000	1.0000
JC69IG	37	-1536.4995	3146.9990	0.0000	1.0000
JC69I	36	-1537.8326	3147.6652	0.0000	1.0000
F81	38	-1538.2021	3152.4043	0.0000	1.0000
JC69	35	-1553.0578	3176.1156	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1455.7698	3181.5995	0.6233	0.6233
GTRI	44	-1456.3659	3182.7919	0.3434	0.9667
GTRIG	45	-1455.7442	3187.6861	0.0297	0.9964
HKYG	40	-1473.5946	3192.6983	0.0024	0.9988
HKYI	40	-1474.4486	3194.4063	0.0010	0.9999
HKYIG	41	-1473.4589	3198.5645	0.0001	1.0000
GTR	43	-1475.9733	3215.8688	0.0000	1.0000
K2PG	37	-1495.6198	3218.3355	0.0000	1.0000
SYMG	41	-1484.0526	3219.7519	0.0000	1.0000
K2PI	37	-1496.7760	3220.6480	0.0000	1.0000
SYMI	41	-1484.7965	3221.2399	0.0000	1.0000
HKY	39	-1491.9686	3223.3086	0.0000	1.0000
K2PIG	38	-1495.5114	3224.2563	0.0000	1.0000
SYMIG	42	-1484.0022	3225.7889	0.0000	1.0000
K2P	36	-1512.2536	3245.4653	0.0000	1.0000
SYM	40	-1500.4282	3246.3654	0.0000	1.0000
F81G	39	-1521.3441	3282.0595	0.0000	1.0000
F81I	39	-1522.3845	3284.1403	0.0000	1.0000
F81IG	40	-1521.1399	3287.7889	0.0000	1.0000
JC69G	36	-1536.6876	3294.3334	0.0000	1.0000
JC69I	36	-1537.8326	3296.6234	0.0000	1.0000
JC69IG	37	-1536.4995	3300.0949	0.0000	1.0000
F81	38	-1538.2021	3309.6379	0.0000	1.0000
JC69	35	-1553.0578	3320.9360	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/chr1_34724.phylip = 1 - 463;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_34724.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/chr1_34724.phylip = 1 - 463;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_34724.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 (GTRG)]

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


-------------------------------------------------------------
End of Output

