X-Git-Url: https://code.communitydata.science/ml_measurement_error_public.git/blobdiff_plain/ce6d12d4e6955ba501c9b7ca6715646cf52a5815..d900230ebc453f44afb78f8bc9b12363788e8bfe:/simulations/measerr_methods.R

diff --git a/simulations/measerr_methods.R b/simulations/measerr_methods.R
index 98ab28d..8f8eb5d 100644
--- a/simulations/measerr_methods.R
+++ b/simulations/measerr_methods.R
@@ -23,7 +23,7 @@ likelihood.logistic <- function(model.params, outcome, model.matrix){
 }
 
 ## outcome_formula <- y ~ x + z; proxy_formula <- w_pred ~ y + x + z + x:z + x:y + z:y 
-measerr_mle_dv <- function(df, outcome_formula, outcome_family=binomial(link='logit'), proxy_formula, proxy_family=binomial(link='logit'),maxit=1e6, method='optim'){
+measerr_mle_dv <- function(df, outcome_formula, outcome_family=binomial(link='logit'), proxy_formula, proxy_family=binomial(link='logit'),maxit=1e6, method='optim',optim_method='L-BFGS-B'){
     df.obs <- model.frame(outcome_formula, df)
     proxy.model.matrix <- model.matrix(proxy_formula, df)
     proxy.variable <- all.vars(proxy_formula)[1]