https://github.com/cran/robCompositions
Tip revision: 54cbcf1d6d3ecbca1153e8f00a804e837237f09d authored by Matthias Templ on 07 February 2014, 00:00:00 UTC
version 1.7.0
version 1.7.0
Tip revision: 54cbcf1
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(as.numeric(x[complete.cases(x)]))
})
gmI <- apply(object$xImp, 2, function(x) {
geometricmean(as.numeric(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
}