andrjohns
diff --git a/‎R/laplace.R‎
Lines changed: 2 additions & 2 deletions b/‎R/laplace.R‎
Lines changed: 2 additions & 2 deletions
@@ -78,7 +78,7 @@ stan_laplace <- function(fn, par_inits = NULL, n_pars = NULL, additional_args =
   mode_file <- paste0(inputs$output_basename, "_mode.json")
   if (!is.null(mode)) {
     if (inherits(mode, "StanOptimize")) {
-      mode_vals <- mode@estimates[, -1] # First estimate is lp__
+      mode_vals <- mode@estimates[, -c(1,2)] # ignore lp__ & converged__
     } else {
       mode_vals <- mode
     }
@@ -98,7 +98,7 @@ stan_laplace <- function(fn, par_inits = NULL, n_pars = NULL, additional_args =
       sig_figs = sig_figs
     )
     opt <- do.call(stan_optimize, c(curr_args, opt_args))
-    mode_vals <- opt@estimates[, -1]
+    mode_vals <- opt@estimates[, -c(1,2)]
   }
   mode_vals <- as.numeric(mode_vals)
   if (length(mode_vals) != length(inputs$inits[[1]])) {
Original file line number	Diff line number	Diff line change
`@@ -78,7 +78,7 @@ stan_laplace <- function(fn, par_inits = NULL, n_pars = NULL, additional_args =`
`78`	`78`	`mode_file <- paste0(inputs$output_basename, "_mode.json")`
`79`	`79`	`if (!is.null(mode)) {`
`80`	`80`	`if (inherits(mode, "StanOptimize")) {`
`81`		`- mode_vals <- mode@estimates[, -1] # First estimate is lp__`
	`81`	`+ mode_vals <- mode@estimates[, -c(1,2)] # ignore lp__ & converged__`
`82`	`82`	`} else {`
`83`	`83`	`mode_vals <- mode`
`84`	`84`	`}`
`@@ -98,7 +98,7 @@ stan_laplace <- function(fn, par_inits = NULL, n_pars = NULL, additional_args =`
`98`	`98`	`sig_figs = sig_figs`
`99`	`99`	`)`
`100`	`100`	`opt <- do.call(stan_optimize, c(curr_args, opt_args))`
`101`		`- mode_vals <- opt@estimates[, -1]`
	`101`	`+ mode_vals <- opt@estimates[, -c(1,2)]`
`102`	`102`	`}`
`103`	`103`	`mode_vals <- as.numeric(mode_vals)`
`104`	`104`	`if (length(mode_vals) != length(inputs$inits[[1]])) {`