pgam (version 0.3.3)

pgam: Poisson-Gamma Additive Models

Description

Fit Poisson-Gamma Additive Models using the roughness penalty approach

Usage

pgam(formula, dataset, omega = 0.5, beta = 0.1, offset = 1, digits = getOption("digits"),
maxit = 100, eps = 1e-06, control = list(trace = 10, REPORT = 10, maxit = 100, abstol = 1e-04,
reltol = 0.001, factr = 1e+07, pgtol = 0), optim.method = "BFGS", partial.resid = "response",
smoother = "spline", bkf.eps = 0.001, bkf.maxit = 100, numerical.se = TRUE, verbose = TRUE)

Arguments

formula

a model formula. See pgam.parser for details

dataset

a data set in the environment search path

omega

initial value for the discount factor

beta

vector of initial values for covariates coefficients. If a sigle value is supplied it is replicated to fill in the whole vector

offset

default is $1$. Other value can be supplied here

digits

number of decimal places for printing information out

maxit

convergence control iterations

eps

convergence control criterion

control

convergence control of optim. See its help for details

optim.method

optimization method for optim. See its help for details

partial.resid

type of partial residual to be used if semiparametric fitting. See residuals.pgam for details

smoother

smoother to be used in backfitting. See pgam.smooth for details

bkf.eps

convergence control criterion for the backfitting algorithm

bkf.maxit

convergence control iterations for the backfitting algorithm

numerical.se

if TRUE numerical standard error of parameters are returned, else analytical ones are returned

verbose

if TRUE information during estimation process is printed out

Value

• List containing an object of class pgam.

Details

There are a lot of details to be written. It will be soon enough.

Specific information can be obtained on functions help.

This algorithm fits fully parametric Poisson-Gamma model also.

References

Junger, W. L. (2004) Modelo Poisson-Gama Semi-Param�trico: Uma Abordagem de Penaliza��o por Rugosidade. MSc Thesis. Rio de Janeiro, PUC-Rio, Departamento de Engenharia El�trica

Harvey, A. C., Fernandes, C. (1989) Time series models for count data or qualitative observations. Journal of Business and Economic Statistics, 7(4):407--417

Campos, E. L., De Leon, A. C. M. P., Fernandes, C. A. C. (2003) Modelo Poisson-Gama para S�ries Temporais de Dados de Contagem - Teoria e Aplica��es. 10a ESTE - Escola de S�ries Temporais e Econometria

Green, P. J., Silverman, B. W. (1994) Nonparametric Regression and Generalized Linear Models: a roughness penalty approach. Chapman and Hall, London

See Also

predict.pgam, pgam.parser, residuals.pgam, backfitting

Examples

library(pgam)
data(aihrio)
attach(aihrio)
form <- ITRESP5~f(WEEK)+HOLIDAYS+rain+PM+g(tmpmax,7)+g(wet,3)
m <- pgam(form,aihrio,omega=.8,beta=.01,maxit=1e2,eps=1e-4,optim.method="BFGS",partial.resid="response")

summary(m)