rmReact: faild if resulting mod had only one gene left.

R/rmReact.R

@@ -169,7 +169,7 @@ rmReact <- function(model, react, rm_met = TRUE) {
       keepGenes <- ifelse(is.na(keepGenes), FALSE, TRUE)
       #print(keepGenes)
 
-      rxnGeneMat(mod_out)   <- rxnGeneMat(mod_out)[, keepGenes]
+      rxnGeneMat(mod_out)   <- rxnGeneMat(mod_out)[, keepGenes, drop = FALSE]
       #print(dim(rxnGeneMat))
   }
   

inst/NEWS.Rd

@@ -8,6 +8,7 @@
 
 \section{Changes in version 1.3.1 2015-08-06}{
   \itemize{
+  	\item \code{rmReact} error, if resulting model had only one gene left, corrected.
     \item \code{deadEndMetabolite()} could miss deadEndMetabolites if reactions differ in stoichiometry e.g. 1 and 2.
   }
 }