fixed bug for MVP.Data

xiaolei-lab · Sep 24, 2024 · f31f5f2 · f31f5f2
1 parent f2dc1c9
commit f31f5f2
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 2 deletions.
diff --git a/R/MVP.Data.r b/R/MVP.Data.r
@@ -183,6 +183,8 @@ MVP.Data <- function(fileMVP = NULL, fileVCF = NULL, fileHMP = NULL, fileBed = N
             sep = sep.kin,
             cpu = ncpus
         )
+    }else{
+        K <- NULL
     }
 
     # get pc

diff --git a/R/MVP.PCA.r b/R/MVP.PCA.r
@@ -60,7 +60,7 @@ function(M=NULL, K=NULL, maxLine=10000, ind_idx=NULL, mrk_idx=NULL, pcs.keep=5,
     if(is.null(K)){
         K <- MVP.K.VanRaden(M=M, ind_idx = ind_idx, mrk_idx = mrk_idx, maxLine = maxLine, cpu = cpu, verbose = verbose)
     }else{
-        K <- K[ind_idx, ind_idx]
+        if(!is.null(ind_idx))    K <- K[ind_idx, ind_idx]
     }
 
     logging.log("Eigen Decomposition on GRM", "\n", verbose = verbose)

diff --git a/src/mvp_omp.h b/src/mvp_omp.h
@@ -17,7 +17,6 @@
 
 #if defined(_OPENMP)
 #include <omp.h>
-// [[Rcpp::plugins(openmp)]]
 #else
 #endif
-Original file line number
+Diff line change
@@ Expand Up / @@ -17,7 +17,6 @@ @@
     #if defined(_OPENMP)
     #include <omp.h>
-    // [[Rcpp::plugins(openmp)]]
     #else
     #endif
@@ Expand Down @@