
Solver for the negative Poisson log-likelihood
sSolver(z, a, extra)
Vector containing the fitted values
Vector with Lagrange multipliers
Value of the target function
Gradient at point x
Vector containing observed response
Matrix with active constraints
List with element y
containing the observed response vector
This function is called internally in activeSet
by setting mySolver = sSolver
.
activeSet
##Minimizing Poisson log-liklihood
set.seed(12345)
yp <- rpois(9,5)
Atot <- cbind(1:8, 2:9) ##Matrix defining isotonicity (total order)
x0 <- 1:9 ##starting values
fit.poi <- activeSet(Atot, sSolver, x0 = x0, y = yp)
Run the code above in your browser using DataLab