Question

我的目标是获取c_y1下的数字区域，其定义如下：

require(expm)

    c_y1 <- function(u){
  pi_initial[1]*(expm(T_initial*u))*1*(expm(T_initial*(Y[v,1]-u)*t_initial))} 
Z1V.YV <- list()
for (v in 1:nrow(Y)){
  Z1V.YV[[v]] <- {(area(c_y1, 0, v))/(expm(T_initial*Y[v,1])*t_initial)}
}

其中expm()是其内容的矩阵幂。此外，T_initial是一个矩阵，t_initial是一个向量，pi_initial[1]是一个实数，Y是一个数据帧，使得Y[,1]是一个向量nrow(Y)维实数空间。因此，目标是在此矩阵的每个元素上进行集成。

我收到以下形式的警告：

> warnings()
Warning messages:
1: In if (abs(a1 - a2) < eps) return(a2) :
  the condition has length > 1 and only the first element will be used

以至于存在nrow(Y)-1条与此完全相同的警告。什么地方出了错？有没有办法在所有v上逐个区域应用面积函数？

在R中的矩阵上积分

0 个答案: