Calculates the negative of the Hessian of the log of the observed data likelihood, obtained via Louis's method, evaluated at the last step of the EM algorithm described in Wang et al. (2014+). This
is a support function for PH.ICsurv.EM.
PH.Louis.ICsurv(b, g, bLi, bRi, d1, d2, d3, Xp)estimates of the regression coefficients obtained at the convergence of the EM algorithm.
estimates of the spline coefficients obtained from at the convergence of the EM algorithm.
an I-spline basis matrix of dimension c(length(knots)+order-2, length(x)), corresponding to the left end points of the observed intervals.
an I-spline basis matrix of dimension c(length(knots)+order-2, length(x),), corresponding to the left end points of the observed intervals.
vector indicating whether an observation is left-censored (1) or not (0).
vector indicating whether an observation is interval-censored (1) or not (0).
vector indicating whether an observation is right-censored (1) or not (0).
design matrix of predictor variables (in columns), should be specified without an intercept term.
Hessian matrix.
To obtain the Hessian matrix of the observed likelihood evaluated at the last step of the EM algorithm.
Louis, T. (1982). Finding the observed information matrix when using the EM algorithm. Journal of the Royal Statistical Society, Series B 44, 226-233.
Wang, L., McMahan, C., and Hudgens, M. (2014+). A flexible and computationally efficient method for fitting the proportional hazards model to interval censored data. Submitted.