
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:35 02/17/2011

Input data from file "183-phylip/chr2_22071.phylip" (29 taxa, 488 characters)

Minimum AIC  model: GTRIG
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	60	-1510.4906	3158.1242	0.4957	0.4957
HKYIG	61	-1510.0313	3159.8185	0.2125	0.7082
GTRIG	65	-1505.3706	3161.0730	0.1135	0.8217
GTRG	64	-1506.9059	3161.4808	0.0925	0.9143
HKYI	60	-1512.7483	3162.6394	0.0519	0.9661
GTRI	64	-1507.9746	3163.6183	0.0318	0.9979
K2PIG	58	-1519.2141	3170.3816	0.0011	0.9990
K2PG	57	-1521.0023	3171.3813	0.0007	0.9996
K2PI	57	-1521.8717	3173.1202	0.0003	0.9999
SYMIG	62	-1517.0908	3176.5627	0.0000	1.0000
SYMG	61	-1518.9053	3177.5665	0.0000	1.0000
SYMI	61	-1519.7803	3179.3165	0.0000	1.0000
F81G	59	-1544.9629	3224.4678	0.0000	1.0000
F81IG	60	-1544.5471	3226.2370	0.0000	1.0000
F81I	59	-1547.2059	3228.9539	0.0000	1.0000
HKY	59	-1548.9088	3232.3596	0.0000	1.0000
JC69IG	57	-1552.1382	3233.6531	0.0000	1.0000
JC69G	56	-1553.8624	3234.5369	0.0000	1.0000
GTR	63	-1544.8111	3234.6411	0.0000	1.0000
JC69I	56	-1556.4758	3239.7636	0.0000	1.0000
K2P	56	-1560.2223	3247.2567	0.0000	1.0000
SYM	60	-1557.2764	3251.6956	0.0000	1.0000
F81	58	-1583.0768	3298.1071	0.0000	1.0000
JC69	55	-1592.4412	3309.1417	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1505.3706	3140.7413	0.3044	0.3044
HKYG	60	-1510.4906	3140.9813	0.2700	0.5744
GTRG	64	-1506.9059	3141.8118	0.1782	0.7527
HKYIG	61	-1510.0313	3142.0627	0.1572	0.9099
GTRI	64	-1507.9746	3143.9493	0.0612	0.9711
HKYI	60	-1512.7483	3145.4965	0.0282	0.9994
K2PIG	58	-1519.2141	3154.4282	0.0003	0.9997
K2PG	57	-1521.0023	3156.0046	0.0001	0.9999
K2PI	57	-1521.8717	3157.7435	0.0001	0.9999
SYMIG	62	-1517.0908	3158.1815	0.0000	1.0000
SYMG	61	-1518.9053	3159.8106	0.0000	1.0000
SYMI	61	-1519.7803	3161.5606	0.0000	1.0000
F81G	59	-1544.9629	3207.9258	0.0000	1.0000
F81IG	60	-1544.5471	3209.0941	0.0000	1.0000
F81I	59	-1547.2059	3212.4118	0.0000	1.0000
GTR	63	-1544.8111	3215.6222	0.0000	1.0000
HKY	59	-1548.9088	3215.8176	0.0000	1.0000
JC69IG	57	-1552.1382	3218.2763	0.0000	1.0000
JC69G	56	-1553.8624	3219.7249	0.0000	1.0000
JC69I	56	-1556.4758	3224.9515	0.0000	1.0000
K2P	56	-1560.2223	3232.4447	0.0000	1.0000
SYM	60	-1557.2764	3234.5527	0.0000	1.0000
F81	58	-1583.0768	3282.1537	0.0000	1.0000
JC69	55	-1592.4412	3294.8824	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1510.4906	3392.4002	0.5980	0.5980
K2PG	57	-1521.0023	3394.8525	0.1755	0.7734
K2PI	57	-1521.8717	3396.5915	0.0735	0.8470
HKYI	60	-1512.7483	3396.9154	0.0625	0.9095
K2PIG	58	-1519.2141	3397.4665	0.0475	0.9570
HKYIG	61	-1510.0313	3397.6719	0.0428	0.9998
GTRG	64	-1506.9059	3409.9920	0.0001	0.9999
GTRI	64	-1507.9746	3412.1295	0.0000	1.0000
GTRIG	65	-1505.3706	3413.1118	0.0000	1.0000
SYMG	61	-1518.9053	3415.4199	0.0000	1.0000
SYMI	61	-1519.7803	3417.1699	0.0000	1.0000
SYMIG	62	-1517.0908	3417.9811	0.0000	1.0000
JC69G	56	-1553.8624	3454.3825	0.0000	1.0000
F81G	59	-1544.9629	3455.1544	0.0000	1.0000
JC69IG	57	-1552.1382	3457.1243	0.0000	1.0000
JC69I	56	-1556.4758	3459.6092	0.0000	1.0000
F81I	59	-1547.2059	3459.6404	0.0000	1.0000
F81IG	60	-1544.5471	3460.5131	0.0000	1.0000
HKY	59	-1548.9088	3463.0462	0.0000	1.0000
K2P	56	-1560.2223	3467.1023	0.0000	1.0000
GTR	63	-1544.8111	3479.6121	0.0000	1.0000
SYM	60	-1557.2764	3485.9717	0.0000	1.0000
F81	58	-1583.0768	3525.1920	0.0000	1.0000
JC69	55	-1592.4412	3525.3497	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 183-phylip/chr2_22071.phylip = 1 - 488;
 Partition Dummy = 1:183-phylip/chr2_22071.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_22071.phylip = 1 - 488;
 Partition Dummy = 1:183-phylip/chr2_22071.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 183-phylip/chr2_22071.phylip = 1 - 488;
 Partition Dummy = 1:183-phylip/chr2_22071.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

