R/vcov.MMLongit.R

In mercaldo/MMLB: Marginalized Models for Longitudinal Binary Data

vcov.MMLongit <- function (object) {
  n_beta     = length(object$beta)
  n_alpha    = length(object$alpha)
  beta_ix    = seq(n_beta)
  alpha_ix   = seq(n_beta+1, n_beta+n_alpha)
  object_vc  = object$mod.cov
  if(object$control[7]) object_vc  = object$rob.cov
  list('beta'=object_vc[beta_ix,beta_ix], 'alpha'=object_vc[alpha_ix,alpha_ix])
}