Моей целью было получить числовую область под 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?