在 Rcpp 中使用 gsl_ran_multinomial 時出現問題

[英]Problem when using gsl_ran_multinomial in Rcpp

我正在嘗試盡快生成多項式隨機變量。 我了解到gsl_ran_multinomial可能是一個不錯的選擇。 但是，我根據這篇帖子的回答嘗試使用： https://stackoverflow.com/a/23100665/21039115 ，結果總是報錯。 詳細來說，我的代碼是 結果是這樣的 這太荒謬了 ...

從新的 package 命名空間中的 Rcpp function 加載 R 包的 function

[英]Loading a R package's function from Rcpp function in a new package namespace

我正在為 R 開發 package，並希望在 Rcpp function 中加載 R package 的功能。 我知道 Rcpp::Environment env("package:package_of_interest")，但這僅在通過庫（"package_of_interest"）從 R 加載 ...

使用 Rcpp 循環進行日期迭代

[英]Date iteration with Rcpp loop

出於緊固目的，我試圖將 R 中的一個簡單的“for 循環”轉換為 Rcpp 循環。 我有一個名為“date_vector”的日期向量，它由 X 個相同的日期組成。 對於 i 的每次迭代，我將 1 分鍾添加到 date_vector 值。 R 'for 循環'（見下文）工作正常，但對於我非常大的數據 ...

有沒有一種簡單的方法可以在 Rcpp 中生成多維數組並將該數組導出到 R

[英]Is there a simple way to generate a multidimensional array in Rcpp and export that array to R

我正在尋找一種使用 Rcpp 構建 3D 數組並使該數組可在 R 中訪問的簡單方法。總體而言，我對 Rcpp 和 c++ 編碼仍然很陌生，所以任何幫助將不勝感激。 這是我的源代碼： 這個想法是我可以將 fit() 分配給 R（比如 D）中的一個元素，並且能夠從該數組中的位置調用。 我意識到我可以讓 ...

使用 Rcpp 從 R 運行 MPI C++ 程序

[英]Running MPI C++ program from R using Rcpp

我知道我可以使用帶有mpirun的 x 處理器運行 MPI C++ 程序： 是否可以使用 Rcpp 庫修改 MPI C++ 程序，以便可以使用 x 處理器從 R 調用它，而無需按照從 R 調用 MPI 現在運行 C 代碼中的建議求助於 Rmpi 庫？ 為了說明的目的，我修改了從 R 調用的 MPI ...

為 std::array 專門化 Rcpp::as()

[英]Specialising Rcpp::as() for std::array

我的一個項目將 C++11 std::array用於固定大小的數組類型，因此我試圖專門Rcpp::as()以方便從SEXP到這些類型的轉換。 由於這需要部分專業化，因此我遵循了 Rcpp-extending 中概述的Exporter路線：#include <RcppCommon.h> ...

將 outer() 從 base R 翻譯成 RcppArmadillo

[英]Translate outer() from base R to RcppArmadillo

有什么方法可以有效地將用於將兩個向量相乘的 outer() 函數從 R base 轉換為 RcppArmadillo？ 我試圖這樣做，但它根本沒有效率。 舉個例子： 使用 R base 結果慢了 2 倍。 我確信這可以通過矩陣乘法來完成，知道怎么做嗎？ ...

使用 Rcpp 從 std::uniform_int_distribution 采樣時出現“來自 C 堆棧溢出的段錯誤”

[英]Getting "segfault from C stack overflow" when using Rcpp to sample from std::uniform_int_distribution

我編寫了這個像 R 的sample.int(..., replace =FALSE)函數一樣工作的 C++ 小函數。 本質上，它從均勻分布的整數中提取並將結果寫入一個集合，直到該集合的大小為size 。 也許我在這里遺漏了一些東西，但我發現以下行為很奇怪。 這是一個代表： 使用Rcpp執行此操作會 ...

使用 Rcpp 遍歷大型多維數組

[英]Looping Through Large, Multidimensional Array Using Rcpp

我正在嘗試創建涉及循環遍歷大型多維數組（例如：維度 = 20 x 1000 x 60）的模型，這些模型運行速度非常慢，就像我在 R 中編寫它們一樣。我下載了 Rcpp 並一直在嘗試實現這樣的模型，因為 C++ 可以很好地處理循環。 通常，我會在 R 中編寫這樣的函數： 正如預期的那樣，函數產生： 這 ...

使用 VCM 提高我的 R function 性能

[英]Improving my R function performance with VCM

