Create fit_dwd.R