Takes a fitted vagam object produced by the main vagam function and produces predictions given a new set of values for the model covariates or the original values used for the model fit.
# S3 method for vagam
predict(object, new.smoothX, new.paraX = NULL, terms = NULL,
alpha = 0.05, type = "link", ...)An object of class "vagam".
A new matrix of covariates, each of which are were entered as additive smooth terms in the fitted GAM.
A new matrix of covariates, each of which were be entered as parametric terms in the fitted GAM. Note the predictions will account for the intercept ONLY if new.paraX is supplied.
If terms = NULL, prediction is made across all smoothing covariates. Else, prediction is made to the specified smoothing covariate, with no intercept added.
Level of significance for the pointwise confidence bands for predictions.
When type = "link" (default) the linear predictor (with associated standard errors) is returned. When type = "response" predictions on the scale of the response are returned (with associated standard errors).
This is currently ignored.
A data frame containing information such as the predicted response, standard errors, and lower and upper bounds of the pointwise confidence bands.
Current implemented a basic method of constructing predictions either for a single smoothing covariate, or across all the smoothing (and parametric if supplied) covariates, based on a GAM fitted using the main vagam function. By default, standard errors and this pointwise confidence bands are also produced based on fitted vagam object. Under the variational approximations framework, the smooths and confidence bands are constructed based on the variational approximation to the posterior distribution of the smoothing coefficients (which are treated as random effects with a normal prior under the mixed model framework). Please see Hui et al., (2018) for more information.
Hui, F. K. C., You, C., Shang, H. L., and Mueller, S. (2018). Semiparametric regression using variational approximations, Journal of the American Statistical Association, forthcoming.
vagam for the main fitting function
# NOT RUN {
## Please see examples in the help file for the vagam function.
# }
Run the code above in your browser using DataLab