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