我试图通过重复采样（蒙特卡洛模拟）来估计参数。每次采样时，我都必须对data.table对象运行特定操作。这是一个例子：

visits.site <- site.visits.map[site.name %in% sample.sites, 
                              .(length(unique(visits))),
                              by = week]

由于采样的数量经常达到数百万，这需要很多时间。我使用该parApply函数使循环在 CPU 的多个内核上并行执行。这明显提高了性能。（我在 8 个内核上运行它）。注意：这是一个令人尴尬的并行问题。

我想通过在 GPU 上运行它来进一步提高性能。我看过gpuR但我不确定它是否具有类似于parApplyGPU 的功能。R中是否有任何具有此功能的软件包？

如果还有其他方法可以提高性能，请告诉我。