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

Input data from file "591-tomodel-phylip/chr1_7424.phylip" (7 taxa, 390 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	-1158.9525	2351.3634	0.4608	0.4608
HKYI	16	-1159.6792	2352.8168	0.2228	0.6836
HKYIG	17	-1158.9020	2353.4492	0.1624	0.8461
GTRG	20	-1156.2786	2354.8336	0.0813	0.9273
GTRI	20	-1156.8755	2356.0274	0.0447	0.9721
GTRIG	21	-1156.2300	2356.9708	0.0279	1.0000
SYMG	17	-1175.9734	2387.5919	0.0000	1.0000
HKY	15	-1178.6765	2388.6365	0.0000	1.0000
GTR	19	-1174.4138	2388.8818	0.0000	1.0000
SYMIG	18	-1175.9080	2389.6596	0.0000	1.0000
SYMI	17	-1177.3899	2390.4250	0.0000	1.0000
K2PG	13	-1183.3680	2393.7041	0.0000	1.0000
K2PIG	14	-1183.2909	2395.7018	0.0000	1.0000
K2PI	13	-1184.4861	2395.9403	0.0000	1.0000
F81G	15	-1185.7465	2402.7764	0.0000	1.0000
F81I	15	-1186.5134	2404.3103	0.0000	1.0000
F81IG	16	-1185.6758	2404.8100	0.0000	1.0000
SYM	16	-1191.7053	2416.8690	0.0000	1.0000
K2P	12	-1200.1141	2425.0558	0.0000	1.0000
F81	14	-1202.4296	2433.9791	0.0000	1.0000
JC69G	12	-1205.6962	2436.2201	0.0000	1.0000
JC69I	12	-1206.6334	2438.0944	0.0000	1.0000
JC69IG	13	-1205.6358	2438.2396	0.0000	1.0000
JC69	11	-1220.8289	2464.3561	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1158.9525	2349.9049	0.4194	0.4194
HKYI	16	-1159.6792	2351.3584	0.2028	0.6221
HKYIG	17	-1158.9020	2351.8041	0.1623	0.7844
GTRG	20	-1156.2786	2352.5572	0.1113	0.8957
GTRI	20	-1156.8755	2353.7510	0.0613	0.9570
GTRIG	21	-1156.2300	2354.4599	0.0430	1.0000
SYMG	17	-1175.9734	2385.9467	0.0000	1.0000
GTR	19	-1174.4138	2386.8277	0.0000	1.0000
HKY	15	-1178.6765	2387.3531	0.0000	1.0000
SYMIG	18	-1175.9080	2387.8160	0.0000	1.0000
SYMI	17	-1177.3899	2388.7798	0.0000	1.0000
K2PG	13	-1183.3680	2392.7360	0.0000	1.0000
K2PIG	14	-1183.2909	2394.5818	0.0000	1.0000
K2PI	13	-1184.4861	2394.9722	0.0000	1.0000
F81G	15	-1185.7465	2401.4930	0.0000	1.0000
F81I	15	-1186.5134	2403.0269	0.0000	1.0000
F81IG	16	-1185.6758	2403.3515	0.0000	1.0000
SYM	16	-1191.7053	2415.4106	0.0000	1.0000
K2P	12	-1200.1141	2424.2282	0.0000	1.0000
F81	14	-1202.4296	2432.8591	0.0000	1.0000
JC69G	12	-1205.6962	2435.3925	0.0000	1.0000
JC69I	12	-1206.6334	2437.2669	0.0000	1.0000
JC69IG	13	-1205.6358	2437.2716	0.0000	1.0000
JC69	11	-1220.8289	2463.6577	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1158.9525	2413.3633	0.6507	0.6507
HKYI	16	-1159.6792	2414.8167	0.3146	0.9652
HKYIG	17	-1158.9020	2419.2286	0.0347	0.9999
GTRG	20	-1156.2786	2431.8801	0.0001	1.0000
GTRI	20	-1156.8755	2433.0740	0.0000	1.0000
GTRIG	21	-1156.2300	2437.7490	0.0000	1.0000
K2PG	13	-1183.3680	2444.2959	0.0000	1.0000
K2PI	13	-1184.4861	2446.5321	0.0000	1.0000
HKY	15	-1178.6765	2446.8453	0.0000	1.0000
K2PIG	14	-1183.2909	2450.1079	0.0000	1.0000
SYMG	17	-1175.9734	2453.3712	0.0000	1.0000
SYMI	17	-1177.3899	2456.2043	0.0000	1.0000
SYMIG	18	-1175.9080	2459.2066	0.0000	1.0000
F81G	15	-1185.7465	2460.9852	0.0000	1.0000
GTR	19	-1174.4138	2462.1845	0.0000	1.0000
F81I	15	-1186.5134	2462.5191	0.0000	1.0000
F81IG	16	-1185.6758	2466.8099	0.0000	1.0000
K2P	12	-1200.1141	2471.8220	0.0000	1.0000
SYM	16	-1191.7053	2478.8689	0.0000	1.0000
JC69G	12	-1205.6962	2482.9863	0.0000	1.0000
JC69I	12	-1206.6334	2484.8606	0.0000	1.0000
F81	14	-1202.4296	2488.3852	0.0000	1.0000
JC69IG	13	-1205.6358	2488.8315	0.0000	1.0000
JC69	11	-1220.8289	2507.2853	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 591-tomodel-phylip/chr1_7424.phylip = 1 - 390;
 Partition Dummy = 1:591-tomodel-phylip/chr1_7424.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 591-tomodel-phylip/chr1_7424.phylip = 1 - 390;
 Partition Dummy = 1:591-tomodel-phylip/chr1_7424.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 591-tomodel-phylip/chr1_7424.phylip = 1 - 390;
 Partition Dummy = 1:591-tomodel-phylip/chr1_7424.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

