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++)

当Rcpp运行时，R studio崩溃了

1 个答案: