search for: maxeval

...library(cubature) ff <- function(x, rho){ mu <- rep(0,3) Sigma <-(1-rho)*diag(3)+matrix(rho,3,3) f <- dmnorm(x, mu, Sigma) f } adaptIntegrate(ff, lower=c(-10, -10, -10), upper=c(3,2,1), rho=0.1, maxEval=10000) adaptIntegrate(ff, lower=rep(-Inf, 3), upper=c(3,2,1), rho=0.1, maxEval=10000) [[alternative HTML version deleted]]

Dear R-users, I am facing problem with integrating in R a likelihood function which is a function of four parameters. It's giving me the result at the end but taking more than half an hour to run. I'm wondering is there any other efficient way deal with. The following is my code. I am ready to provide any other description of my function if you need to move forward.

...eval_f0 <- function(x) { romberg(f, 0.5, 0.5001)} ARL1 <- nloptr( x0=c(0.653), eval_f=eval_f0, lb = c(0), ub = c(6), eval_g_ineq = eval_g0, opts = list("algorithm"="NLOPT_LN_COBYLA", "maxeval"=1000), ) print( ) Thanks [[alternative HTML version deleted]]

I have a frame which contains 3 columns: "date" "defectnum" "state" And I want to get the most recent state change for a given defect number. date is POSIXct. I have tried: aggregate(ev$date, by=list(ev$defectnum), max) Which appears to be working except that the dates seem to come back as integers (presumably the internal representation of POSIXct). When I

...1), 2,2) + mu <- rep(0,2) + f <- pmnorm(c((z-a*x1)/b, (z-a*x2)/b), mu, Sigma)*dmnorm(c(0,0), mu, diag(2)) + f + } > > adaptIntegrate(integrand, lower=rep(-Inf, 2), upper=c(2,2), + rho=0.1, a=0.6, b=0.3, z=3, maxEval=10000) Error in if (any(lower > upper)) stop("lower>upper integration limits") : missing value where TRUE/FALSE needed [[alternative HTML version deleted]]

...t;- function(x) romberg(function(s) f(x, s), 0, 6) of <- nloptr( x0=c(0.653,0.09), eval_f=eval_f0, lb = c(0,0), ub = c(6,1), eval_g_ineq = eval_g0, opts = list("algorithm"="NLOPT_LN_COBYLA", "maxeval"=1000), ) print(of) Regards, Aya [[alternative HTML version deleted]]

..." COM(6)=G(6,1) " COM(7)=HH(1,1) " COM(8)=HH(2,1) $THETA (0, 1.7) ;KA (0, 4.7) ;CL (0, 38) ;V (0.0, .5) ; EMAX (0.0, 100) ; EC50 (0.0, 1, 10) ; N $OMEGA 0.04 0.04 0.04 0.04 0.04 0.04 $SIGMA 0.04 5 $ESTIMATION MAXEVAL=9999 NOABORT METHOD=0 POSTHOC PRINT=5 $COVARIANCE $TABLE ID TIME AGE BW GNDR TVKA TVV TVCL ETA1 ETA2 ETA3 EVID NOPRINT ONEHEADER FILE=TABLE.PAR $TABLE NOPRINT FILE=../3033.TAB ONEHEADER ID TIME AGE BW GNDR KA V CL TVKA TVCL TVV FLAG FLA2 EVID EFF CONC IPRED FLAG $TABLE NOPRINT ONEHEADER FILE=../30...

...)/6),0,6, tol = .Machine$double.eps^(2/3))$value obj <- nloptr( x0=c(0.653,0.09), eval_f=eval_f0, lb = c(0,0), ub = c(6,1), eval_g_ineq = eval_g0, opts = list("algorithm"="NLOPT_LN_COBYLA", "maxeval"=600)) obj(ARL1) Thanks, Emara [[alternative HTML version deleted]]

Signed-off-by: Jan Beulich <jbeulich@novell.com> Index: head-2007-02-27/drivers/char/mem.c =================================================================== --- head-2007-02-27.orig/drivers/char/mem.c 2007-03-05 10:00:18.000000000 +0100 +++ head-2007-02-27/drivers/char/mem.c 2007-02-27 16:27:37.000000000 +0100 @@ -787,7 +787,7 @@ static const struct file_operations mem_ .open =

...ach variable in vector x for objective function local_opts=list("algorithm" = "NLOPT_LN_AUGLAG", "xtol_rel"=0.01) ## some algorithms in nloptr require a local algorithm too opts=list("algorithm" = "NLOPT_GN_ISRES", "xtol_rel" = 0.01, "maxeval"=100000 ,"local_opts" = local_opts, "print_level"=2) nloptr(x0=x0, eval_f = eval_f, lb=lb, ub=ub, eval_g_eq = eval_g_eq, opts=opts) -- Alicia Ellis Postdoc Gund Institute for Ecological Economics University of Vermont 617 Main Street Burlington, VT 05405 (802) 656-1046...

depends on "autoconf: regenerate configure scripts with 4.4 version" This series removes some of the really old deadwood from the tools build and makes some other things which are on their way out configurable at build time with a default depending on how far down the slope I judge them to be. * nuke in tree copy of libaio * nuke obsolete tools: xsview, miniterm, lomount & sv *

I'm currently working on gobject bindings for libguestfs. I haven't got as far as compiling anything yet, but I've attached the C header for initial review. Matt -- Matthew Booth, RHCA, RHCSS Red Hat Engineering, Virtualisation Team GPG ID: D33C3490 GPG FPR: 3733 612D 2D05 5458 8A8A 1600 3441 EA19 D33C 3490 -------------- next part -------------- An embedded and

This includes two kernel patches, scsi-target and IO_CMD_EPOLL_WAIT. The former is a modified version of the scsi target infrastructure in mainline. The latter enables applications to handle AIO and non-AIO fds in the same loop. blktap uses the different patch, AIO event queue patch for the same aim. The IO_CMD_EPOLL_WAIT patch will be merged into mainline (and the AIO event queue will not) so

[This email is either empty or too large to be displayed at this time]

This is a follow-up to the other generator changes in: https://www.redhat.com/archives/libguestfs/2017-February/msg00217.html Rich.

Split up the huge generator/actions.ml into several smaller files. Rich.

The existing mechanism was clunky, slow and used ~ 10 MB of local disk. Rich.