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

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

Minimum AIC  model: SYMI
Minimum AICc model: SYMI
Minimum BIC  model: K2PI


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
SYMI	17	-1574.0411	3183.0152	0.4620	0.4620
SYMIG	18	-1573.5084	3184.0611	0.2739	0.7358
SYMG	17	-1575.6782	3186.2894	0.0899	0.8257
GTRI	20	-1572.7930	3186.8725	0.0671	0.8928
GTRIG	21	-1572.3363	3188.0897	0.0365	0.9294
K2PI	13	-1581.0906	3188.7327	0.0265	0.9559
K2PIG	14	-1580.2860	3189.2094	0.0209	0.9767
GTRG	20	-1574.4888	3190.2640	0.0123	0.9891
K2PG	13	-1582.5354	3191.6224	0.0062	0.9953
HKYI	16	-1580.3753	3193.5786	0.0023	0.9977
HKYIG	17	-1579.5971	3194.1271	0.0018	0.9994
HKYG	16	-1581.8078	3196.4436	0.0006	1.0000
JC69I	12	-1616.5275	3257.5270	0.0000	1.0000
JC69IG	13	-1616.0974	3258.7463	0.0000	1.0000
JC69G	12	-1617.7189	3259.9099	0.0000	1.0000
F81I	15	-1615.6256	3261.9808	0.0000	1.0000
F81IG	16	-1615.2159	3263.2598	0.0000	1.0000
SYM	16	-1615.6409	3264.1098	0.0000	1.0000
F81G	15	-1616.8244	3264.3783	0.0000	1.0000
K2P	12	-1621.2233	3266.9187	0.0000	1.0000
GTR	19	-1614.7592	3268.6806	0.0000	1.0000
HKY	15	-1620.4649	3271.6592	0.0000	1.0000
JC69	11	-1655.5622	3333.5233	0.0000	1.0000
F81	14	-1654.7002	3338.0377	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMI	17	-1574.0411	3182.0823	0.4475	0.4475
SYMIG	18	-1573.5084	3183.0168	0.2805	0.7280
SYMG	17	-1575.6782	3185.3565	0.0871	0.8150
GTRI	20	-1572.7930	3185.5861	0.0776	0.8926
GTRIG	21	-1572.3363	3186.6725	0.0451	0.9377
K2PI	13	-1581.0906	3188.1812	0.0212	0.9589
K2PIG	14	-1580.2860	3188.5721	0.0174	0.9764
GTRG	20	-1574.4888	3188.9776	0.0142	0.9906
K2PG	13	-1582.5354	3191.0709	0.0050	0.9956
HKYI	16	-1580.3753	3192.7506	0.0022	0.9978
HKYIG	17	-1579.5971	3193.1942	0.0017	0.9995
HKYG	16	-1581.8078	3195.6156	0.0005	1.0000
JC69I	12	-1616.5275	3257.0549	0.0000	1.0000
JC69IG	13	-1616.0974	3258.1948	0.0000	1.0000
JC69G	12	-1617.7189	3259.4379	0.0000	1.0000
F81I	15	-1615.6256	3261.2513	0.0000	1.0000
F81IG	16	-1615.2159	3262.4318	0.0000	1.0000
SYM	16	-1615.6409	3263.2818	0.0000	1.0000
F81G	15	-1616.8244	3263.6488	0.0000	1.0000
K2P	12	-1621.2233	3266.4466	0.0000	1.0000
GTR	19	-1614.7592	3267.5185	0.0000	1.0000
HKY	15	-1620.4649	3270.9297	0.0000	1.0000
JC69	11	-1655.5622	3333.1245	0.0000	1.0000
F81	14	-1654.7002	3337.4003	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1581.0906	3246.8532	0.7546	0.7546
K2PG	13	-1582.5354	3249.7429	0.1779	0.9325
K2PIG	14	-1580.2860	3251.7573	0.0650	0.9975
SYMI	17	-1574.0411	3258.8072	0.0019	0.9994
SYMG	17	-1575.6782	3262.0814	0.0004	0.9998
SYMIG	18	-1573.5084	3264.2549	0.0001	0.9999
HKYI	16	-1580.3753	3264.9623	0.0001	1.0000
HKYG	16	-1581.8078	3267.8272	0.0000	1.0000
HKYIG	17	-1579.5971	3269.9191	0.0000	1.0000
GTRI	20	-1572.7930	3275.8507	0.0000	1.0000
GTRG	20	-1574.4888	3279.2422	0.0000	1.0000
GTRIG	21	-1572.3363	3281.4504	0.0000	1.0000
JC69I	12	-1616.5275	3311.2137	0.0000	1.0000
JC69G	12	-1617.7189	3313.5966	0.0000	1.0000
JC69IG	13	-1616.0974	3316.8668	0.0000	1.0000
K2P	12	-1621.2233	3320.6054	0.0000	1.0000
F81I	15	-1615.6256	3328.9497	0.0000	1.0000
F81G	15	-1616.8244	3331.3473	0.0000	1.0000
F81IG	16	-1615.2159	3334.6435	0.0000	1.0000
SYM	16	-1615.6409	3335.4935	0.0000	1.0000
HKY	15	-1620.4649	3338.6282	0.0000	1.0000
GTR	19	-1614.7592	3353.2699	0.0000	1.0000
JC69	11	-1655.5622	3382.7700	0.0000	1.0000
F81	14	-1654.7002	3400.5856	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_23648.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23648.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_23648.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23648.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_23648.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23648.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

