パッケージRcppSimpleTensor(https://github.com/tlamadon/RcppSimpleTensor)のtensorFunctionを使用して、次のような3次元関数を計算します。
DGP <- function( x,y,z ) {
return((x + y - 5)^2 + (z-5)^2)
}
x、y、zの値の3つのベクトルを作成します。
xdata <- seq(0,10,length=10)
ydata <- seq(0,10,length=10)
zdata <- seq(0,5,length=15)
tensorFunctionを定義します。
library(RcppSimpleTensor)
RcppVals <- tensorFunction( R[i,j,k] ~ (X[i] + Y[j] - 5)^2 + (Z[k] - 5)^2 )
RcppArray <- RcppVals(xdata,ydata,zdata)
ただし、この関数を実行すると、「tmpfun(X、Y、Z、k、j、i)のエラー:単一の値が必要です」というエラーが発生します。
私はここで何が間違っているのですか?