Create dushoff_model.R