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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	44	-1578.0480	3250.3224	0.7853	0.7853
GTRIG	45	-1578.2205	3252.9606	0.2100	0.9952
GTRG	44	-1583.1535	3260.5335	0.0048	1.0000
SYMIG	42	-1591.2812	3272.2237	0.0000	1.0000
SYMI	41	-1593.2759	3273.9415	0.0000	1.0000
HKYIG	41	-1593.4079	3274.2054	0.0000	1.0000
HKYI	40	-1595.8772	3276.8793	0.0000	1.0000
SYMG	41	-1599.0439	3285.4774	0.0000	1.0000
K2PI	37	-1603.5567	3285.4866	0.0000	1.0000
K2PIG	38	-1604.0705	3288.7578	0.0000	1.0000
HKYG	40	-1604.8477	3294.8203	0.0000	1.0000
K2PG	37	-1612.5007	3303.3746	0.0000	1.0000
F81IG	40	-1616.7194	3318.5639	0.0000	1.0000
F81I	39	-1618.9678	3320.8029	0.0000	1.0000
JC69IG	37	-1623.1318	3324.6368	0.0000	1.0000
JC69I	36	-1625.4528	3327.0423	0.0000	1.0000
F81G	39	-1627.7808	3338.4290	0.0000	1.0000
JC69G	36	-1634.0416	3344.2199	0.0000	1.0000
GTR	43	-1635.5233	3362.9870	0.0000	1.0000
SYM	40	-1650.2376	3385.6003	0.0000	1.0000
HKY	39	-1654.7647	3392.3969	0.0000	1.0000
K2P	36	-1662.3530	3400.8427	0.0000	1.0000
F81	38	-1676.6194	3433.8557	0.0000	1.0000
JC69	35	-1682.9077	3439.7223	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	44	-1578.0480	3244.0960	0.7601	0.7601
GTRIG	45	-1578.2205	3246.4409	0.2353	0.9954
GTRG	44	-1583.1535	3254.3071	0.0046	1.0000
SYMIG	42	-1591.2812	3266.5623	0.0000	1.0000
SYMI	41	-1593.2759	3268.5519	0.0000	1.0000
HKYIG	41	-1593.4079	3268.8157	0.0000	1.0000
HKYI	40	-1595.8772	3271.7543	0.0000	1.0000
SYMG	41	-1599.0439	3280.0877	0.0000	1.0000
K2PI	37	-1603.5567	3281.1134	0.0000	1.0000
K2PIG	38	-1604.0705	3284.1410	0.0000	1.0000
HKYG	40	-1604.8477	3289.6953	0.0000	1.0000
K2PG	37	-1612.5007	3299.0013	0.0000	1.0000
F81IG	40	-1616.7194	3313.4389	0.0000	1.0000
F81I	39	-1618.9678	3315.9355	0.0000	1.0000
JC69IG	37	-1623.1318	3320.2635	0.0000	1.0000
JC69I	36	-1625.4528	3322.9056	0.0000	1.0000
F81G	39	-1627.7808	3333.5616	0.0000	1.0000
JC69G	36	-1634.0416	3340.0833	0.0000	1.0000
GTR	43	-1635.5233	3357.0467	0.0000	1.0000
SYM	40	-1650.2376	3380.4753	0.0000	1.0000
HKY	39	-1654.7647	3387.5295	0.0000	1.0000
K2P	36	-1662.3530	3396.7060	0.0000	1.0000
F81	38	-1676.6194	3429.2389	0.0000	1.0000
JC69	35	-1682.9077	3435.8154	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	44	-1578.0480	3443.1328	0.8877	0.8877
K2PI	37	-1603.5567	3448.4852	0.0611	0.9488
GTRIG	45	-1578.2205	3450.0012	0.0286	0.9775
HKYI	40	-1595.8772	3452.6968	0.0074	0.9849
GTRG	44	-1583.1535	3453.3438	0.0054	0.9903
SYMI	41	-1593.2759	3454.0179	0.0038	0.9941
HKYIG	41	-1593.4079	3454.2818	0.0034	0.9975
K2PIG	38	-1604.0705	3456.0363	0.0014	0.9989
SYMIG	42	-1591.2812	3456.5519	0.0011	1.0000
SYMG	41	-1599.0439	3465.5538	0.0000	1.0000
K2PG	37	-1612.5007	3466.3731	0.0000	1.0000
HKYG	40	-1604.8477	3470.6378	0.0000	1.0000
JC69I	36	-1625.4528	3485.7539	0.0000	1.0000
JC69IG	37	-1623.1318	3487.6353	0.0000	1.0000
F81I	39	-1618.9678	3492.3544	0.0000	1.0000
F81IG	40	-1616.7194	3494.3814	0.0000	1.0000
JC69G	36	-1634.0416	3502.9315	0.0000	1.0000
F81G	39	-1627.7808	3509.9806	0.0000	1.0000
GTR	43	-1635.5233	3551.5599	0.0000	1.0000
K2P	36	-1662.3530	3559.5543	0.0000	1.0000
SYM	40	-1650.2376	3561.4178	0.0000	1.0000
HKY	39	-1654.7647	3563.9484	0.0000	1.0000
JC69	35	-1682.9077	3594.1400	0.0000	1.0000
F81	38	-1676.6194	3601.1342	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4994.phylip = 1 - 681;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4994.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=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4994.phylip = 1 - 681;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4994.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=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4994.phylip = 1 - 681;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4994.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=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

