Как интегрировать только переменную в многомерную функцию в R?

Question

У меня есть следующая функция:

f <- function(x,y,z,w) {dpois(x, lambda = w*cos(y-z) ) }

Я хочу построить функцию g (y, z, w), определенную интегралом от f по x. Я могу сделать это, когда у меня есть две переменные, например, такой код:

margin <- function(y) { sapply(y, function(y) { integrate(function(x) f(x,y), llimx, ulimx)$value }) }

Но я не знаю, как это сделать для более чем двух переменных.

Answer 1

Может быть, вы можете попробовать код ниже

f <- Vectorize(function(x,y,z,w) dpois(x, w*cos(y-z)),"x") 
g <- function(y,z,w) integrate(f, lower = 0, upper =Inf, y, z, w)