update real data examples code and rerun project.

[ml_measurement_error_public.git] / simulations / pl_methods.R

diff --git a/simulations/pl_methods.R b/simulations/pl_methods.R
index b3007d12579dde9d358fe850fb2b86def89ed7b2..f014eec1b8f00479417ad9181d6c6bc5127a6af0 100644 (file)
--- 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)
         (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),
         return(-ll)
     }    
     mlefit <- mle2(minuslogl = nll, control=list(maxit=1e6),method='L-BFGS-B',lower=c(B0=-Inf, Bxy=-Inf, Bzy=-Inf),