https://github.com/cran/robCompositions
Raw File
Tip revision: 0c00ba572d6b53d50c6deb186fd64191d28834cd authored by Matthias Templ on 15 December 2008, 00:00:00 UTC
version 1.0.1
Tip revision: 0c00ba5
adjust.R
adjust <- function(x){
	# x ... object from class "imp"
	if(class(x) != "imp") stop("object x must be from class imp")
	xneu=xi$xImp
	s1 <- rowSums(xi$xOrig, na.rm=TRUE)
	for(i in 1:nrow(xi$xImp)){
		s <- sum(xi$xImp[i, !xi$wind[i,]])
		s2 <- sum(xi$xImp[i, xi$wind[i,]])
		fac <- s / (s + s2)
		s1[i] <-  s1[i] / fac
	}
	impS <- s1/rowSums(xi$xImp)
	for(i in 1:ncol(xi$xImp)){
		xneu[,i] <- xi$xImp[,i] * impS
	}
	x$xImp <- xneu
	invisible(x)
}
back to top