r - モデルに確率論を追加するのに問題がある

Question

この質問が既に尋ねられている場合は事前にお詫び申し上げますが、オンラインで役立つ情報が見つからないようです!

ロジスティック成長モデルを含む、生成した関数に確率論を追加しようとしています。

これは私の機能です：

    ricker <- function(K, R, B0, t, E, Q)
    {
  # create empty matrix with only biomass of first year; this is where all         other values will be added
  B <- matrix(nrow = t, ncol = length(Q))
  B[1,1:length(Q)] <- B0

  for (i in 1:t) {
    for (j in 1:length(Q)) {
      B[i + 1,j]  <- 
{
  B[i,j] + R * B[i,j] * (1 - B[i,j]/K) - Q[j] * E * B[i,j]

  }
}
}
return(B)
}

上で作成した行列には、34 行と 34 列があります。以下は、データに確率論を追加するために使用しているコードです。

# create variables/parameters
k <- 420000 
b0 <- as.matrix(3363) 
T <- 34 # number of years in data 
e <- fp 
r=1.2

pb <- txtProgressBar(style=3)

for (loop in 1:10) {

  B <- ricker(K=k, R=r, B0=b0, t=T, E=0.5, Q=fp)

  B[loop] <- B[loop] - rnorm(1,0,0.25)

  if (B[loop] <0) {

   B [loop] <- 0
  }
  setTxtProgressBar(pb,value=l/T)
}
#'
matplot(B,type="l", col=1)

私はこのようなグラフになりたい：

ただし、出力は最初の出力と何ら変わりはなく、何が欠けているのかわかりません。この問題に関するヘルプは大歓迎です!!

ありがとう！！