Skip to content

Commit

Permalink
added a forgotten )
Browse files Browse the repository at this point in the history
  • Loading branch information
borishejblum authored Jan 10, 2025
1 parent 79e870b commit a43f7f5
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion R/cit_gsa.R
Original file line number Diff line number Diff line change
Expand Up @@ -429,7 +429,7 @@ cit_gsa <- function(M,
temp <- indi_pi_gs_tab - matrix(prop_gs_vec, nrow=n_gs_vec, ncol=n_g_t, byrow=TRUE)

# new prop/new pi computation = the one of the gene set, here it's a matrix
new_prop <- apply(temp, MARGIN=2, function(colvec){colvec%*%temp}/n_gs_vec +
new_prop <- apply(temp, MARGIN=2, function(colvec){colvec%*%temp})/n_gs_vec +
sapply(prop_gs_vec, function(s){s*prop_gs_vec})

Sigma2 <- 1/n * tcrossprod(H) %x% (new_prop - prop_gs_vec %x% t(prop_gs_vec))
Expand Down

0 comments on commit a43f7f5

Please sign in to comment.