Hans W Borchers wrote :
> First define a function from those points:
>
> fx <- approxfun(x, f_x)
> fy <- approxfun(y, f_y)
> f <- function(x) abs(fx(x)-fy(x))
>
> and now you can apply integrate() or trapz():
>
> xx <- sort(c(x, y))
> yy <- f(xx)
> trapz(xx, yy)
>
> tra
Le 25.01.2011 15:23, Rmh a écrit :
> g <- function(x) abs(f1(x)-f2(x))
>
> now you have one function and you can integrate it.
Thank you Rich.
Unfortunately I have no f1 and f2 functions, only a set of observed
points on two lines - and no idea about the underlying distribution to
create a funct
Xavier Robin unige.ch> writes:
> Hello,
>
> I need to integrate the absolute difference between two lines measured
> on different points.
>
> # For example :
> x <- seq(0, 1, 1/100)
> f_x <- runif(101) + x
> y <- seq(0, 1, 1/23)
> f_y <- runif(24) + (1 - y)
>
> plot(x, f_x, type="l")
> lines(y
g <- function(x) abs(f1(x)-f2(x))
now you have one function and you can integrate it.
Rich
Sent from my iPhone
On Jan 25, 2011, at 7:32, Xavier Robin wrote:
> Hello,
>
> I need to integrate the absolute difference between two lines measured
> on different points.
>
> # For example :
> x <-
Hello,
I need to integrate the absolute difference between two lines measured
on different points.
# For example :
x <- seq(0, 1, 1/100)
f_x <- runif(101) + x
y <- seq(0, 1, 1/23)
f_y <- runif(24) + (1 - y)
plot(x, f_x, type="l")
lines(y, f_y)
Then I would like to compute Integral( | f_x - f_y
5 matches
Mail list logo