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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-1176.5704	2394.8175	0.2424	0.2424
GTRG	20	-1176.5973	2394.8713	0.2360	0.4784
HKYI	16	-1181.1265	2395.3303	0.1876	0.6660
HKYG	16	-1181.1294	2395.3360	0.1871	0.8531
GTRIG	21	-1176.5704	2396.9887	0.0819	0.9350
HKYIG	17	-1181.1175	2397.4493	0.0650	1.0000
HKY	15	-1199.1933	2429.3353	0.0000	1.0000
GTR	19	-1195.1851	2429.8842	0.0000	1.0000
K2PI	13	-1206.3499	2439.4164	0.0000	1.0000
K2PG	13	-1206.3527	2439.4218	0.0000	1.0000
K2PIG	14	-1206.3527	2441.5338	0.0000	1.0000
SYMI	17	-1203.8389	2442.8920	0.0000	1.0000
SYMG	17	-1203.8400	2442.8943	0.0000	1.0000
SYMIG	18	-1203.8400	2445.0399	0.0000	1.0000
F81G	15	-1209.2459	2449.4404	0.0000	1.0000
F81I	15	-1209.2825	2449.5136	0.0000	1.0000
F81IG	16	-1209.0225	2451.1222	0.0000	1.0000
K2P	12	-1223.2239	2471.0607	0.0000	1.0000
SYM	16	-1220.2922	2473.6616	0.0000	1.0000
F81	14	-1223.5015	2475.8315	0.0000	1.0000
JC69G	12	-1230.6167	2485.8463	0.0000	1.0000
JC69I	12	-1230.6378	2485.8886	0.0000	1.0000
JC69IG	13	-1230.4770	2487.6706	0.0000	1.0000
JC69	11	-1244.5790	2511.6756	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1176.5704	2393.1408	0.2703	0.2703
GTRG	20	-1176.5973	2393.1946	0.2631	0.5334
HKYI	16	-1181.1265	2394.2530	0.1550	0.6885
HKYG	16	-1181.1294	2394.2588	0.1546	0.8430
GTRIG	21	-1176.5704	2395.1407	0.0994	0.9425
HKYIG	17	-1181.1175	2396.2350	0.0575	1.0000
GTR	19	-1195.1851	2428.3703	0.0000	1.0000
HKY	15	-1199.1933	2428.3867	0.0000	1.0000
K2PI	13	-1206.3499	2438.6998	0.0000	1.0000
K2PG	13	-1206.3527	2438.7053	0.0000	1.0000
K2PIG	14	-1206.3527	2440.7054	0.0000	1.0000
SYMI	17	-1203.8389	2441.6777	0.0000	1.0000
SYMG	17	-1203.8400	2441.6800	0.0000	1.0000
SYMIG	18	-1203.8400	2443.6801	0.0000	1.0000
F81G	15	-1209.2459	2448.4918	0.0000	1.0000
F81I	15	-1209.2825	2448.5650	0.0000	1.0000
F81IG	16	-1209.0225	2450.0450	0.0000	1.0000
K2P	12	-1223.2239	2470.4477	0.0000	1.0000
SYM	16	-1220.2922	2472.5844	0.0000	1.0000
F81	14	-1223.5015	2475.0031	0.0000	1.0000
JC69G	12	-1230.6167	2485.2334	0.0000	1.0000
JC69I	12	-1230.6378	2485.2757	0.0000	1.0000
JC69IG	13	-1230.4770	2486.9541	0.0000	1.0000
JC69	11	-1244.5790	2511.1580	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1181.1265	2462.3757	0.4897	0.4897
HKYG	16	-1181.1294	2462.3815	0.4883	0.9780
HKYIG	17	-1181.1175	2468.6154	0.0216	0.9997
GTRI	20	-1176.5704	2478.2942	0.0002	0.9998
GTRG	20	-1176.5973	2478.3480	0.0002	1.0000
GTRIG	21	-1176.5704	2484.5518	0.0000	1.0000
HKY	15	-1199.1933	2492.2517	0.0000	1.0000
K2PI	13	-1206.3499	2494.0495	0.0000	1.0000
K2PG	13	-1206.3527	2494.0550	0.0000	1.0000
K2PIG	14	-1206.3527	2500.3127	0.0000	1.0000
GTR	19	-1195.1851	2509.2660	0.0000	1.0000
F81G	15	-1209.2459	2512.3568	0.0000	1.0000
F81I	15	-1209.2825	2512.4300	0.0000	1.0000
SYMI	17	-1203.8389	2514.0581	0.0000	1.0000
SYMG	17	-1203.8400	2514.0604	0.0000	1.0000
F81IG	16	-1209.0225	2518.1676	0.0000	1.0000
SYMIG	18	-1203.8400	2520.3181	0.0000	1.0000
K2P	12	-1223.2239	2521.5398	0.0000	1.0000
F81	14	-1223.5015	2534.6104	0.0000	1.0000
JC69G	12	-1230.6167	2536.3254	0.0000	1.0000
JC69I	12	-1230.6378	2536.3677	0.0000	1.0000
SYM	16	-1220.2922	2540.7071	0.0000	1.0000
JC69IG	13	-1230.4770	2542.3037	0.0000	1.0000
JC69	11	-1244.5790	2557.9923	0.0000	1.0000



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

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

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


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

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


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

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


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

