print.jointmeta1SE: Print function for 'jointmeta1SE' objects

In mesudell/joineRmeta: Joint Modelling for Meta-Analytic (Multi-Study) Data

Description

This function extracts the results of the bootstrapping function jointmetaSE applied to the results of a one stage joint meta model fitted using jointmeta1, namely jointmeta1SE objects. The bootstrapping function returns not only results but also the covariance matrix for the estimated parameters and the results of each bootstrap. This print function allows just the results of the bootstraps to easily be displayed.

Usage

## S3 method for class 'jointmeta1SE'
print(x, ...)

Arguments

x	a jointmeta1SE object, the result of applying jointmetaSE to the joint model fit obtained by applying jointmeta1 to a multi-study joint data set.
...	additional arguments; currently none are used.

Value

a data frame containing the estimates, standard errors and 95% confidence intervals for the parameters from the model and any overall effects requested in the jointmetaSE function call.

See Also

jointmeta1, jointmetaSE, jointmeta1SE.object, jointmeta1.object

Examples

   #change example data to jointdata object
   jointdat2<-tojointdata(longitudinal = simdat2$longitudinal,
   survival = simdat2$survival, id = 'id',longoutcome = 'Y',
   timevarying = c('time','ltime'),
   survtime = 'survtime', cens = 'cens',time = 'time')

   #set variables to factors
   jointdat2$baseline$study <- as.factor(jointdat2$baseline$study)
   jointdat2$baseline$treat <- as.factor(jointdat2$baseline$treat)

   #fit multi-study joint model
   #note: for demonstration purposes only - max.it restricted to 5
   #model would need more iterations to truely converge
   onestagefit<-jointmeta1(data = jointdat2, long.formula = Y ~ 1 + time +
                           + treat + study, long.rand.ind = c('int', 'time'),
                           long.rand.stud = c('treat'),
                           sharingstrct = 'randprop',
                           surv.formula = Surv(survtime, cens) ~ treat,
                           study.name = 'study', strat = TRUE, max.it=5)

    ## Not run: 
        #calculate the SE
        onestagefitSE <- jointmetaSE(fitted = onestagefit, n.boot = 200)

        #print the results of the bootstrap function
        print(onestagefitSE)
    
## End(Not run)

mesudell/joineRmeta documentation built on Jan. 24, 2020, 6:06 p.m.