This is an example of implementation of the damage function proposed by Estrada F, Tol RSJ and Botzen WJW in EXTENDING INTEGRATED ASSESSMENT MODELS' DAMAGE FUNCTIONS TO INCLUDE ADAPTATION AND DYNAMIC
We present a new computational and statistical approach for fitting isotonic models under convex differentiable loss functions through recursive partitioning. Models along the partitioning path are al