https://github.com/cran/robCompositions
Raw File
Tip revision: 37555063454f14ae40f041b3a8f0eb2d3dc3be4c authored by Matthias Templ on 10 January 2011, 00:00:00 UTC
version 1.4.4
Tip revision: 3755506
summary.imp.R
summary.imp <- function(object, ...){
  geometricmean <- function (x) {
    if (any(na.omit(x == 0)))
        0
    else exp(mean(log(unclass(x)[is.finite(x) & x > 0])))
  }
gm <- apply(object$xOrig, 2, function(x) {
  geometricmean(x[complete.cases(x)])
})
gmI <- apply(object$xImp, 2, function(x) {
  geometricmean(x[complete.cases(x)]) ## gewichten!
})

  d <- data.frame(orig=gm,
                  imp=gmI)
  cat("\n geometric mean of the original data and the imputed data: \n")

  ## Einfluss der Imputation mittels bootstrap schaetzen
  ## laut script sim

  d
}
back to top