Description Usage Arguments Value Author(s) References See Also

A Kiefer-Wolfowitz NPMLE procedure for estimation of a Gaussian model with independent mean and variance prior components with weighted longitudinal data. This version iterates back and forth from Gamma and Gaussian forms of the likelihood.

1 |

`y` |
A vector of observations |

`id` |
A strata indicator vector indicating grouping of y |

`w` |
A vector of weights corresponding to y |

`u` |
A vector of bin boundaries for the mean effects |

`v` |
A vector of bin boundaries for the variance effects |

`eps` |
Convergence tolerance for iterations |

`maxit` |
A limit on the number of allowed iterations |

`...` |
optional parameters to be passed to KWDual to control optimization |

A list consisting of the following components:

`u` |
midpoints of the mean bin boundaries |

`fu` |
the function values of the mixing density of the means |

`v` |
midpoints of the variance bin boundaries |

`fv` |
the function values of the mixing density of the variances. |

`logLik` |
vector of log likelihood values for each iteration |

`du` |
Bayes rule estimate of the mixing density means. |

`dv` |
Bayes rule estimate of the mixing density variances. |

`status` |
Mosek convergence status for each iteration |

J. Gu and R. Koenker

Gu, J. and R. Koenker (2015) Empirical Bayesball Remixed: Empirical
Bayes Methods for Longitudinal Data, *J. Applied Econometrics*, 32, 575-599.

Koenker, R. and J. Gu, (2017) REBayes: An R Package for Empirical Bayes Mixture Methods,
*Journal of Statistical Software*, 82, 1–26.

WGLVmix for a more general bivariate mixing distribution version and WTLVmix for an alternative estimator exploiting a Student/Gamma decomposition

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.