https://github.com/cran/robCompositions
Tip revision: e0e46e2d22f0b04308b896fd000cd7f493f04383 authored by Matthias Templ on 16 August 2010, 00:00:00 UTC
version 1.4.3
version 1.4.3
Tip revision: e0e46e2
adjust.R
adjust <- function(x){
# x ... object from class "imp"
if(class(x) != "imp") stop("object x must be from class imp")
xneu=x$xImp
s1 <- rowSums(x$xOrig, na.rm=TRUE)
for(i in 1:nrow(x$xImp)){
s <- sum(x$xImp[i, !x$wind[i,]])
s2 <- sum(x$xImp[i, x$wind[i,]])
fac <- s / (s + s2)
s1[i] <- s1[i] / fac
}
impS <- s1/rowSums(x$xImp)
for(i in 1:ncol(x$xImp)){
xneu[,i] <- x$xImp[,i] * impS
}
x$xImp <- xneu
invisible(x)
}