Patch a bug in is.subcladeof

src/robitax.c

@@ -622,7 +622,12 @@ SEXP R_is_under_taxon(SEXP Rtaxonomy, SEXP Rtaxid, SEXP Rparent)
               prep[i]=NA_LOGICAL;
             else {
               taxon = eco_findtaxonbytaxid(ptax, ptaxid[i]);
-              prep[i] = eco_isundertaxon(taxon, parent);            
+              if (taxon) {
+                prep[i] = eco_isundertaxon(taxon, parent);
+              }
+              else {
+                prep[i]=NA_LOGICAL;
+              }
             }  
           }          
         }