similar to: Integration with "adapt"

R-listers, In using the adapt function, I am getting the following warning: Ifail=2, lenwrk was too small. -- fix adapt() ! Check the returned relerr! in: adapt(ndim = 2, lower = lower.limit, upper = upper.limit, functn = pr.set, Would someone explain what the 'lenwrk' value indicates in order to help diagnose this issue. Also, what are the possible codes for Ifail, so I can set

Hi all, I found a problem on R that looks like a bug to me. I am working with R for windows verson 1.9.0, but the problem also happens on the linux version 1.9.0. I hope that the problem is on my reasoning, and I would be grateful if someone could explain me why this happens. There goes the example: > x<-seq(-2,2,0.001) > y<-rep(0,length(x)) > for (i in x){y[1000*i+2001]<-i^2}

Hi all I was wondering if someone could help me. I have to estimate some parameters, so I am using the function nlm. Inside this function I have to integrate, hence I am using the function adapt. I don't understand why it is giving the following warnings: At the beginning: Warning: a final empty element has been omitted the part of the args list of 'c' being evaluated was:

Dear People, I'm trying to use the function adapt, from the adapt library package, which does multidimensional numerical integration. I think I must be using the wrong syntax or something, because even a simple example does not work. Consider foo <- function(x){x[1]*x[2]} and adapt(2, lo = c(-1,-1), up = c(1,1), functn = foo) This simply hangs. A more complicated example crashes R,

Hello Bug people, I have an unexpected behavior and am unsure whether the problem is in my thinking, my implementation or the program R. Basically I get two different answers depending on how I call a function which takes other functions as arguments as indicated below. To me it should make no difference if f is a function that returns the function g then z(f(x)) whould give the same as y<-

I'm having trouble with adapt. I'm trying to use it in a Bayesian setting, to integrate the posterior distribution, and to find posterior means. I tried using the following script, and things went ok: data = rnorm(100,0.2,1.1) data = c(data,rnorm(10,3,1)) data = data[abs(data)<2*sd(data)] prior = function(x){ dgamma(x[2],shape=2,scale=1)*dnorm(x[1],0,.5) } liklihood =

I tried to integrate numerically a function wich is similar to the following: > dummy <- function(x) { exp(-1*x) * dnorm(x) } > dummy(-100) [1] 0 > dummy(-1000) [1] NaN > dummy(-10000) [1] NaN If I choose the lower boundary to be too small integrate causes a segmentation fault: > library(integrate) > integrate(dummy, -100, 0)$value [1] 1.387143 > integrate(dummy, -1000,

Dear r-users, Is there any command in R allowing to evaluate a double integral? for instance let say I want to evaluate the following integral: integrate[lo=(0,1),up=(2,3)] f(x,y)=x^2+y^2 where lo is the vector of lower bounds and up that of upper bounds. I thaught the function "adapt" would work but it did not. Many thanks, Dominique K.

how does one numerically intergate the following: A=function(x,y) { xy } over the range: 2<x<0 4<y<10 say. ie how would one set up the integrate function? i forgot!

Hi, I need to integrate a 2D function over range where the limits depend on the other e.g integrate f(x,y)=x*y over {x,0,1} and {y,x,1}. i.e \int_0^1 \int_x^1 xy dydx I checked adapt but it doesn't seem to help here. Are they any packages for this sort of thing? I tried RSitesearch but couldn't find the answer to this. Many thanks for you help. Regards Saptarshi Saptarshi Guha

Hi all, Running R 0.90.1 on a RH 6.1 system. Installation of the integrate_2.1-2 package went smoothly. My code contains a loop in which integrate() is called several times in each pass. I get a segmentation fault after what seems to be a random number of calls to integrate(). Debug output shows: Program received signal SIGSEGV, Segmentation fault. promiseArgs (el=0x40276414,

Hi all, Does anyone know of a function that return the best lambda in box-cox transformation? many thanks for any help. -- Cristiane S. Rocha Laboratorio Genoma Funcional - Bioinform?tica Centro de Biologia Molecular e Engenharia Genetica Universidade Estadual de Campinas Campinas - SP - Brasil Tel:(19)3788-1119

Hi all. Running 2.0.4b on a (now) solaris 2.6 box. Since the upgrade the smbd & nmbd daemons won't start. I tried starting smbd by hand, running it through truss. Here's the last few lines: getuid() = 0 [0] fstat(3, 0xEFFFF858) = 0 write(3, " [", 1) = 1 write(3, " 1 9 9 9 / 0

I get the following errors on the console of my newly installed 2.0.4b server on a Sun E4000: nmbd/nmbd_packets.c: (1412) find_subnet_for_nmb_packet: response record not found for response id <num>. nmbd/nmbd_responserecordsdb.c: (240) find_response_record: response packet id <num> received with no matching record I don't know the internals well enough to decipher these, can

Hello, I'm having a error message when I try to load the som library from a cgi. When a run my script via web the following error appears in the apache log: Error: package som was built for i686-pc-linux-gnu, referer: https:// .... I call the library function setting the path to the library directory: >library(som, lib.loc="/usr/local/lib/R/library/") The script runs fine

Hi all. Is it possible to have a generic share like [homes] that hides all user home dir shares by default by appending a "$" to them? I'd search the archives but they have not yet been tx'd. I realise there is "browse vs nobrowse" but I specifically need to have end users map to \\srv\user$ rather than \\srv\user. -- Mark

I’ve encountered this kind of architecture before, a long time ago (academically). In a flat-address-space machine such as X64, there is still an instruction/data distinction, but usually only down at the level of I-cache versus D-cache (instruction fetch versus data fetch). A Harvard architecture machine exposes that to the programmer, which effectively doubles the available address space.

On Wed, 15 Jan 2020, Reid Kleckner via llvm-dev wrote: > On Wed, Jan 15, 2020 at 11:14 AM Luke Drummond <luke.drummond at codeplay.com> > wrote: > On Wed Jan 15, 2020 at 6:58 PM, Reid Kleckner wrote: > > Does the same limitation exist in VS 2017? I think that's our > support > > floor > > these days: >

I'm not across most of this debug info stuff but I'll stomp in here to confirm that AVR is a Harvard architecture, with separate addressing for the data and program buses via specialized instructions which will load from either one, or the other, but never both. It makes sense that this particular problem would also affect AVR - the backend does have some issues with debug info

I have no strong opinion on this topic personally (I always set the build type to something if I’m not generating a Visual Studio project), but this “non-optimized, non-asserting, non-debug-info-containing” build that you get when not specifying the build type to CMake is a quirk that I feel most CMake users expect at this point. I don’t think there’s any real value in providing a useful default