
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:25 04/04/2011

Input data from file "591-tomodel-phylip/chr13_5324.phylip" (7 taxa, 659 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1632.5370	3306.3905	0.4182	0.4182
GTRI	20	-1632.9017	3307.1199	0.2904	0.7086
GTRIG	21	-1632.5279	3308.5064	0.1452	0.8538
HKYG	16	-1638.5891	3310.0256	0.0679	0.9217
HKYI	16	-1638.8126	3310.4726	0.0543	0.9760
HKYIG	17	-1638.5866	3312.1280	0.0237	0.9998
GTR	19	-1641.2039	3321.5972	0.0002	1.0000
HKY	15	-1647.7802	3326.3070	0.0000	1.0000
F81G	15	-1674.4338	3379.6142	0.0000	1.0000
F81I	15	-1674.5627	3379.8719	0.0000	1.0000
F81IG	16	-1674.4374	3381.7222	0.0000	1.0000
F81	14	-1683.8233	3396.2988	0.0000	1.0000
K2PG	13	-1701.9422	3430.4487	0.0000	1.0000
K2PI	13	-1701.9725	3430.5094	0.0000	1.0000
K2PIG	14	-1701.9389	3432.5299	0.0000	1.0000
SYMG	17	-1700.5249	3436.0045	0.0000	1.0000
SYMI	17	-1700.6418	3436.2383	0.0000	1.0000
SYMIG	18	-1700.5130	3438.0948	0.0000	1.0000
K2P	12	-1712.2081	3448.8993	0.0000	1.0000
SYM	16	-1710.2196	3453.2865	0.0000	1.0000
JC69G	12	-1729.4123	3483.3076	0.0000	1.0000
JC69I	12	-1729.4526	3483.3883	0.0000	1.0000
JC69IG	13	-1729.4174	3485.3991	0.0000	1.0000
JC69	11	-1739.6150	3501.6380	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1632.5370	3305.0739	0.4265	0.4265
GTRI	20	-1632.9017	3305.8033	0.2962	0.7227
GTRIG	21	-1632.5279	3307.0559	0.1583	0.8810
HKYG	16	-1638.5891	3309.1782	0.0548	0.9358
HKYI	16	-1638.8126	3309.6252	0.0438	0.9796
HKYIG	17	-1638.5866	3311.1733	0.0202	0.9998
GTR	19	-1641.2039	3320.4079	0.0002	1.0000
HKY	15	-1647.7802	3325.5605	0.0000	1.0000
F81G	15	-1674.4338	3378.8677	0.0000	1.0000
F81I	15	-1674.5627	3379.1254	0.0000	1.0000
F81IG	16	-1674.4374	3380.8749	0.0000	1.0000
F81	14	-1683.8233	3395.6466	0.0000	1.0000
K2PG	13	-1701.9422	3429.8844	0.0000	1.0000
K2PI	13	-1701.9725	3429.9450	0.0000	1.0000
K2PIG	14	-1701.9389	3431.8778	0.0000	1.0000
SYMG	17	-1700.5249	3435.0497	0.0000	1.0000
SYMI	17	-1700.6418	3435.2836	0.0000	1.0000
SYMIG	18	-1700.5130	3437.0261	0.0000	1.0000
K2P	12	-1712.2081	3448.4163	0.0000	1.0000
SYM	16	-1710.2196	3452.4392	0.0000	1.0000
JC69G	12	-1729.4123	3482.8246	0.0000	1.0000
JC69I	12	-1729.4526	3482.9053	0.0000	1.0000
JC69IG	13	-1729.4174	3484.8347	0.0000	1.0000
JC69	11	-1739.6150	3501.2300	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1638.5891	3381.0298	0.5426	0.5426
HKYI	16	-1638.8126	3381.4768	0.4339	0.9765
HKYIG	17	-1638.5866	3387.5156	0.0212	0.9977
HKY	15	-1647.7802	3392.9213	0.0014	0.9991
GTRG	20	-1632.5370	3394.8884	0.0005	0.9996
GTRI	20	-1632.9017	3395.6178	0.0004	1.0000
GTRIG	21	-1632.5279	3401.3611	0.0000	1.0000
GTR	19	-1641.2039	3405.7316	0.0000	1.0000
F81G	15	-1674.4338	3446.2285	0.0000	1.0000
F81I	15	-1674.5627	3446.4863	0.0000	1.0000
F81IG	16	-1674.4374	3452.7265	0.0000	1.0000
F81	14	-1683.8233	3458.5167	0.0000	1.0000
K2PG	13	-1701.9422	3488.2638	0.0000	1.0000
K2PI	13	-1701.9725	3488.3244	0.0000	1.0000
K2PIG	14	-1701.9389	3494.7479	0.0000	1.0000
K2P	12	-1712.2081	3502.3050	0.0000	1.0000
SYMG	17	-1700.5249	3511.3920	0.0000	1.0000
SYMI	17	-1700.6418	3511.6259	0.0000	1.0000
SYMIG	18	-1700.5130	3517.8591	0.0000	1.0000
SYM	16	-1710.2196	3524.2907	0.0000	1.0000
JC69G	12	-1729.4123	3536.7133	0.0000	1.0000
JC69I	12	-1729.4526	3536.7940	0.0000	1.0000
JC69IG	13	-1729.4174	3543.2141	0.0000	1.0000
JC69	11	-1739.6150	3550.6279	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 591-tomodel-phylip/chr13_5324.phylip = 1 - 659;
 Partition Dummy = 1:591-tomodel-phylip/chr13_5324.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 591-tomodel-phylip/chr13_5324.phylip = 1 - 659;
 Partition Dummy = 1:591-tomodel-phylip/chr13_5324.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_5324.phylip = 1 - 659;
 Partition Dummy = 1:591-tomodel-phylip/chr13_5324.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

