Create optim.likRatePhylo.R