R/adjLibrarySize.R

In SCORPION: Single Cell Oriented Reconstruction of PANDA Individual Optimized Networks

adj_library_size <- function(X, librarySize = colSums(X)) {
  # Design matrix
  H <- model.matrix(~librarySize)

  # Solve coefficients: beta = (H'H)^(-1) H' X'
  HtH_inv <- solve(crossprod(H))
  HtX <- crossprod(H, t(X))
  beta <- HtH_inv %*% HtX

  # Correction: H * beta, then subtract
  X <- t(t(X) - H %*% beta)

  # Clear memory
  gc()
  return(X)
}