vcov.jointmeta1SE: Extract the variance covariance matrix from the bootstrapped...

Description Usage Arguments Value See Also Examples

View source: R/vcov.jointmeta1SE.R

Description

Function applied to a jointmeta1SE object, the result of the jointmetaSE function to extract the variance covariance matrix for the estimated model parameters

Usage

1
2
## S3 method for class 'jointmeta1SE'
vcov(object, ...)

Arguments

object

an object of class jointmeta1SE

...

additional arguments; currently none are used.

Value

a variance covariance matrix for the fixed effects from the longitudinal sub-model, the time-to-event sub-model, the association parameters, the random effects and the error term.

See Also

jointmeta1, jointmetaSE, jointmeta1SE.object

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
   #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)

        #extract the variance covariance matrix
        vcov(onestagefitSE)
    
## End(Not run)

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