Compute Hessian of Unlinked Monotone Regression objective function from Balabdaoui, Doss, and Durot
UMRhess_generic(mm, ww_m, yy, ww_y = rep(1/length(yy), length(yy)), dens, BBp)
Current (unsorted) estimate/iterate at which to compute gradient. (Length is <= than the number of X observations in the problem).
Weights (nonnegative, sum to 1) corresponding to mm. Same length as mm.
Y (response) observation vector (numeric)
Weights (nonnegative, sum to 1) corresponding to yy. Same length as yy. Default is just 1/length(yy) for each value.
This is the error density, a function object (Balabdaoui, Doss, Durot (2020+). Function accepting vector or matrix arguments.
This is derivative of "B" function ("B prime"), where B is defined in the paper. Function accepting vector or matrix arguments.
See paper for derivations.