predict.singleRStaticCountData function

Predict method for singleRStaticCountData class

A method for predict function, works analogous to predict.glm

but gives the possibility to get standard errors of mean/distribution parameters and directly get pop size estimates for new data.

## S3 method for class 'singleRStaticCountData'
predict(
  object,
  newdata,
  type = c("response", "link", "mean", "popSize", "contr"),
  se.fit = FALSE,
  na.action = NULL,
  weights,
  cov,
  ...
)

Arguments

• object: an object of singleRStaticCountData class.

• newdata: an optional data.frame containing new data.

• type: the type of prediction required, possible values are:

  • "response"-- For matrix containing estimated distributions parameters.
  • "link" -- For matrix of linear predictors.
  • "mean" -- For fitted values of both \mjseqn Y and \mjseqn Y|Y>0.
  • "contr" -- For inverse probability weights (here named for observation contribution to population size estimate).
  • "popSize" -- For population size estimation. Note this results in a call to redoPopEstimation and it is usually better to call this function directly.

  by default set to "response".

• se.fit: a logical value indicating whether standard errors should be computed. Only matters for type in "response", "mean", "link".

• na.action: does nothing yet.

• weights: optional vector of weights for type in "contr", "popSize".

• cov: optional matrix or function or character specifying either a covariance matrix or a function to compute that covariance matrix. By default vcov.singleRStaticCountData can be set to e.g. vcovHC.

• ...: arguments passed to other functions, for now this only affects vcov.singleRStaticCountData method and cov function.

Returns

Depending on type argument if one of "response", "link", "mean"

a matrix with fitted values and possibly standard errors if se.fit

argument was set to TRUE, if type was set to "contr"

a vector with inverses of probabilities, finally for "popSize"

an object of class popSizeEstResults with its own methods containing population size estimation results.

Details

\loadmathjax

Standard errors are computed with assumption of regression coefficients being asymptotically normally distributed, if this assumption holds then each of linear predictors i.e. each row of \mjseqn \boldsymbol \eta=\boldsymbol X_vlm\boldsymbol \beta

is asymptotically normally distributed and their variances are expressed by well known formula. The mean \mjseqn \mu and distribution parameters are then differentiable functions of asymptotically normally distributed variables and therefore their variances can be computed using (multivariate) delta method.

See Also

redoPopEstimation() stats::summary.glm() estimatePopsize()