Fixed error with primersearch and primersearch_raw when there is a single match to a single query

Author: zachary-foster

NEWS.md

@@ -3,6 +3,8 @@
 ## Current
 
 * When an ordered factor is supplied to the `groups` argument of `compare_groups`, the order of levels is used to arrange the results. This can be used to change the order of groups when the output is used with `heat_tree_matrix` (issue [#323](https://github.com/grunwaldlab/metacoder/issues/323)).
+* Added `calc_diff_abund_deseq2` function to use DESeq2 to produce output like `compare_group` useful for `heat_tree_matrix`
+* Fixed error with `primersearch` and `primersearch_raw` when there is a single match to a single query (issue [#326](https://github.com/grunwaldlab/metacoder/issues/326))
 
 ## metacoder 0.3.5 
 

R/primersearch.R

@@ -89,7 +89,7 @@ parse_primersearch <- function(file_path) {
                    "\tAmplimer length: ([0-9]+) bp", sep = '\n')
   primer_data <- stringr::str_match_all(primer_chunks, pattern)
   primer_data <- as.data.frame(cbind(rep(names(primer_chunks), vapply(primer_data, nrow, numeric(1))),
-                                     do.call(rbind, primer_data)[, -1]), stringsAsFactors = FALSE)
+                                     do.call(rbind, primer_data)[, -1, drop = FALSE]), stringsAsFactors = FALSE)
   # Reformat amplicon data
   colnames(primer_data) <- c("pair_name", "amplimer", "input", "name", "f_primer", "f_start",
                              "f_mismatch",  "r_primer", "r_start", "r_mismatch")

src/RcppExports.cpp

@@ -5,6 +5,11 @@
 
 using namespace Rcpp;
 
+#ifdef RCPP_USE_GLOBAL_ROSTREAM
+Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
+Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
+#endif
+
 // euclid
 double euclid(NumericVector a, NumericVector b);
 RcppExport SEXP _metacoder_euclid(SEXP aSEXP, SEXP bSEXP) {