运行 Rcpp 时 R Studio 崩溃

Question

我知道 RO 是一个非调试平台，但我是 Rcpp 的新手，所以希望有人能帮助我为什么会发生此错误，因为我不确定它是否与 Rstudio 或 Rccp 函数有关。

这里有一个示例数据。我放置原始数据的原因是为了确保人们可以重现我遇到的错误。

https://www.dropbox.com/s/s9xgu5h0qfojy9q/test.csv?dl=0

test <- read.csv("test.csv")
test <- data.table(test)

我使用 Rcpp 从我的一篇帖子中获得的功能如下：

Rcpp::cppFunction('
              List WRSI(NumericVector daily_tawc,
              NumericVector TAWC,
              NumericVector Precp,
              NumericVector swc,
              NumericVector PETc) {

              int n = daily_tawc.length();
              NumericVector SW(n);
              NumericVector PAW(n);
              NumericVector AETc(n);

              double SW_ini_glob = daily_tawc[0];
              double SW_max_glob = TAWC[0];

              SW[0] = SW_ini_glob;

              for (int i = 1; i <= n; i++) {

              PAW[i] = SW[i-1] + Precp[i];

              if (PAW[i] >= swc[i]) {
              AETc[i] = PETc[i];
              } else {
              AETc[i] = PAW[i]/swc[i]*PETc[i];
              }

              if (AETc[i] > PAW[i]) {
              AETc[i] = PAW[i];
              }

              SW[i] = SW[i-1] + Precp[i] - AETc[i];

              if(SW[i] > SW_max_glob) {
              SW[i] = SW_max_glob;
              }

              if(SW[i] < 0) {
              SW[i] = 0;
              }
              }
              return Rcpp::List::create(Rcpp::Named("SW") = SW,
              Rcpp::Named("PAW") = PAW,
              Rcpp::Named("AETc") = AETc);
              }')

当我为测试数据运行上述函数时，

  test[, c("SW","PAW","AETc") := WRSI(daily_tawc,
                                       TAWC,
                                       Precp,
                                       swc,
                                       PETc), keyby = .(state.id,meso.id,micro.id,mun.id,year)]

每次 Rstudio 崩溃时。我想知道是否有人可以告诉我这是怎么回事？ Rstudio 在运行上述功能时崩溃的可能原因是什么。我不确定函数中是否有什么东西让我的 Rstudio 崩溃。

谢谢

Answer 1

无需详述或搜索重复项：

for (int i = 1; i <= n; i++)

这看起来错了。应该是这样的

for (int i = 0; i < n; i++)