R/print.outlierTest.R

In robmixglm: Robust Generalized Linear Models (GLM) using Mixtures

print.outlierTest <- function(x, ...) {
  if (!inherits(x, "outlierTest"))
    stop("Use only with 'outlierTest' objects.\n")
  thestat <- x$t0[1]
  thestat <- ifelse(thestat<0.0,0.0,thestat)
  nullstat <- x$t[,1]
  nullstat <- ifelse(nullstat<0.0,0.0,nullstat)
  pos <- sum(nullstat>thestat)
  test <- (pos+1)/(length(x$t[,1])+1)
  out <- sprintf("p value %1.4f", test)
  cat(out)
  invisible(out)
}

summary.outlierTest <- function(object,...) {
  if (!inherits(object, "outlierTest"))
    stop("Use only with 'outlierTest' objects.\n")
  class(object) <- "summary.outlierTest"
  return(object)
}

print.summary.outlierTest <- function(x,...) {
  if (!inherits(x, "summary.outlierTest"))
    stop("Use only with 'summary.outlierTest' objects.\n")
  thestat <- x$t0[1]
  thestat <- ifelse(thestat<0.0,0.0,thestat)
  nullstat <- x$t[,1]
  nullstat <- ifelse(nullstat<0.0,0.0,nullstat)
  pos <- sum(nullstat>thestat)
  test <- (pos+1)/(length(x$t[,1])+1)
  out <- sprintf("p value %1.4f", test)
  cat(out)
  invisible(out)
}