search for: lowerlimit

...length(dose)) { psi0<-1/((1+exp(x[1]+x[2]*dose[i]))*(1+exp(x[3]+x[4]*dose[i]))) psi1<-exp(x[1]+x[2]*dose[i])/((1+exp(x[1]+x[2]*dose[i]))*(1+exp(x[3]+x[4]*dose[i]))) v<-v*(psi0^y0[i])*(psi1^y1[i])*((1-psi0-psi1)^y2[i]) } return(v) } adaptIntegrate(lf, lowerLimit = c(-20, 0,-20,0), upperLimit = c(0,10,0,10)) --------------------------------------------------------------------------------------------------------------------------------- Thanks for your attention. Kind regards, Jamil. [[alternative HTML version deleted]]

Hello Dear List members, as you can see (and guess) from the code below adaptIntegrate(f,lowerLimit=c(-1,-1),upperLimit=c(.9999,.9999)) $integral [1] 9.997e-09 $error [1] 1.665168e-16 $functionEvaluations [1] 17 $returnCode [1] 0 > adaptIntegrate(f,lowerLimit=c(-1,-1),upperLimit=c(1,1)) the last command runs for 45 mins now. -this one takes only less than sec: adaptIntegrate(f,lowerLim...

...om Cubature to do numerical integration on a double integral with a 1 x 2 vector x. Say the function is something simple to start like f(x)=x1*x2 and I wish to integrate x1 over (0,365-x2) and x2 over (0,365) f <- function(x) {(x[2])*(x[1])} # "x" is vector int1<-adaptIntegrate(f, lowerLimit = c(0, 0), upperLimit = c(365-x[2], 365)) I recieve the following error: Error in adaptIntegrate(f, lowerLimit = c(0, 0), upperLimit = c(365 - : object 'x' not found The problem is that this code works fine on other machines and I wonder If I am missing some referenced package somew...

...f(x1,x2)=2/3*(x1+x2) in the interval 0<x1<x2<7. To be sure I tried it by hand and got 114.33, but the following R code is giving me 102.6667. ------------------------------------------------------------------- library(cubature) f<-function(x) { 2/3 * (x[1] + x[2] ) } adaptIntegrate(f, lowerLimit = c(0, 0), upperLimit = c(x[2],7)) ------------------------------------------------------------------- I guess something is wrong with the way I have assigned limits in my codes. May I seek some advice from you. Many thanks. Regards, Jamil. [[alternative HTML version deleted]]

I have sent a similar email about this problem. But I did not describe carefully, so I would like to explain it in more details. I am using 'virsh setmem' to ajust vm memory online. However, I don't know what is the lowerlimit that can be set to. And I try to use 'virsh dommemstat' to get 'unused' memory so that I can calculate the lowerlimit memory with this value, but it doesn't work when it comes to windows systems as I cannot set vm memory down to 'available - unused' in windows system. I...

...the integrand also calls the function adaptIntegrate. In particular I want \int \hat{f}(x,y) - f(x,y) dx dy where \hat{f}(x,y) = \int K(a,b, x, y) da db and in this simulation study I know what the true value of f(x,y) is. I can use adaptIntegrate to calculate \hat{f}(x1,y1) = adaptIntegrate( K, lowerLimit=c(a1,b1), upperLimit=c(a2,b2), x=x1, y=y1)$integral which gives a real number. However if I define the following integrand and try to use adaptIntegrate on it I run into problems: integrand <- function(x1,y1){ return(adaptIntegrate( K, lowerLimit=c(a1,b1), upperLimit=c(a2,b2), x=x1, y=y1)$...

I am using 'virsh setmem' to adjust the vm memory. But I don't know what is the lowerlimit memory that can be adjusted to, so I try to update the version of libvirt&amp;qemu to get caches of the vm memory so that I can calculate the lowerlimit. And when I upgrade libvirt to 4.10.0, qemu upgrade to 4.1.0, still, I cannot get 'stat-disk-caches'.&nbsp; How can I solve it ? T...

Dear list, I'm seeking some advice regarding a particular numerical integration I wish to perform. The integrand f takes two real arguments x and y and returns a vector of constant length N. The range of integration is [0, infty) for x and [a,b] (finite) for y. Since the integrand has values in R^N I did not find a built-in function to perform numerical quadrature, so I wrote my own after

Hi All, I'm trying to use one of the (2D) numerical integration functions, which is not where the problem is. The function definition is as follows: adaptIntegrate(f, lowerLimit, upperLimit, ...) The problem is that I want to integrate a 3D function which has been parametrised such that it is a 2D function: eg. I want to integrate f(x_start+gradient_x*t1, y_start+gradient_y*t2) for t1 in [0,1] and t2 in [0,1]. To explain the example in another way I only have access to f...