### ### for (z in 1:2){ source(datasets[z] + "/root_time.Rev") data <- readDiscreteCharacterData(datasets[z] + "/400concatenated.nexus") start_tree <- readTrees(datasets[z] + "/start.tre")[1] source("simple_clock_script.Rev") source("simple_heterogeneity_script.Rev") clear(data) } q()