DA/raaql-paper-experiments Files · test.R · Centrum Wiskunde & Informatica (CWI)

Files @ 5cabbae4fdd5

Branch filter:

Location: DA/raaql-paper-experiments/test.R

5cabbae4fdd5 387 B text/S-plus Show Annotation Show as Raw Download as Raw

Hannes Muehleisen

submitted

s <- 6

wts <- readRDS(paste0("../paper-old/wts", s, ".rds"))
print(str(wts))
x <- readRDS(paste0("../paper-old/x", s, ".rds"))
dim(x) <- c(length(x), 1)
pw <- 1L

print(system.time({
	repmeans<-matrix(data=as.numeric(NA), ncol=NCOL(x), nrow=ncol(wts))
	for(i in 1:ncol(wts)){
	  repmeans[i,]<-t(colSums(wts[,i]*x*pw)/sum(pw*(wts[,i])))
	}
}))

print(system.time({
	print(repmeans)
}))