Description Usage Arguments Details Value Author(s) References

Finds the maximum likelihood estimate of an identity-link Poisson GLM using an EM algorithm, where each of the coefficients is restricted to be non-negative.

1 2 3 |

`y` |
non-negative integer response vector. |

`x` |
non-negative covariate matrix. |

`standard` |
standardising vector, where each element is a positive constant that (multiplicatively) standardises the fitted value of the corresponding element of the response vector. The default is a vector of ones. |

`offset` |
non-negative additive offset vector. The default is a vector of zeros. |

`start` |
starting values for the parameter estimates. Each element must be
greater than |

`control` |
an |

`accelerate` |
a character string that determines the acceleration
algorithm to be used, (partially) matching one of |

`control.method` |
a list of control parameters for the acceleration algorithm. See |

This is a workhorse function for `addreg`

, and runs the EM algorithm to find the
constrained non-negative MLE associated with an identity-link Poisson GLM. See Marschner (2010)
for full details.

A list containing the following components

`coefficients` |
the constrained non-negative maximum likelihood estimate of the parameters. |

`residuals` |
the residuals at the MLE, that is |

`fitted.values` |
the fitted mean values. |

`rank` |
the number of parameters in the model (named “ |

`family` |
included for compatibility — will always be |

`linear.predictors` |
included for compatibility — same as |

`deviance` |
up to a constant, minus twice the maximised log-likelihood. |

`aic` |
a version of Akaike's |

`aic.c` |
a small-sample corrected
version of Akaike's |

`null.deviance` |
the deviance for the null model, comparable with |

`iter` |
the number of iterations of the EM algorithm used. |

`weights` |
included for compatibility — a vector of ones. |

`prior.weights` |
included for compatibility — a vector of ones. |

`standard` |
the |

`df.residual` |
the residual degrees of freedom. |

`df.null` |
the residual degrees of freedom for the null model. |

`y` |
the |

`converged` |
logical. Did the EM algorithm converge
(according to |

`boundary` |
logical. Is the MLE on the boundary of the parameter
space — i.e. are any of the |

`loglik` |
the maximised log-likelihood. |

`nn.design` |
the non-negative |

Mark W. Donoghoe markdonoghoe@gmail.com.

This function is based on code from Marschner, Gillett and O'Connell (2012) written by Alexandra Gillett.

Hurvich, C. M., J. S. Simonoff and C.-L. Tsai (1998). Smoothing parameter
selection in nonparametric regression using an improved Akaike
information criterion.
*Journal of the Royal Statistical Society: Series B (Statistical Methodology)* 60(2): 271–293.

Marschner, I. C. (2010). Stable computation of maximum likelihood estimates
in identity link Poisson regression. *Journal of Computational and
Graphical Statistics* 19(3): 666–683.

Marschner, I. C., A. C. Gillett and R. L. O'Connell (2012).
Stratified additive Poisson models: Computational methods
and applications in clinical epidemiology. *Computational
Statistics and Data Analysis* 56(5): 1115–1130.

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.