Minor update - fix a potential error using R 4.3.0+ where arguments a…

…round "&&" must have the same length
VanLoo-lab · Apr 27, 2023 · c4c10f8 · c4c10f8
1 parent 1cca6cc
commit c4c10f8
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/ASCAT/R/ascat.aspcf.R b/ASCAT/R/ascat.aspcf.R
@@ -173,8 +173,9 @@ ascat.aspcf = function(ASCATobj, selectsamples = 1:length(ASCATobj$samples), asc
         if(length(homsegs)==3) {
           homsegs=t(as.matrix(homsegs))
         }
-        if(!is.null(homsegs)&&!is.na(homsegs)&&dim(homsegs)[1]!=0) {
+        if(!is.null(homsegs)&&dim(homsegs)[1]!=0) {
           for (i in 1:dim(homsegs)[1]) {
+            if (anyNA(homsegs[i,])) next
             # note that only the germline homozygous segment is resegmented, plus a bit extra (but this is NOT replaced)
             startpos = max(homsegs[i,2],startchr)
             endpos = min(homsegs[i,3],endchr)