similar to: double integral

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

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,

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:

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

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, I need to calculate a multidimensional integration on R. I am using the command "adapt" (from library adapt), although sometimes I get the following error message: Ifail=2, lenwrk was too small. -- fix adapt() ! Check the returned relerr! in: adapt(3, linf, lsup, functn = Integrando1) I guess it happens because the domain of integration is too small, although I tried a

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,

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,

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 =

Hi all I hope to get some help on the separator. I am sure that I hava a wrong understanding what the separator is. I thaught the separator is the delimiter within the imap folder structure. For example: Imap folder structure: folder Flights, subfolder Lufthansa If the delimiter is ?.? the ls ?a command should result: .Flights .Flights.Lufthansa If delimiter is

Hi! I'm learning Mathematical Statistics at a hungarian university. My teacher, who teaches these lessons, using the R to introduce some topics. I want to pass on the exam, and this exam contains some R questions too, but the teacher doesn't thaught us too much about R, so I would like to ask for some help from you in solving some example questions in R. So, here are some

On Fri, 22 Feb 2019 at 17:46, Rowland Penny via samba <samba at lists.samba.org> wrote: > On Fri, 22 Feb 2019 17:41:50 +0100 > Hans Schou via samba <samba at lists.samba.org> wrote: > > > On Fri, 22 Feb 2019 at 17:27, Rowland Penny via samba > > <samba at lists.samba.org> wrote: > > > > > > > > [global] > > > > log

Hi , I picked up a x100p the other day and thaught I'd havea go at getting the driver going for linux 2.6, things have gone pretty, two basic problems. 1. makefiles, with 2.6 you can't get away with using the old makefile to build the kernel modules, they will build but you'll get an error along the lines of "no module found in object". This is due to not using the new

Hi Folks, The recent correspondence about "strange fisher.test result", and especially Peter Dalgaard's reply on Tue 03 April 2007 (which I want to investigate further) led me to take a close look at the code for binom.test(). I now have a query! The code for the two-sided case computes the p-value as follows: if (p == 0) (x == 0) else if (p == 1) (x == n)

Dear helpeRs, I seem to be a little bit confused on the result I am getting from the few codes below: > u=v=seq(0,1,length=30) > u [1] 0.00000000 0.03448276 0.06896552 0.10344828 0.13793103 0.17241379 [7] 0.20689655 0.24137931 0.27586207 0.31034483 0.34482759 0.37931034 [13] 0.41379310 0.44827586 0.48275862 0.51724138 0.55172414 0.58620690 [19] 0.62068966 0.65517241 0.68965517 0.72413793

_________________________________________________________________ http://redirect.gimas.net/?n=M0906xWLM2009 Neu: Messenger 2009! Hier kostenlos downloaden!

Can someone help interprete the error message below? i was trying to load the package "copula" from the R command prompt. > Error in loadNamespace(package, c(which.lib.loc, lib.loc), keep.source = keep.source) : in 'copula' methods specified for export, but none defined: show, summary, persp, contour Error: package/namespace load failed for 'copula' Thanks, Dominique

Dear R helper, does anyone have an idea on why R.2.4.0 draws the surface for the two command lines below and the next time it renders the error message below for exactly the same command lines: > norm.cop <- normalCopula(0.5) > persp(norm.cop, dcopula) Error in ceiling(length.out) : Non-numeric argument to mathematical function. I will appreciate any help from anyone thanks, Dominique

On Fri, 22 Feb 2019 at 17:27, Rowland Penny via samba <samba at lists.samba.org> wrote: > > [global] > > log level = 3 > > unix charset = UTF8 > > dos charset = ISO-8859-1 > > workgroup = ACME.DOM > > Sorry, but you cannot use a workgroup with a dot in it. > I tried to remove ".DOM" but then I can not start