diff --git a/DESCRIPTION b/DESCRIPTION
index b2a14bf4c69fa8fe95ec9e08fcd039658514b78d..237c3bec3679302ea79e73c46f04026962614159 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,8 +1,8 @@
 Package: splatter
 Type: Package
 Title: Simple Simulation of Single-cell RNA Sequencing Data
-Version: 1.9.3.9006
-Date: 2019-08-13
+Version: 1.9.3.9007
+Date: 2019-08-14
 Author: Luke Zappia
 Authors@R:
     c(person("Luke", "Zappia", role = c("aut", "cre"),
diff --git a/NEWS.md b/NEWS.md
index eec7a01dfc4c2fc8d7177e54ddc9ad32c364a46f..325e91f285453e3e86ecefbbd889e1ad468b2486 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,3 +1,7 @@
+### Version 1.9.3.9007 (2019-08-14)
+
+* Fix bug in selectFit
+
 ### Version 1.9.3.9006 (2019-08-13)
 
 * Add splotchEstimate function
diff --git a/R/splotch-estimate.R b/R/splotch-estimate.R
index 93f7996f40e8fa82f5b9d08d61912ec8dfef2461..7e3230f6e2b5cec2090593ecd84415c6be3e27bf 100644
--- a/R/splotch-estimate.R
+++ b/R/splotch-estimate.R
@@ -164,10 +164,14 @@ selectFit <- function(data, distr, weights = NULL, verbose = TRUE) {
     }
 
     aics <- fitdistrplus::gofstat(fits)$aic
-    selected <- which(aics == min(aics, na.rm = TRUE))
+    # Flatten in case aics is a list
+    aics.flat <- unlist(aics)
+    selected <- which(aics.flat == min(aics.flat, na.rm = TRUE))
 
     if (verbose) {
-        message("Selected ", names(fits)[selected], " fit using AIC")
+        # Work around to get name in case aics is a list
+        name <- names(fits)[names(aics) == names(aics.flat)[selected]]
+        message("Selected ", name, " fit using AIC")
     }
 
     return(fits[[selected]])