R：在多边形

Question

我想整合多边形，我已经尝试使用mvtnorm包中的polyCub.SV（）函数：

#defining the covariate function:

elev.fn<-function(s){
  s<- matrix(s,ncol=2)
  (s[,1]+s[,2]-100)/40.8
}


mu <- function(s,beta0,beta1) exp(beta0+beta1*elev.fn(s=s))


polyCub.SV(Mypolygon,f=mu,beta0=beta0,beta1=beta1,nGQ=2,rotation=FALSE)

结果是NAN

我做错了什么？