X-Git-Url: https://code.communitydata.science/ml_measurement_error_public.git/blobdiff_plain/82fe7b0f482a71c95e8ae99f7e6d37b79357506a..c45ea9dfebca86dfddc1e9237aa74866c5166519:/simulations/pl_methods.R diff --git a/simulations/pl_methods.R b/simulations/pl_methods.R index b3007d1..f014eec 100644 --- a/simulations/pl_methods.R +++ b/simulations/pl_methods.R @@ -31,8 +31,8 @@ zhang.mle.dv <- function(df){ (1-w_pred) * (log(1-fpr) - exp(log(1-fnr-fpr)+pi.y.1))))) ll <- ll + sum(lls) - print(paste0(B0,Bxy,Bzy)) - print(ll) +# print(paste0(B0,Bxy,Bzy)) +# print(ll) return(-ll) } mlefit <- mle2(minuslogl = nll, control=list(maxit=1e6),method='L-BFGS-B',lower=c(B0=-Inf, Bxy=-Inf, Bzy=-Inf),