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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-4223.6656	8541.0813	0.8618	0.8618
GTRG	44	-4226.5795	8544.7428	0.1381	1.0000
GTRI	44	-4236.3047	8564.1931	0.0000	1.0000
HKYIG	41	-4241.0373	8567.1828	0.0000	1.0000
HKYG	40	-4243.9498	8570.8573	0.0000	1.0000
SYMIG	42	-4247.6260	8582.5148	0.0000	1.0000
SYMG	41	-4250.9379	8586.9840	0.0000	1.0000
K2PIG	38	-4254.4947	8587.6572	0.0000	1.0000
HKYI	40	-4252.7164	8588.3904	0.0000	1.0000
K2PG	37	-4259.9366	8596.4020	0.0000	1.0000
SYMI	41	-4259.6324	8604.3732	0.0000	1.0000
K2PI	37	-4266.5574	8609.6435	0.0000	1.0000
F81IG	40	-4322.2725	8727.5026	0.0000	1.0000
F81G	39	-4324.9145	8730.6398	0.0000	1.0000
F81I	39	-4332.9232	8746.6573	0.0000	1.0000
JC69IG	37	-4336.8907	8750.3102	0.0000	1.0000
JC69G	36	-4339.6486	8753.6908	0.0000	1.0000
JC69I	36	-4347.6754	8769.7443	0.0000	1.0000
GTR	43	-4362.3400	8814.1013	0.0000	1.0000
HKY	39	-4378.5883	8837.9875	0.0000	1.0000
SYM	40	-4388.5563	8860.0702	0.0000	1.0000
K2P	36	-4397.9198	8870.2331	0.0000	1.0000
F81	38	-4460.3243	8999.3164	0.0000	1.0000
JC69	35	-4476.5273	9025.3168	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-4223.6656	8537.3313	0.8714	0.8714
GTRG	44	-4226.5795	8541.1591	0.1285	1.0000
GTRI	44	-4236.3047	8560.6094	0.0000	1.0000
HKYIG	41	-4241.0373	8564.0745	0.0000	1.0000
HKYG	40	-4243.9498	8567.8996	0.0000	1.0000
SYMIG	42	-4247.6260	8579.2520	0.0000	1.0000
SYMG	41	-4250.9379	8583.8757	0.0000	1.0000
K2PIG	38	-4254.4947	8584.9894	0.0000	1.0000
HKYI	40	-4252.7164	8585.4328	0.0000	1.0000
K2PG	37	-4259.9366	8593.8732	0.0000	1.0000
SYMI	41	-4259.6324	8601.2649	0.0000	1.0000
K2PI	37	-4266.5574	8607.1147	0.0000	1.0000
F81IG	40	-4322.2725	8724.5450	0.0000	1.0000
F81G	39	-4324.9145	8727.8290	0.0000	1.0000
F81I	39	-4332.9232	8743.8464	0.0000	1.0000
JC69IG	37	-4336.8907	8747.7814	0.0000	1.0000
JC69G	36	-4339.6486	8751.2973	0.0000	1.0000
JC69I	36	-4347.6754	8767.3508	0.0000	1.0000
GTR	43	-4362.3400	8810.6799	0.0000	1.0000
HKY	39	-4378.5883	8835.1767	0.0000	1.0000
SYM	40	-4388.5563	8857.1125	0.0000	1.0000
K2P	36	-4397.9198	8867.8396	0.0000	1.0000
F81	38	-4460.3243	8996.6486	0.0000	1.0000
JC69	35	-4476.5273	9023.0547	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-4226.5795	8763.2498	0.6238	0.6238
GTRIG	45	-4223.6656	8764.4696	0.3390	0.9628
HKYG	40	-4243.9498	8769.8003	0.0236	0.9863
HKYIG	41	-4241.0373	8771.0227	0.0128	0.9991
K2PIG	38	-4254.4947	8776.7950	0.0007	0.9999
K2PG	37	-4259.9366	8780.6314	0.0001	1.0000
GTRI	44	-4236.3047	8782.7001	0.0000	1.0000
HKYI	40	-4252.7164	8787.3335	0.0000	1.0000
SYMG	41	-4250.9379	8790.8239	0.0000	1.0000
SYMIG	42	-4247.6260	8791.2477	0.0000	1.0000
K2PI	37	-4266.5574	8793.8729	0.0000	1.0000
SYMI	41	-4259.6324	8808.2131	0.0000	1.0000
F81G	39	-4324.9145	8924.6821	0.0000	1.0000
F81IG	40	-4322.2725	8926.4456	0.0000	1.0000
JC69G	36	-4339.6486	8933.0079	0.0000	1.0000
JC69IG	37	-4336.8907	8934.5396	0.0000	1.0000
F81I	39	-4332.9232	8940.6996	0.0000	1.0000
JC69I	36	-4347.6754	8949.0614	0.0000	1.0000
GTR	43	-4362.3400	9027.7232	0.0000	1.0000
HKY	39	-4378.5883	9032.0298	0.0000	1.0000
K2P	36	-4397.9198	9049.5502	0.0000	1.0000
SYM	40	-4388.5563	9059.0132	0.0000	1.0000
F81	38	-4460.3243	9188.4542	0.0000	1.0000
JC69	35	-4476.5273	9199.7178	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10174.phylip = 1 - 1150;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10174.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10174.phylip = 1 - 1150;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10174.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.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/chr5_10174.phylip = 1 - 1150;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10174.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

