Extract Coefficients from robcov_vglm Objects
Usage
# S3 method for class 'robcov_vglm'
coef(object, ...)Examples
trial <- sim_markov("actt2", n_patients = 40, follow_up_time = 6, seed = 1)
markov_data <- prepare_markov_data(trial, absorbing_state = 8)
fit <- vglm_markov(
ordered(y) ~ time + tx + yprev,
family = VGAM::cumulative(reverse = TRUE, parallel = TRUE),
data = markov_data
)
#> Warning: `id_var` was not supplied to vglm_markov(); stored patient-ID metadata, cluster-robust standard errors for frequentist wrappers, and automatic FWB refits will not be available from the fitted object.
#> Warning: 12 diagonal elements of the working weights variable 'wz' have been replaced by 1.819e-12
#> Warning: 76 diagonal elements of the working weights variable 'wz' have been replaced by 1.819e-12
#> Warning: 77 diagonal elements of the working weights variable 'wz' have been replaced by 1.819e-12
#> Warning: 128 diagonal elements of the working weights variable 'wz' have been replaced by 1.819e-12
#> Warning: 151 diagonal elements of the working weights variable 'wz' have been replaced by 1.819e-12
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: 151 diagonal elements of the working weights variable 'wz' have been replaced by 1.819e-12
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: 151 diagonal elements of the working weights variable 'wz' have been replaced by 1.819e-12
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: 151 diagonal elements of the working weights variable 'wz' have been replaced by 1.819e-12
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: fitted values close to 0 or 1
#> Warning: It seems that the nonparallelism assumption has resulted in intersecting linear/additive predictors. Try propodds() or fitting a partial nonproportional odds model or choosing some other link function, etc.
#> Warning: some quantities such as z, residuals, SEs may be inaccurate due to convergence at a half-step
robust_fit <- robcov_vglm(fit, cluster = markov_data$id)
coef(robust_fit)
#> (Intercept):1 (Intercept):2 (Intercept):3 (Intercept):4 (Intercept):5
#> -0.2507381 -21.5699735 -21.8089291 -24.2885588 -28.7619954
#> (Intercept):6 time tx yprev2 yprev3
#> -33.3092327 -0.2626019 -0.1647912 -0.8884239 11.8641090
#> yprev4 yprev5 yprev6 yprev7
#> 23.9845226 26.5162037 32.2500878 37.2189589