View source: R/fast.PH.Louis.ICsurv.R

fast.PH.Louis.ICsurv | R Documentation |

This function is a fast version of `PH.Louis.ICsurv`

. It 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 `fast.PH.ICsurv.EM`

.

fast.PH.Louis.ICsurv(b, g, bLi, bRi, d1, d2, d3, Xp)

`b` |
estimates of the regression coefficients obtained at the convergence of the EM algorithm. |

`g` |
estimates of the spline coefficients obtained from at the convergence of the EM algorithm. |

`bLi` |
an I-spline basis matrix of dimension c(length(knots)+order-2, length(x)), corresponding to the left end points of the observed intervals. |

`bRi` |
an I-spline basis matrix of dimension c(length(knots)+order-2, length(x),), corresponding to the left end points of the observed intervals. |

`d1` |
vector indicating whether an observation is left-censored (1) or not (0). |

`d2` |
vector indicating whether an observation is interval-censored (1) or not (0). |

`d3` |
vector indicating whether an observation is right-censored (1) or not (0). |

`Xp` |
design matrix of predictor variables (in columns), should be specified without an intercept term. |

To obtain the Hessian matrix of the observed likelihood evaluated at the last step of the EM algorithm.

Hessian matrix.

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.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.