Skip to content
Snippets Groups Projects
Commit a3d19d22 authored by Luke Zappia's avatar Luke Zappia
Browse files

Fix bug in selectFit

Handle missing (not NA) AIC values
parent dd9803d3
No related branches found
No related tags found
No related merge requests found
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"),
......
### Version 1.9.3.9007 (2019-08-14)
* Fix bug in selectFit
### Version 1.9.3.9006 (2019-08-13)
* Add splotchEstimate function
......
......@@ -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]])
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment