diff options
Diffstat (limited to 'Code/knn_classify.R')
| -rw-r--r-- | Code/knn_classify.R | 1 |
1 files changed, 1 insertions, 0 deletions
diff --git a/Code/knn_classify.R b/Code/knn_classify.R index 46df992..995ae7b 100644 --- a/Code/knn_classify.R +++ b/Code/knn_classify.R @@ -39,6 +39,7 @@ rowsum.tau <- dim(X)[2] # the gene's TPM value is at least 1 on average sd.val <- apply(X, 1, sd) sd.tau <- summary(sd.val)[3] # genes whose gene expression varies least are to be filtered index <- rowSums(X) > rowsum.tau & sd.val > 10 +index[is.na(index)] <- FALSE n.train <- dim(X)[2] X.3 <- log(cbind(X[index,], X.2[index,]) + 1) |