
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:02 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/chr6_3033.phylip" (19 taxa, 372 characters)

Minimum AIC  model: SYMI
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	37	-1401.4578	2885.3347	0.2271	0.2271
K2PG	37	-1401.6683	2885.7558	0.1839	0.4110
SYMI	41	-1396.7504	2885.9372	0.1680	0.5790
SYMG	41	-1396.8992	2886.2348	0.1448	0.7238
K2PIG	38	-1401.4008	2887.7024	0.0695	0.7933
SYMIG	42	-1396.5330	2888.0447	0.0586	0.8519
GTRI	44	-1394.1000	2888.3101	0.0513	0.9032
GTRG	44	-1394.1585	2888.4272	0.0484	0.9515
GTRIG	45	-1393.6970	2890.0933	0.0210	0.9726
HKYI	40	-1400.5775	2891.0645	0.0129	0.9855
HKYG	40	-1400.7783	2891.4660	0.0106	0.9961
HKYIG	41	-1400.5123	2893.4610	0.0039	1.0000
K2P	36	-1419.8445	2919.6411	0.0000	1.0000
SYM	40	-1415.1266	2920.1625	0.0000	1.0000
GTR	43	-1412.5878	2922.7122	0.0000	1.0000
HKY	39	-1418.9687	2925.3349	0.0000	1.0000
JC69I	36	-1442.8007	2965.5536	0.0000	1.0000
JC69G	36	-1443.4719	2966.8961	0.0000	1.0000
JC69IG	37	-1442.8089	2968.0369	0.0000	1.0000
F81I	39	-1442.3079	2972.0134	0.0000	1.0000
F81G	39	-1442.9947	2973.3869	0.0000	1.0000
F81IG	40	-1442.3169	2974.5432	0.0000	1.0000
JC69	35	-1460.1293	2997.7586	0.0000	1.0000
F81	38	-1459.6236	3004.1481	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMI	41	-1396.7504	2875.5009	0.1890	0.1890
SYMG	41	-1396.8992	2875.7984	0.1628	0.3518
GTRI	44	-1394.1000	2876.2000	0.1332	0.4850
GTRG	44	-1394.1585	2876.3171	0.1256	0.6107
K2PI	37	-1401.4578	2876.9155	0.0932	0.7038
SYMIG	42	-1396.5330	2877.0659	0.0864	0.7902
K2PG	37	-1401.6683	2877.3367	0.0755	0.8657
GTRIG	45	-1393.6970	2877.3939	0.0733	0.9390
K2PIG	38	-1401.4008	2878.8015	0.0363	0.9753
HKYI	40	-1400.5775	2881.1551	0.0112	0.9865
HKYG	40	-1400.7783	2881.5567	0.0091	0.9956
HKYIG	41	-1400.5123	2883.0247	0.0044	1.0000
SYM	40	-1415.1266	2910.2531	0.0000	1.0000
GTR	43	-1412.5878	2911.1756	0.0000	1.0000
K2P	36	-1419.8445	2911.6889	0.0000	1.0000
HKY	39	-1418.9687	2915.9373	0.0000	1.0000
JC69I	36	-1442.8007	2957.6013	0.0000	1.0000
JC69G	36	-1443.4719	2958.9438	0.0000	1.0000
JC69IG	37	-1442.8089	2959.6177	0.0000	1.0000
F81I	39	-1442.3079	2962.6158	0.0000	1.0000
F81G	39	-1442.9947	2963.9893	0.0000	1.0000
F81IG	40	-1442.3169	2964.6338	0.0000	1.0000
JC69	35	-1460.1293	2990.2586	0.0000	1.0000
F81	38	-1459.6236	2995.2472	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-1401.4578	3021.9146	0.5356	0.5356
K2PG	37	-1401.6683	3022.3358	0.4339	0.9694
K2PIG	38	-1401.4008	3027.7195	0.0294	0.9988
SYMI	41	-1396.7504	3036.1755	0.0004	0.9993
SYMG	41	-1396.8992	3036.4730	0.0004	0.9996
HKYI	40	-1400.5775	3037.9109	0.0002	0.9998
HKYG	40	-1400.7783	3038.3124	0.0001	1.0000
SYMIG	42	-1396.5330	3041.6595	0.0000	1.0000
HKYIG	41	-1400.5123	3043.6993	0.0000	1.0000
GTRI	44	-1394.1000	3048.6313	0.0000	1.0000
GTRG	44	-1394.1585	3048.7484	0.0000	1.0000
K2P	36	-1419.8445	3052.7691	0.0000	1.0000
GTRIG	45	-1393.6970	3053.7442	0.0000	1.0000
SYM	40	-1415.1266	3067.0089	0.0000	1.0000
HKY	39	-1418.9687	3068.7742	0.0000	1.0000
GTR	43	-1412.5878	3079.6881	0.0000	1.0000
JC69I	36	-1442.8007	3098.6815	0.0000	1.0000
JC69G	36	-1443.4719	3100.0240	0.0000	1.0000
JC69IG	37	-1442.8089	3104.6168	0.0000	1.0000
F81I	39	-1442.3079	3115.4527	0.0000	1.0000
F81G	39	-1442.9947	3116.8262	0.0000	1.0000
F81IG	40	-1442.3169	3121.3896	0.0000	1.0000
JC69	35	-1460.1293	3127.4199	0.0000	1.0000
F81	38	-1459.6236	3144.1652	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3033.phylip = 1 - 372;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3033.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3033.phylip = 1 - 372;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3033.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (K2PI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3033.phylip = 1 - 372;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3033.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

