frailtypack (version 2.2-13)

print.frailtyPenal: Print a Short Summary of parameter estimates of a shared gamma frailty model

Description

Prints a short summary of parameter estimates of a 'frailtyPenal' object

Usage

## S3 method for class 'frailtyPenal':
print(x, digits = max(options()$digits - 4, 3), ...)

Arguments

x

the result of a call to the frailtyPenal function

digits

number of digits to print

...

other unusued arguments

Value

• n
  the number of observations used in the fit.
• n.groups
  the maximum number of groups used in the fit
• n.events
  the number of events observed in the fit
• theta
  variance of frailty parameter
• coef
  the coefficients of the linear predictor, which multiply the columns of the model matrix.
• SE(H)
  the standard error of the estimates deduced from the variance matrix of theta and of the coefficients.
• SE(HIH)
  the standard error of the estimates deduced from the robust estimation of the variance matrix of theta and of the coefficients.
• p
  p-value

See Also

summary.frailtyPenal, frailtyPenal for Shared frailty models, plot.frailtyPenal

Examples

data(readmission)

# /*** Shared frailty model ***/

modSha<-frailtyPenal(Surv(time,event)~as.factor(dukes)+cluster(id),
             n.knots=10,kappa1=10000,data=readmission,Frailty=TRUE)

print(modSha,type="surv",conf=FALSE)

# /*** Cox proportional hazard model using Penalized likelihood ***/

modCox<-frailtyPenal(Surv(time,event)~as.factor(dukes)+cluster(id),
             n.knots=10,kappa1=10000,data=readmission,Frailty=FALSE)
print(modCox)