similar to: How to calculate a double integral ...?

Dear R-users, I'm wondering how to calculate this double integral in R: int_a^b int_c^y g(x, y) dx dy where g(x,y) = exp(- alpha (y - x)) * b Thanks for answering! Cheers, Alui [[alternative HTML version deleted]]

Dear all, I found that the plotmath seems won't work on xlab in plot (and other plot functions including matplot, contour, ...): > plot(1:10, 1:10, xlab=expression(alpha)) the plot shows nothing for x-label. My R-version is R.1.1.0 > version _ platform Windows arch x86 os Win32 system x86, Win32 status major 1

Hello, I''m trying to compile Xen unstable on NetBSD, I''ve aplied the patches from http://cvsweb.netbsd.org/bsdweb.cgi/pkgsrc/sysutils/xentools41/patches/?only_with_tag=MAIN and copied the blk files from http://cvsweb.netbsd.org/bsdweb.cgi/pkgsrc/sysutils/xentools41/files/?only_with_tag=MAIN (I don''t know which of these patches have been applied to xen unstable, so

In rw1050, I found that > pweibull(3:10, 2) [1] 0.9998766 0.9999999 1.0000000 1.0000000 NaN NaN [7] NaN NaN Warning message: NaNs produced in: pweibull(q, shape, scale, lower.tail, log.p) more surprisingly, > pweibull(3:10, 2.1) [1] 0.9999566 1.0000000 1.0000000 -Inf NaN NaN [7] NaN NaN Warning message: NaNs produced in: pweibull(q,

Dear R and ESS users, I got the following error when trying to load a source file (smart.R in directory E:\R\) into XEmacs through the menu iESS and "Load source file": Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file `E:Rsmart.R' It seems the slash "/" or "\" being missing, but I don't know

Dear R-friends, For the following data: > xy x y i 1 731 0.313 2 2 739 0.340 2 3 790 0.373 2 4 855 0.451 2 5 980 0.608 2 6 575 0.156 1 7 608 0.207 1 8 630 0.249 1 9 670 0.332 1 10 838 0.377 1 11 964 0.466 1 > coplot(y ~ x|i, data=xy) coplot gives 3 panels, rather than 2, namely one for i=1 and one for i=2. Futhermore, when I extand data fram xy to have i=3 as follows:

Dear R-Friends, Sorry to bother everyone with my earlier question: "Do we have similar R function to work like Minitab's INDICATOR?" I found the way to make things work: e.g., > x <- c(2,2,5,3,6,5,NA) > model.matrix(~ factor(x) - 1) factor(x)2 factor(x)3 factor(x)5 factor(x)6 1 1 0 0 0 2 1 0 0 0

Hello all, I'm hoping to find a way to evaluate the following sort of integral in R. \int_a^b \int_{g(y)}^Inf f(x,y) dx dy. The integral has no closed form and so must be evaluated numerically. The "adapt" package provides for multidimensional integration but does not appear to allow the limits of integration to be a function. I need to evaluate a number of integrals of this

Dear R gurus, To start, let me confess to not being an experienced programmer, although I have used R fairly extensively in my work as a graduate student in statistics. I wish to find the root of a function of two variables that is defined by an integral which must be evaluated numerically. So the problem I want to solve is of the form: Find k such that f(k)=0, where f(y) = int_a^b g(x,y)

Hi, I read the 2000 thread on a MetaPost device in R. If I understand correctly, the main problem with the concept is that R wants the device driver to give back information on the size of strings/labels. To the bet of my knowledge, MetaPost _does_ make it possible to measure the bounding box of text (see section 7.3: Measuring text in the MetaPost manual). For example, one could get the size of

As I have previously asked, in response to a similar question: Is this a homework problem? cheers, Rolf Turner rolf at math.unb.ca

Dear list, [cross-posting from Stack Overflow where this question has remained unanswered for two weeks] I'd like to perform a numerical integration in one dimension, I = int_a^b f(x) dx where the integrand f: x in IR -> f(x) in IR^p is vector-valued. integrate() only allows scalar integrands, thus I would need to call it many (p=200 typically) times, which sounds suboptimal. The

Hi, I am wondering if anybody ever come across any implementation of product integral in R? As far as I googled, I haven't come across any package. Is there any? Thank you. http://en.wikipedia.org/wiki/Product_integral Regards, Robert

Dear R-friends, When I selected "Save as postscript" via mouse left-button on the active "R Graphics" window, those colored lines I made in the plot changed in the saved postscript file, e.g., color blue becomes green, red becomes blue, ... I have not changed any setting/options. What's going on? My version of R is > version _ platform i386-pc-mingw32 arch

Hi all, I realize this is OT. I just got a Cisco AS5300, and I need to configure it like such: Asterisk -----(H323/SIP)------> Cisco ----- (E1/PRI)------->Telco So calls originate from the Asterisk side (registered users on SIP or just ZAP phones), and they go out H323 or SIP to Cisco, where they go out PRI. I have the Asterisk side sorted :) (either H323 or SIP), I need help in the

Dear R-friends, As the binary windows version of R being distributed, what is the correct way to install contributed packages without using winst.exe as before? The answer in "R for Windows FAQ" says: "You can install pre-compiled packages either from CRAN or from a local .zip file by using install.packages: see its help page." But the information available in

Dear R users, This is a KDE beginner's question. I have this distribution: > length(cap) [1] 200 > summary(cap) Min. 1st Qu. Median Mean 3rd Qu. Max. 459.9 802.3 991.6 1066.0 1242.0 2382.0 I need to compute the sum of the values times their probability of occurence. The graph is fine, den <- density(cap, from=min(cap), to=max(cap), give.Rkern=F)

Spencer Thanks for your thoughts on this. I did a bit of work and did end up with a method (more a trick), but it did work. I am certain there are better ways to do this, but here is how I resolved the issue. The integral I need to evaluate is \begin{equation} \frac{\int_c^{\infty} p(x|\theta)f(\theta)d\theta} {\int_{-\infty}^{\infty} p(x|\theta)f(\theta)d\theta} \end{equation} Where

Hi.. i was trying to integrate the indicator funtion but had problems when limits where negative or equal to the indicator condition my function is ________________________ fun1<-function(x){ as.numeric(x>=2) } _________________________ which should be Ind(x>=2)*x seems to work for the following two cases

I thought I was imaging it but... I have a friend (honestly) that likes to clear his inbox every month by copying older into a folder (e.g. June-2009). Every month he complains that all the mail in the new folder has disappeared. This has happened three times now, and he's sort-of-right. I now believe it's a bug rather than him. The setup is Dovecot 1.1.3 on FreeBSD 7.0 and the client is