Browse code

na.rm=T

guangchuang yu authored on 16/08/2017 11:34:23
Showing 1 changed files

... ...
@@ -144,12 +144,12 @@ get_clade_position_ <- function(data, node) {
144 144
     y <- sp.df$y
145 145
 
146 146
     if ("branch.length" %in% colnames(data)) {
147
-        xmin <- min(x)-data[i, "branch.length"]/2
147
+        xmin <- min(x, na.rm=TRUE)-data[i, "branch.length"]/2
148 148
     } else {
149
-        xmin <- min(sp.df$branch)
149
+        xmin <- min(sp.df$branch, na.rm=TRUE)
150 150
     }
151 151
     data.frame(xmin=xmin,
152
-               xmax=max(x),
153
-               ymin=min(y)-0.5,
154
-               ymax=max(y)+0.5)
152
+               xmax=max(x, na.rm=TRUE),
153
+               ymin=min(y, na.rm=TRUE) - 0.5,
154
+               ymax=max(y, na.rm=TRUE) + 0.5)
155 155
 }