我有一个功能，我想采取它的双重积分。问题是函数内部有一些向量。我不知道如何对函数进行向量化，然后将每个函数整合到非向量的变量上。
beta=c(1,2,3)
sigmay=0.1
sigmax1=0.1
sigmay1=0.1
# each of stt, sezz and smm is an identified n element vector
func1 = function(smtheo, sttheo) {
Xtheo=c(1, log(smtheo), log(sezz))
mu <- beta %*% Xtheo
ypri=dnorm(log(sttheo), mu, sigmay)
joipri=dnorm(log(stt), log(sttheo), sigmay1)*dnorm(log(smm),log(smtheo), sigmax1)
return(ypri*joipri)}
int=integral2(func1(smtheo, sttheo), a1, a2, b1, b2)
实际上我有n个func1，每个都应该在smtheo和stttheo限制上集成。有谁知道吗？