### ### for (z in 1:6){ data <- readDiscreteCharacterData(datasets[z]) source("simple_heterogeneity_script.Rev") clear(data) } q()