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

Input data from file "683-tomodel-phylip/chr4_6701.phylip" (7 taxa, 667 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	16	-1899.3671	3831.5711	0.3686	0.3686
HKYI	16	-1899.8438	3832.5245	0.2288	0.5974
HKYIG	17	-1899.3021	3833.5472	0.1372	0.7346
GTRG	20	-1896.1654	3833.6310	0.1316	0.8662
GTRI	20	-1896.5622	3834.4246	0.0885	0.9547
GTRIG	21	-1896.1661	3835.7647	0.0453	1.0000
K2PG	13	-1915.9213	3858.3999	0.0000	1.0000
K2PI	13	-1916.3078	3859.1731	0.0000	1.0000
K2PIG	14	-1915.8913	3860.4268	0.0000	1.0000
SYMG	17	-1913.4791	3861.9012	0.0000	1.0000
SYMI	17	-1913.7565	3862.4561	0.0000	1.0000
SYMIG	18	-1913.4902	3864.0359	0.0000	1.0000
HKY	15	-1930.8560	3892.4493	0.0000	1.0000
F81G	15	-1931.5904	3893.9181	0.0000	1.0000
F81I	15	-1932.0627	3894.8627	0.0000	1.0000
GTR	19	-1928.1366	3895.4479	0.0000	1.0000
F81IG	16	-1931.5395	3895.9159	0.0000	1.0000
JC69G	12	-1945.2917	3915.0605	0.0000	1.0000
JC69I	12	-1945.7006	3915.8783	0.0000	1.0000
K2P	12	-1945.8194	3916.1160	0.0000	1.0000
JC69IG	13	-1945.2661	3917.0897	0.0000	1.0000
SYM	16	-1943.0598	3918.9566	0.0000	1.0000
F81	14	-1962.0843	3952.8127	0.0000	1.0000
JC69	11	-1974.9795	3972.3620	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1899.3671	3830.7342	0.3411	0.3411
HKYI	16	-1899.8438	3831.6875	0.2118	0.5529
GTRG	20	-1896.1654	3832.3307	0.1535	0.7064
HKYIG	17	-1899.3021	3832.6043	0.1339	0.8403
GTRI	20	-1896.5622	3833.1243	0.1032	0.9436
GTRIG	21	-1896.1661	3834.3321	0.0564	1.0000
K2PG	13	-1915.9213	3857.8425	0.0000	1.0000
K2PI	13	-1916.3078	3858.6156	0.0000	1.0000
K2PIG	14	-1915.8913	3859.7826	0.0000	1.0000
SYMG	17	-1913.4791	3860.9582	0.0000	1.0000
SYMI	17	-1913.7565	3861.5131	0.0000	1.0000
SYMIG	18	-1913.4902	3862.9803	0.0000	1.0000
HKY	15	-1930.8560	3891.7120	0.0000	1.0000
F81G	15	-1931.5904	3893.1808	0.0000	1.0000
F81I	15	-1932.0627	3894.1254	0.0000	1.0000
GTR	19	-1928.1366	3894.2733	0.0000	1.0000
F81IG	16	-1931.5395	3895.0790	0.0000	1.0000
JC69G	12	-1945.2917	3914.5834	0.0000	1.0000
JC69I	12	-1945.7006	3915.4012	0.0000	1.0000
K2P	12	-1945.8194	3915.6389	0.0000	1.0000
JC69IG	13	-1945.2661	3916.5323	0.0000	1.0000
SYM	16	-1943.0598	3918.1197	0.0000	1.0000
F81	14	-1962.0843	3952.1685	0.0000	1.0000
JC69	11	-1974.9795	3971.9590	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1899.3671	3902.7788	0.6009	0.6009
HKYI	16	-1899.8438	3903.7322	0.3731	0.9740
HKYIG	17	-1899.3021	3909.1517	0.0248	0.9988
K2PG	13	-1915.9213	3916.3788	0.0007	0.9995
K2PI	13	-1916.3078	3917.1519	0.0005	0.9999
GTRG	20	-1896.1654	3922.3865	0.0000	0.9999
K2PIG	14	-1915.8913	3922.8217	0.0000	1.0000
GTRI	20	-1896.5622	3923.1801	0.0000	1.0000
GTRIG	21	-1896.1661	3928.8907	0.0000	1.0000
SYMG	17	-1913.4791	3937.5057	0.0000	1.0000
SYMI	17	-1913.7565	3938.0605	0.0000	1.0000
SYMIG	18	-1913.4902	3944.0305	0.0000	1.0000
HKY	15	-1930.8560	3959.2539	0.0000	1.0000
F81G	15	-1931.5904	3960.7227	0.0000	1.0000
F81I	15	-1932.0627	3961.6672	0.0000	1.0000
F81IG	16	-1931.5395	3967.1237	0.0000	1.0000
JC69G	12	-1945.2917	3968.6169	0.0000	1.0000
JC69I	12	-1945.7006	3969.4347	0.0000	1.0000
K2P	12	-1945.8194	3969.6724	0.0000	1.0000
JC69IG	13	-1945.2661	3975.0686	0.0000	1.0000
GTR	19	-1928.1366	3979.8263	0.0000	1.0000
SYM	16	-1943.0598	3990.1643	0.0000	1.0000
F81	14	-1962.0843	4015.2076	0.0000	1.0000
JC69	11	-1974.9795	4021.4897	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 683-tomodel-phylip/chr4_6701.phylip = 1 - 667;
 Partition Dummy = 1:683-tomodel-phylip/chr4_6701.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 683-tomodel-phylip/chr4_6701.phylip = 1 - 667;
 Partition Dummy = 1:683-tomodel-phylip/chr4_6701.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 683-tomodel-phylip/chr4_6701.phylip = 1 - 667;
 Partition Dummy = 1:683-tomodel-phylip/chr4_6701.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

