Description Usage Arguments Value Note Author(s) References Examples

View source: R/find.shapeMLE.R

Compute the maximum likelihood estimator (MLE) of a shape-constrained hazard function under IID sampling. We assume that the data are continuous and allow for right censoring. The function 'find.shapeMLE' allows for four different shape constraints: increasing, decreasing, unimodal, and u-shaped.

1 |

`x` |
vector of length n containing the data |

`delta` |
logical vector containing the (right) censoring information. If delta_i=1 then the observation was not censored. The default is delta_i=1 for all i, that is, assuming that no observations were censored. |

`type` |
string indicating type of shape constraint. Options are "increasing", "decreasing", "unimodal", and "ushaped". |

`plot` |
logical, if TRUE, the graphical representation of the MLE is shown |

A list containing the following elements:

`h.range` |
endpoints for the values of the hazard MLE |

`h.val` |
values of the hazard MLE in between the endpoints |

`phi` |
the criterion function |

`H` |
the cumulative hazard MLE evaluated at the data points |

`mode` |
location of the mode (for unimodal) or antimode (for u-shaped). Note that the antimode is not unique, and the midpoint of all possible values is returned. |

`type` |
string indicating type of shape constraint used |

The MLE can be found in different ways. We use the graphical representation via convex minorants or concave majorants of appropriate functions. Also, for the increasing and unimodal setting, the MLE takes the value of infinity at the mode (the largest observation for `type="increasing"`

). In such situations, this value is removed from the likelihood in the maximization process. A similar approach was taken in Gernander (1956).

Rihong Hui and Hanna Jankowski <[email protected]>

Hui, R. and Jankowski, H. (2012). Maximum likelihood estimation of a shape-constrained hazard in the proportional hazard model. Technical Report. http://www.math.yorku.ca/~hkj/

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 |

CPHshape documentation built on May 30, 2017, 4:32 a.m.

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.