estRegSmooth: Discrete Survival Additive Regression Model Fitting

In discSurv: Discrete Time Survival Analysis

Discrete Survival Additive Regression Model Fitting

Description

Wrapper for estimation of discrete survival generalized additive models. Several preprocessing options such as time-dependent covariates are available.

Usage

estRegSmooth(
  dataShort,
  dataTransform = "dataLong",
  formulaVariable = ~s(timeInt),
  timeColumn,
  eventColumn,
  idColumn = NULL,
  timeAsFactor = FALSE,
  family = stats::binomial,
  storeAugData = TRUE,
  ...
)

Arguments

dataShort	Original data set in short format with each row corresponding to one independent observation (class "data.frame").
dataTransform	Specification of the data transformation function from short to long format (class "character"). There are two available options: Without time dependent covariates ("dataLong") and with time dependent covariates ("dataLongTimeDep"). The default is set to the former.
formulaVariable	Specifies the right hand side of the regression formula (class "formula"). The default is to use the discrete time variable, which corresponds to a covariate free hazard. It is recommended to always include the discrete time variable "timeInt".
timeColumn	Character giving the column name of the observed times. It is required that the observed times are discrete (class "integer").
eventColumn	Column name of the event indicator (class "character"). It is required that this is a binary variable with 1=="event" and 0=="censored".
idColumn	Name of column of identification number of persons (class "character"). Default is set to use function dataLong, that does not need this argument.
timeAsFactor	Should the time intervals be coded as factor (class "logical")? Default is FALSE. In the default settings the column is treated as quantitative variable (class "numeric").
family	Specifies the assumption about the response distribution and the link function. Default value is the discrete survival continuation ratio model with logit-link (for comparison see family).
storeAugData	Should the augmented data set be saved (class "logical")? Defaults is TRUE. The data set is available as attribute "augData".
...	Specification of additional arguments in function gam.

Details

Variables in argument formulaVariable need to be separated by "+ ". For example if the two variables timeInt and X1 should be included as smooth functions the formula would be "~ s(timeInt) + s(X1)". The variable timeInt is constructed before estimation of the model.

Value

Returns an object of class "gam".

Note

Splines can also be used to include heterogeneity of mixed models with penalized fixed effects. An advantage of this approach is that one does not have to assume a specific distribution for the random effects.

Model is rescaled to exclude the intercept.

Author(s)

Thomas Welchowski [e-mail]

References

\insertRef

bergerTutorialdiscSurv

\insertReftutzModelDiscdiscSurv

See Also

dataLong, dataLongTimeDep, gam, estRegFrailty

Examples

# Example with unemployment data
library(Ecdat)
data(UnempDur)

# Select subsample
SubUnempDur <- UnempDur[1:100, ]

# Estimate discrete survival continuation ratio model
estRegSmoothModel <- estRegSmooth(dataShort = SubUnempDur, dataTransform = "dataLong", 
formulaVariable =~ s(timeInt) + s(age) + ui + logwage * ui, 
eventColumn = "censor1", timeColumn = "spell")
summary(estRegSmoothModel)

# Estimated smooth functions
plot(estRegSmoothModel)

discSurv documentation built on April 29, 2026, 9:07 a.m.