如何在R中的多元函數中僅積分一個變量？

Question

我有下一個功能：

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

我想構造一個由 f 對 x 的積分定義的函數 g(y,z,w)。 當我有兩個變量時我可以這樣做，例如代碼是這樣的

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)