R/or.table.R

In descriptio: Descriptive Statistical Analysis

# or.table(x0, y0, weights = w0, na.rm = FALSE)
# or.table(x0, y0, weights = w0, na.rm = TRUE)
# or.table(x0, y0, weights = w1, na.rm = FALSE)
# or.table(x0, y0, weights = w1, na.rm = TRUE)
# or.table(x1, y0, weights = w0, na.rm = FALSE)
# or.table(x1, y0, weights = w0, na.rm = TRUE)

or.table <- function(x, y, weights = NULL, na.rm = FALSE, na.value = "NA", digits = 3) {
  t0 <- weighted.table(x, y, weights, stat = "freq", mar = FALSE, na.rm = na.rm, na.value = na.value)
  OR <- t0
  for(i in 1:nrow(t0)) {
    for(j in 1:ncol(t0)) {
      a <- t0[i,j]
      b <- rowSums(t0)[i] - a
      c <- colSums(t0)[j] - a
      d <- sum(t0) - (a+b+c)
      OR[i,j] <- (a*d) / (b*c)
    }
  }
  if(!is.null(digits)) OR <- round(OR,digits)
  return(OR)
}