我遇到了這個奇怪的段錯誤並且沒有線索如何解決它。我正在運行一些馬爾可夫鏈蒙特卡洛算法(一種近似分佈的順序算法)。我並行化該算法的每次迭代。所以它就像 for (iter in 1:T){
res[iter] = mclapply(fun)
}
現在奇怪的是,當我的數據集的大小相對適中時,該算法可以運行沒有問題。然後我增加數據集的大小(80,000個觀察值,不是超大的),該算法適用於
我使用的是doSNOW包,更具體地說是parLapply函數對大柵格數據集(操作系統:Windows x64)列表執行重分類(以及後續的其他操作)。 代碼看起來有點像這樣簡約的例子: library(raster)
library(doSNOW)
#create list containing test rasters
x <- raster(ncol=10980,nrow=10900)
我試圖做到這一點: nmf.sub <- function(n){
sub.data.matrix <- data.matrix[, (index[n, ])] ## the index is a permutation of the original matrix at a 0.8 resampling proportion (doesn't really matter)
temp.resu