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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	40	-1567.9716	3222.4382	0.5750	0.5750
HKYIG	41	-1567.6076	3224.0486	0.2570	0.8320
HKYI	40	-1569.5575	3225.6101	0.1177	0.9497
GTRG	44	-1566.1974	3228.2989	0.0307	0.9804
GTRIG	45	-1565.8291	3229.9382	0.0135	0.9939
GTRI	44	-1567.8345	3231.5731	0.0060	0.9999
K2PG	37	-1580.8737	3241.2828	0.0000	1.0000
K2PIG	38	-1580.5684	3242.9829	0.0000	1.0000
K2PI	37	-1582.2955	3244.1264	0.0000	1.0000
SYMG	41	-1579.2250	3247.2834	0.0000	1.0000
SYMIG	42	-1578.9450	3249.0708	0.0000	1.0000
SYMI	41	-1580.5719	3249.9771	0.0000	1.0000
HKY	39	-1585.1955	3254.5571	0.0000	1.0000
GTR	43	-1583.3103	3260.1585	0.0000	1.0000
K2P	36	-1597.7609	3272.7557	0.0000	1.0000
F81G	39	-1595.0099	3274.1858	0.0000	1.0000
F81IG	40	-1594.7400	3275.9750	0.0000	1.0000
F81I	39	-1596.5036	3277.1732	0.0000	1.0000
SYM	40	-1595.8878	3278.2707	0.0000	1.0000
JC69G	36	-1605.8431	3288.9200	0.0000	1.0000
JC69IG	37	-1605.5922	3290.7199	0.0000	1.0000
JC69I	36	-1607.2636	3291.7610	0.0000	1.0000
F81	38	-1610.8785	3303.6031	0.0000	1.0000
JC69	35	-1621.6309	3318.2030	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1567.9716	3215.9432	0.5207	0.5207
HKYIG	41	-1567.6076	3217.2152	0.2756	0.7963
HKYI	40	-1569.5575	3219.1151	0.1066	0.9029
GTRG	44	-1566.1974	3220.3947	0.0562	0.9591
GTRIG	45	-1565.8291	3221.6582	0.0299	0.9890
GTRI	44	-1567.8345	3223.6689	0.0109	1.0000
K2PG	37	-1580.8737	3235.7473	0.0000	1.0000
K2PIG	38	-1580.5684	3237.1368	0.0000	1.0000
K2PI	37	-1582.2955	3238.5909	0.0000	1.0000
SYMG	41	-1579.2250	3240.4500	0.0000	1.0000
SYMIG	42	-1578.9450	3241.8899	0.0000	1.0000
SYMI	41	-1580.5719	3243.1438	0.0000	1.0000
HKY	39	-1585.1955	3248.3911	0.0000	1.0000
GTR	43	-1583.3103	3252.6207	0.0000	1.0000
K2P	36	-1597.7609	3267.5219	0.0000	1.0000
F81G	39	-1595.0099	3268.0198	0.0000	1.0000
F81IG	40	-1594.7400	3269.4800	0.0000	1.0000
F81I	39	-1596.5036	3271.0072	0.0000	1.0000
SYM	40	-1595.8878	3271.7757	0.0000	1.0000
JC69G	36	-1605.8431	3283.6862	0.0000	1.0000
JC69IG	37	-1605.5922	3285.1845	0.0000	1.0000
JC69I	36	-1607.2636	3286.5272	0.0000	1.0000
F81	38	-1610.8785	3297.7569	0.0000	1.0000
JC69	35	-1621.6309	3313.2618	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1567.9716	3388.0480	0.7647	0.7647
HKYI	40	-1569.5575	3391.2198	0.1566	0.9213
HKYIG	41	-1567.6076	3393.6226	0.0471	0.9684
K2PG	37	-1580.8737	3394.9442	0.0243	0.9927
K2PI	37	-1582.2955	3397.7878	0.0059	0.9986
K2PIG	38	-1580.5684	3400.6363	0.0014	1.0000
GTRG	44	-1566.1974	3409.7100	0.0000	1.0000
GTRI	44	-1567.8345	3412.9842	0.0000	1.0000
GTRIG	45	-1565.8291	3415.2760	0.0000	1.0000
HKY	39	-1585.1955	3416.1932	0.0000	1.0000
SYMG	41	-1579.2250	3416.8574	0.0000	1.0000
SYMI	41	-1580.5719	3419.5512	0.0000	1.0000
K2P	36	-1597.7609	3422.4161	0.0000	1.0000
SYMIG	42	-1578.9450	3422.5999	0.0000	1.0000
F81G	39	-1595.0099	3435.8220	0.0000	1.0000
GTR	43	-1583.3103	3437.6333	0.0000	1.0000
JC69G	36	-1605.8431	3438.5805	0.0000	1.0000
F81I	39	-1596.5036	3438.8093	0.0000	1.0000
JC69I	36	-1607.2636	3441.4215	0.0000	1.0000
F81IG	40	-1594.7400	3441.5847	0.0000	1.0000
SYM	40	-1595.8878	3443.8804	0.0000	1.0000
JC69IG	37	-1605.5922	3444.3814	0.0000	1.0000
F81	38	-1610.8785	3461.2564	0.0000	1.0000
JC69	35	-1621.6309	3463.8534	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/chr6_3293.phylip = 1 - 546;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3293.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 (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3293.phylip = 1 - 546;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3293.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 BIC model (HKYG)]

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


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

