R/wiqid_methods.R

Defines functions nobs.wiqid logLik.wiqid print.wiqid vcov.wiqid coef.wiqid

# Various methods for class wigid, ie. output for wiqid MLE functions.

# coef and vcov added 2017-02-09

coef.wiqid <- function(object, ...) {
  object$beta[,1]
}

vcov.wiqid <- function(object, ...) {
  object$beta.vcv
}

print.wiqid <- function(x, digits=4, ...)  {
  cat("Call: ")
  print(x$call)
 if(anyDuplicated(x$real) > 0) {
    cat("\nReal values (duplicates omitted):\n")
    print(unique(x$real), digits=digits, ...)
  } else {
    cat("\nReal values:\n")
    print(x$real, digits=digits, ...)
  }
  cat("\nAIC:", AIC(x), "\n")
}

logLik.wiqid <- function(object, ...)  {
  tmp <- as.vector(object$logLik)
  ll <- tmp[1]
  attr(ll, 'df') <- tmp[2]
  attr(ll, 'nobs') <- tmp[3]
  class(ll) <- "logLik"
  return(ll)
}

nobs.wiqid <- function(object, ...)  {
  nobs <- as.numeric(object$logLik)[3]
  return(if(is.null(nobs)) NA_integer_ else nobs)
}

Try the wiqid package in your browser

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

wiqid documentation built on Nov. 18, 2022, 1:07 a.m.