R/print.LORgee.R

In multgee: GEE Solver for Correlated Nominal or Ordinal Multinomial Responses

#' @export
print.LORgee <- function(x, ...) {
  cat("Call:\n")
  print(x$call)
  cat("\nNumber of clusters:", x$max.id, "\n")
  cat("Maximum cluster size:", max(x$clusz), "\n")
  cat("\nNumber of categories:", max(x$categories), "\n")
  cat("Number of observations:", x$nobs, "\n")
  cat("\nCoefficients:\n")
  print(x$coefficients)
  cat("\nNumber of iterations:", x$convergence$niter, "\n")
  cat("Algorithm converged:", x$convergence$conv, "\n")
  if (6 <= nrow(x$local.odds.ratios$theta)) {
    cat("\nLocal Odds Ratios Structure[1:6,1:6]\n")
    print(x$local.odds.ratios$theta[1:6, 1:6])
  } else {
    cat("\nLocal Odds Ratios Estimates[1:6,1:6]\n")
    print(x$local.odds.ratios$theta)
  }
  if (!is.null(x$pvalue)) {
    p_value <- format.pval(x$pvalue, eps = 0.0001, scientific = FALSE)
    cat("\np-value of Null model:", p_value, "\n")
    }
}