R/logitor.R

Defines functions logitor

Documented in logitor

logitor <-
function(formula, data, robust = FALSE, clustervar1 = NULL, 
                   clustervar2 = NULL, start = NULL, control = list()){
  res = logitorest(formula, data, robust, clustervar1, clustervar2, start, control)
  
  est = NULL
  zstat = log(res$or$oddsratio)*res$or$oddsratio/res$or$se
  
  est$oddsratio = cbind(OddsRatio = res$or$oddsratio,
                     StdErr = res$or$se,
                     z.value = zstat,
                     p.value = 2*pt(-abs(zstat), df = Inf))

  colnames(est$oddsratio) = c("OddsRatio","Std. Err.","z","P>|z|")
  rownames(est$oddsratio) =  rownames(res$or)
  est$fit = res$fit
  est$call = match.call() 
  class(est) = "logitor"
  est
}

Try the mfx package in your browser

Any scripts or data that you put into this service are public.

mfx documentation built on May 2, 2019, 12:46 p.m.