如何模拟状态转换矩阵？

Question

假设我想模拟5乘5状态转移矩阵，即如果i和j是状态，则i = j = 0（1）5。对于给定的＆＃34;我＆＃34;我想使用下面给出的关系生成大小为10的随机样本（j）（最终构成的过渡数将变为N = 50）

j=ln(1-U) -1; i=0
j=ln(1-U)+i-2;  i>=1

其中U~Uniform（0,1）。

我尝试将我的R代码编写为：

 nchain <- 100 ## no of simulation##

 i <- 0:4

 U <- runif(50,0,1)

 j <- rep(0,nchain)

 for(k in 1:nchain){
  if (i=0) {
      j=log(1-U)-1
  } else {
    j=log(1-U)+i-2
  }
}

但我无法做到。 如何模拟100种这样的矩阵？