R/mcov.R

In gim: Generalized Integration Model

mcov <- function(para, map, family, data, ref, model, sample.info, V, bet0, outcome, type){
  
  if(family == 'gaussian'){
    mat <- cov.lm(para, map, data, ref, model, sample.info, V, bet0, outcome)
  }
  
  if(family == 'binomial'){
    mat <- cov.lo(para, map, data, ref, model, sample.info, V, bet0, outcome)
  }
  
  if(family == 'case-control'){
    if(type == 'cc-ref'){
      mat <- cov.ccr(para, map, data, ref, sample.info, V, bet0, outcome)
    }else{
      mat <- cov.cc(para, map, data, ref, sample.info, V, bet0, outcome)
    }
  }
  
  if(family == 'cml'){
    mat <- cov.cml(para, map, data, ref, sample.info, V, bet0, outcome)
  }
  
  mat
  
}