R/computeContCols.R

"computeContCols" <-
function(x,CL,vec.ncl,n.obs){
	x<-x%*%CL
	r<-rowSums(x)
	if(is.matrix(vec.ncl))
		tmp <- r * vec.ncl / n.obs
	else 
		tmp <- r %*% t(vec.ncl) / n.obs
	rowSums(x*x/tmp)
}

Try the siggenes package in your browser

Any scripts or data that you put into this service are public.

siggenes documentation built on May 6, 2019, 2:38 a.m.