Re: [R] Integration in R

On Jan 8, 2013, at 1:31 PM, David Winsemius wrote: On Jan 8, 2013, at 1:07 PM, Berend Hasselman wrote: On 08-01-2013, at 22:00, Berend Hasselman wrote: …... David implemented the condition by multiplying by x[1]in a numeric context is 0 when x[1]=x[2]. OOPS!! Reverse the 0 and the 1 in

Re: [R] Integration in R

On Jan 8, 2013, at 1:07 PM, Berend Hasselman wrote: On 08-01-2013, at 22:00, Berend Hasselman wrote: …... David implemented the condition by multiplying by x[1]in a numeric context is 0 when x[1]=x[2]. OOPS!! Reverse the 0 and the 1 in that sentence (TRUE becomes 1 and FALSE becomes 0)

Re: [R] Integration in R

On 08-01-2013, at 22:00, Berend Hasselman wrote: > …... > David implemented the condition by multiplying by x[1] numeric context is 0 when x[1]=x[2]. OOPS!! Reverse the 0 and the 1 in that sentence (TRUE becomes 1 and FALSE becomes 0) Berend __ R-h

Re: [R] Integration in R

On 08-01-2013, at 19:51, Naser Jamil wrote: > Thanks. But then how to implement condition like 0 happy to know that. > David implemented the condition by multiplying by x[1]=x[2]. That is what your requirement does. The condition 0x2 which is what your inequality implies. Berend > On 8 Jan

Re: [R] Integration in R

Thanks. But then how to implement condition like 0 wrote: > Please reply on list. > > > On Jan 8, 2013, at 10:27 AM, Naser Jamil wrote: > > Hi David, >> x[2] is the second variable, x2. It comes from the condition 0> > > No, it doesn't come from those conditions. It is being grabbed from some > "

Re: [R] Integration in R

On Jan 8, 2013, at 10:51 AM, Naser Jamil wrote: Thanks. But then how to implement condition like 0be happy to know that. Multiply the function by the conditional expression: > f<-function(x) { 2/3 * (x[1] + x[2] )*(x[1] < x[2]) } > adaptIntegrate(f, lowerLimit = c(0, 0), upperLimit = c(7,7))

Re: [R] Integration in R

Please reply on list. On Jan 8, 2013, at 10:27 AM, Naser Jamil wrote: Hi David, x[2] is the second variable, x2. It comes from the condition 0 No, it doesn't come from those conditions. It is being grabbed from some "x"-named object that exists in your workspace. If your limits were 7 in

Re: [R] Integration in R

On Jan 8, 2013, at 9:43 AM, Naser Jamil wrote: Hi R-users. I'm having difficulty with an integration in R via the package "cubature". I'm putting it with a simple example here. I wish to integrate a function like: f(x1,x2)=2/3*(x1+x2) in the interval 0by hand and got 114.33, but the follow

Re: [R] Integration in R

Send us a reproducible R code that shows what you actually tried. Ravi Varadhan, Ph.D. Assistant Professor The Center on Aging and Health Division of Geriatric Medicine & Gerontology Johns Hopkins University rvarad...@jhmi.edu 410-502-2619 [[alternative HTML ve

Re: [R] Integration in R

wdunlap tibco.com > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf > Of Rolf Turner > Sent: Wednesday, November 21, 2012 12:57 PM > To: Rehena Sultana > Cc: r-help@r-project.org > Subject: Re: [R] Integration in R &g

Re: [R] Integration in R

On 21/11/12 22:26, Rehena Sultana wrote: Dear R - Experts, I am trying to integrate lognormal distribution (mu = -0.3 and sigma2 = 0.00041.. ) based on the some hypothetical data. But I am getting 0 as the result. I have checked that my R-code is correct as code is giving me result for some o

Re: [R] Integration in R

On 02-10-2012, at 20:50, Dereje Bacha wrote: > Hi > > I am facing a problem of restricting an intercept of systems of equations. > Y1=f(X1,X2,X3) > Y2=(X1,X2,X4) > > I want to restrict an intercept of equation 2 equal to coefficient of X2 of > equation 1. > Please do not hijack a thread

Re: [R] Integration in R

Hi I am facing a problem of restricting an intercept of systems of equations.   Y1=f(X1,X2,X3) Y2=(X1,X2,X4)   I want to restrict an intercept of equation 2 equal to coefficient of X2 of equation 1.  Please help Dereje   From: Naser Jamil To: r-help@r-project

Re: [R] Integration in R

ginal Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf > Of Berend Hasselman > Sent: Tuesday, October 02, 2012 11:19 AM > To: Rui Barradas > Cc: r-help@r-project.org; Naser Jamil > Subject: Re: [R] Integration in R > > > O

Re: [R] Integration in R

Hello, Em 02-10-2012 19:18, Berend Hasselman escreveu: On 02-10-2012, at 20:01, Rui Barradas wrote: Hello, Yes, it's possible to remove the loop. Since the loop is used to compute a running product and all we want is the final result, use the vectorized behavior of R and a final ?prod(). S

Re: [R] Integration in R

On 02-10-2012, at 20:01, Rui Barradas wrote: > Hello, > > Yes, it's possible to remove the loop. Since the loop is used to compute a > running product and all we want is the final result, use the vectorized > behavior of R and a final ?prod(). > Seedup: another 2x. And 4x2 == 8 == 1 [decimal]

Re: [R] Integration in R

Hello, Yes, it's possible to remove the loop. Since the loop is used to compute a running product and all we want is the final result, use the vectorized behavior of R and a final ?prod(). Seedup: another 2x. And 4x2 == 8 == 1 [decimal] order of magnitude. lf2 <-function (x) { v<-1 x1

Re: [R] Integration in R

On 02-10-2012, at 17:23, Naser Jamil wrote: > 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 wa

Re: [R] Integration in R

You are missing basic algebra skills! You had: myfunc<- function(x) {0.25*(9*x^4 + 6*x^2 + 1)} This should be: myfunc<- function(x) {0.25*(9*x^4 - 6*x^2 + 1)} Ravi. --- Ravi Varadhan, Ph.D. Assistant Professor, Division of Geriatric Medici