我正在使用變系數模型運行模擬，但有一些調整。 沒有 R package 可以滿足我的需求。 我的代碼運行速度不夠快。 我期待讓vcm function 運行得更快 到目前為止我所做的是剖析代碼並查看緩慢的部分使用sapply和lapply而不是for循環，但沒有顯着差異如何使用代碼？ 這是一個使用代 ...

rcpp function [實際上不是] package 內部比外部慢很多

[英]Rcpp function [actually not] much slower inside of package than outside

我最近在 Rcpp 中寫了一個計算密集型 function。 現在，我想將這段代碼移植到 R package。但是，我注意到當在 R package 中運行時，代碼要慢很多（~100 倍）。 我已經在這里讀到，這可能與 function 的調用方式有關。 但是，這不是一次性成本。 相反，它與 Rcp ...

使用特征值反轉稀疏矩陣

[英]Inverting a sparse matrix using eigen

我正在嘗試使用稀疏求解器作為 SimplicialLLT 來求逆對稱正定矩陣並返回它。 我從 R 使用 Rcpp 連接 R 和 cpp 得到一個矩陣，我將這個矩陣作為 function cpp_sparse_solver 的參數，使用 sparseView() 將其轉換為 SparseMatrix， ...

在 Rcpp 中將 arma:vec 轉換為 NumericVector

[英]Convert arma:vec to NumericVector in Rcpp

我嘗試在 RcppArmadillo 中使用分位數 function，但是，在教程中，分位數 function 的 output 是 arma::vec（output 的 class in 8814664"82"arraytrix""） arma 中的分位數 function 我想將它轉換為 N ...

在 Rcpp 中檢查 boolean/LogicalVector

[英]Check boolean/LogicalVector in Rcpp

我寫這個 Rcpp function： 這個function的目標是： 1）當沒有輸入時，返回“null”； 2）當輸入TRUE時，返回“true”； 3) 當輸入FALSE時，返回“FALSE”。 然而，當我嘗試這個 function 時，我發現輸入 FALSE 仍然給我“TRUE”。 截屏 ...

這在 Rcpp 中是否正確？

[英]Is this correct in Rcpp?

我想比較每一列，並在計算后返回所有結果。 我嘗試編寫代碼，但結果並不合理。 因為如果矩陣中有 5 列，結果的數量將是 5*4/2=10 而不是 5。我認為問題是代碼中的m 。 我不知道這是否正確。 謝謝。 ...

使用 Rcpp，如何將 RawVector 來回轉換為字符串？

[英]Using Rcpp, how to convert RawVector back/forth to a string?

在Rccp-dev上有一個討論，討論如何從字符串結果轉換為 RawVector（參見https://lists.r-forge.r-project.org/pipermail/rcpp-devel/2010-October/001179.ZFC35FDC70D5FC69D269883A822AZ5 ...

使用 Rcpp 和 Kronecker 產品在 R 中計算 nnet::multinom 多項回歸的 Hessian / Fisher 信息矩陣的更快方法

[英]Faster way to calculate the Hessian / Fisher Information Matrix of a nnet::multinom multinomial regression in R using Rcpp & Kronecker products

似乎對於較大的nnet::multinom多項式回歸模型（具有幾千個系數），計算 Hessian（負對數似然的二階導數矩陣，也稱為觀察到的 Fisher 信息矩陣）變得非常慢，然后阻止我計算方差-協方差矩陣並允許我計算 model 預測的置信區間。 It seems the culprit is ...

在 rcpp 中處理浮點數：R 中的值 -0 使用 Rcpp 在 C++ 中變為 0

[英]handling of floating points in rcpp: a value of -0 in R becomes 0 in C++ using Rcpp

使用 Rcpp 時，R 中值為 -0 的參數在 C++ 世界中變為 0。 這是一個演示這一點的最小示例。 這是 output： 簡而言之，-0 和 0 的差異消失了。 這是一個錯誤嗎？ 有什么解決辦法嗎？ ...

如何檢查 IntegerVector 是否在 Rcpp 中包含 NA 值？

[英]How do I check whether an IntegerVector contains NA values in Rcpp?

我希望檢查提供給 C++ function 的 Rcpp IntegerVector 是否不包含NA值。 在另一個答案之后，我寫了以下內容：IntegerMatrix stop_if_na(const IntegerVector I) { if (Rcpp::is_true(Rcpp::any ...

制作 R package 時使用 rcpp armadillo 時出錯

[英]Error using rcpp armadillo when making R package

我正在嘗試在我的 R 包中使用 C++ 代碼。 我創建了一個測試 package 使用 然后我使用 C++ 然后我創建了R/test-package.R並粘貼在 在sum_func.cpp然后我把 然后我創建了R/sum__func.R並放入 但是當我運行devtools::load_all("." ...