
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:19 04/04/2011

Input data from file "683-tomodel-phylip/chr9_3587.phylip" (7 taxa, 596 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-1586.1899	3213.8407	0.4575	0.4575
GTRI	20	-1586.4486	3214.3581	0.3532	0.8107
GTRIG	21	-1586.2002	3216.0101	0.1546	0.9653
HKYG	16	-1593.8081	3220.5558	0.0159	0.9812
HKYI	16	-1594.0638	3221.0672	0.0123	0.9936
HKYIG	17	-1593.8167	3222.6922	0.0055	0.9990
SYMG	17	-1596.7773	3228.6135	0.0003	0.9993
SYMI	17	-1597.0309	3229.1206	0.0002	0.9995
K2PG	13	-1601.4765	3229.5785	0.0002	0.9997
K2PI	13	-1601.7635	3230.1524	0.0001	0.9998
SYMIG	18	-1596.7960	3230.7774	0.0001	0.9999
K2PIG	14	-1601.4868	3231.6965	0.0001	1.0000
GTR	19	-1600.1744	3239.6681	0.0000	1.0000
HKY	15	-1608.8746	3248.5767	0.0000	1.0000
SYM	16	-1611.0702	3255.0799	0.0000	1.0000
K2P	12	-1616.3995	3257.3342	0.0000	1.0000
F81G	15	-1618.2917	3267.4109	0.0000	1.0000
F81I	15	-1618.5808	3267.9892	0.0000	1.0000
F81IG	16	-1618.2999	3269.5393	0.0000	1.0000
JC69G	12	-1625.8371	3276.2093	0.0000	1.0000
JC69I	12	-1626.0921	3276.7194	0.0000	1.0000
JC69IG	13	-1625.8506	3278.3265	0.0000	1.0000
F81	14	-1631.2865	3291.2958	0.0000	1.0000
JC69	11	-1639.3576	3301.1672	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1586.1899	3212.3798	0.4555	0.4555
GTRI	20	-1586.4486	3212.8972	0.3517	0.8072
GTRIG	21	-1586.2002	3214.4004	0.1659	0.9731
HKYG	16	-1593.8081	3219.6162	0.0122	0.9853
HKYI	16	-1594.0638	3220.1277	0.0095	0.9948
HKYIG	17	-1593.8167	3221.6334	0.0045	0.9993
SYMG	17	-1596.7773	3227.5547	0.0002	0.9995
SYMI	17	-1597.0309	3228.0618	0.0002	0.9997
K2PG	13	-1601.4765	3228.9530	0.0001	0.9998
K2PI	13	-1601.7635	3229.5269	0.0001	0.9999
SYMIG	18	-1596.7960	3229.5919	0.0001	1.0000
K2PIG	14	-1601.4868	3230.9736	0.0000	1.0000
GTR	19	-1600.1744	3238.3487	0.0000	1.0000
HKY	15	-1608.8746	3247.7491	0.0000	1.0000
SYM	16	-1611.0702	3254.1403	0.0000	1.0000
K2P	12	-1616.3995	3256.7991	0.0000	1.0000
F81G	15	-1618.2917	3266.5834	0.0000	1.0000
F81I	15	-1618.5808	3267.1616	0.0000	1.0000
F81IG	16	-1618.2999	3268.5998	0.0000	1.0000
JC69G	12	-1625.8371	3275.6741	0.0000	1.0000
JC69I	12	-1626.0921	3276.1842	0.0000	1.0000
JC69IG	13	-1625.8506	3277.7011	0.0000	1.0000
F81	14	-1631.2865	3290.5729	0.0000	1.0000
JC69	11	-1639.3576	3300.7152	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1601.4765	3286.0262	0.4854	0.4854
K2PI	13	-1601.7635	3286.6001	0.3643	0.8497
HKYG	16	-1593.8081	3289.8601	0.0714	0.9211
HKYI	16	-1594.0638	3290.3715	0.0553	0.9764
K2PIG	14	-1601.4868	3292.4369	0.0197	0.9961
HKYIG	17	-1593.8167	3296.2675	0.0029	0.9990
GTRG	20	-1586.1899	3300.1846	0.0004	0.9994
GTRI	20	-1586.4486	3300.7021	0.0003	0.9997
SYMG	17	-1596.7773	3302.1888	0.0002	0.9999
SYMI	17	-1597.0309	3302.6959	0.0001	1.0000
GTRIG	21	-1586.2002	3306.5954	0.0000	1.0000
SYMIG	18	-1596.7960	3308.6163	0.0000	1.0000
K2P	12	-1616.3995	3309.4819	0.0000	1.0000
HKY	15	-1608.8746	3313.6028	0.0000	1.0000
GTR	19	-1600.1744	3321.7633	0.0000	1.0000
SYM	16	-1611.0702	3324.3842	0.0000	1.0000
JC69G	12	-1625.8371	3328.3570	0.0000	1.0000
JC69I	12	-1626.0921	3328.8671	0.0000	1.0000
F81G	15	-1618.2917	3332.4370	0.0000	1.0000
F81I	15	-1618.5808	3333.0152	0.0000	1.0000
JC69IG	13	-1625.8506	3334.7742	0.0000	1.0000
F81IG	16	-1618.2999	3338.8436	0.0000	1.0000
JC69	11	-1639.3576	3349.0078	0.0000	1.0000
F81	14	-1631.2865	3352.0363	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 683-tomodel-phylip/chr9_3587.phylip = 1 - 596;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3587.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 683-tomodel-phylip/chr9_3587.phylip = 1 - 596;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3587.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 (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr9_3587.phylip = 1 - 596;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3587.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

