Message-ID: <52657A00-C3D9-4E96-AF81-1F06F49D4402@comcast.net>
Date: 2013-11-21T04:39:42Z
From: David Winsemius
Subject: integrate
In-Reply-To: <CA+1Z_YvpWgOyOUmCrz-gu3MEcRCVh+Mskpm7+az9GxtzREU01w@mail.gmail.com>
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"?
--
David Winsemius
Alameda, CA, USA