search for: maxev

...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=../...

...)/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]]

...nsigned int aio_queue_fd_poll(st + return pollflags; +} + -+static struct file_operations aioq_fops = { ++static const struct file_operations aioq_fops = { + .release = aio_queue_fd_close, + .poll = aio_queue_fd_poll +}; @@ -201,7 +201,7 @@ static int ep_events_transfer(struct eve int maxevents, long timeout); static int eventpollfs_delete_dentry(struct dentry *dentry); -static struct inode *ep_eventpoll_inode(void); -+static struct inode *ep_eventpoll_inode(struct file_operations *fops); ++static struct inode *ep_eventpoll_inode(const struct file_operations *fops); static int ev...

...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-10...

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

...mp;ep->lock, flags); ++ goto out; ++ } + /* + * We don''t have any available event to return to the caller. + * We need to sleep here, and we will be wake up by +@@ -1562,7 +1713,7 @@ retry: + if (!res && eavail && + !(res = ep_events_transfer(ep, events, maxevents)) && jtimeout) + goto retry; +- ++out: + return res; + } + +diff --git a/include/linux/aio_abi.h b/include/linux/aio_abi.h +index 30fdcc8..bb67d5b 100644 +--- a/include/linux/aio_abi.h ++++ b/include/linux/aio_abi.h +@@ -41,6 +41,7 @@ enum { + * IOCB_CMD_POLL = 5, + */ + IOCB_C...

[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.