# Looking for the interval for the clock model
setwd("C:/Users/mrnvd/Desktop/Analysis/clock_models/data")
age_data = read.table("atrypids_intervals_FA.txt", header =TRUE)
par(new=TRUE)
hist(age_data$max_age, xlim=c(0, 15))
hist(age_data$min_age, xlim=c(0, 15))