多个索引上的Rcpp和

我想计算一个函数的形式：

$m_{jl}(x) = x + \gamma[j]*zeta[j,l] + sum_{k \neq j} zeta[j,k]$

---

使用

Rcpp

。我的问题是关于sum{k\neq j}zeta[j，k]$。。我希望能够做些事情

zeta[j，-j]

。可能吗？我试过

zeta（j，-）-zeta（j，j）

，但是

Rcpp

不喜欢

zeta（j，-）

，你可以使用sugar函数：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
double myFun(NumericMatrix x) {
   return sum(x)-sum(diag(x));
}

#包括
使用名称空间Rcpp；
//[[Rcpp:：导出]]
双myFun（数字矩阵x）{
返回和（x）-和（diag（x））；
}

R:

---

A-1，因为缺乏研究。此外，当您正确拼写包名时，包作者更有可能做出响应。。。
A <- matrix(1:9,3)

sum(A)-sum(diag(A))
#[1] 26

myFun(A)
#[1] 26