0

forループ環境で関数「cop.theta」を使用して、いくつかの計算から行列を作成するためのこのコードがあります

Mat.corr <- matrix(0,6,5,byrow=F)
for (i in 1:6){
Mat.corr[i,]=cop.theta(index,EXPR,SURV=survp[,i])
}

上記のコードが生成したものと同様の結果を得るために、doParallel パッケージで「foreach」を使用して R コードを作成しました。私のコードは次のとおりです

library(doParallel)
cl <- makeCluster(3) 
registerDoParallel(cl)
getDoParWorkers()
clusterExport(cl, list("QT","EXPR","cop.theta.i"))
clusterEvalQ(cl, library(copula))

foreach(i=1:6,.combine=matrix(0,6,5,byrow=F) %dopar%
Mat.corr[i,]=cop.theta(index,EXPR,QT=survp[,i])

しかし、私はこのエラーが発生しています

Error: unexpected '=' in "foreach(i=1:6,.combine=matrix(0,6,5,byrow=F)
%dopar%  Mat.corr[i,]="

私はどこで間違っていますか？

4

2 に答える 2

1

byrow = Fの横に ")" がありません。

foreach(i=1:6,.combine=matrix(0,6,5,byrow=F)) %dopar%
  Mat.corr[i,]=cop.theta(index,EXPR,QT=survp[,i])

このエラーを修正する必要があります。ただし、ループを並行して実行している限り、ループ内では操作の結果の値がこの特定の行に書き込まれないことに注意してください。これは single-core でのみmat.corr[i,]=...可能foreachです。つまり、内側またはループが終了した後にマージ手順を処理する必要があります。これが私の要点を明確にするためのコードスニペットです。cop.thetaforeach.combineforeach()

mat <- matrix(nrow = 3, ncol = 3)

### multi-core -----
foreach(i = 1:nrow(mat)) %dopar% {
  mat[i, ] <- matrix(rep(i, ncol(mat)), nrow = 1)
}

mat
     [,1] [,2] [,3]
[1,]   NA   NA   NA
[2,]   NA   NA   NA
[3,]   NA   NA   NA


### single core -----
foreach(i = 1:nrow(mat)) %do% {
  mat[i, ] <- matrix(rep(i, ncol(mat)), nrow = 1)
}

mat
     [,1] [,2] [,3]
[1,]    1    1    1
[2,]    2    2    2
[3,]    3    3    3
于 2016-04-06T10:07:06.437 に答える
0

コメントありがとうございます。以下のコードは実際に私の質問に答えています。

foreach(i=1:6,.combine=rbind) %dopar%  cop.theta(index,EXPR,SURV=survp[,i])
于 2016-04-06T16:07:03.437 に答える