Message-ID: <0903A0CD-EF4D-49EE-A0F6-5F05431079C0@comcast.net>
Date: 2013-11-21T04:52:15Z
From: David Winsemius
Subject: integrate
In-Reply-To: <52657A00-C3D9-4E96-AF81-1F06F49D4402@comcast.net>
On Nov 20, 2013, at 8:39 PM, David Winsemius wrote:
>
> On Nov 20, 2013, at 7:12 PM, dan wang wrote:
>
>> Hi all,
>>
>> Can anyone help me with below integrate function?
>>
>> Basically, I want to calculate the integral of the sum of two kernel
>> density functions.
>> But the error shows that:
>> In x - a :
>> longer object length is not a multiple of shorter object length
>>
>> y1 = rnorm(10)
>> y2 = rnorm(10)
>> fhat <- function(x,a){h=density(a)$bw;sum(dnorm((x-a)/h)/h/length(a))}
>>
>> integrand1 <- function(p) {
>> fhat(p,y1)+fhat(p,y2)
>> }
>> integrate(integrand1,lower = -Inf, upper = Inf)
>
> So, .... where is "a"?
Ignore that person. Try this:
Vintegrand <- Vectorize(integrand1)
integrate( Vintegrand, lower = -Inf, upper = Inf)
#----------------
2 with absolute error < 2.4e-05
David Winsemius
Alameda, CA, USA