Create measure_abundance.R