神経筋接合シュミレーション−３

１回の単位時間辺りに、一個、二個、三個の神経細胞がやってくると仮定してみる。

N<- 10000

n<-100

time<-1000

kinnniku<-matrix(0 , time , N)

sinnkei<-matrix(0 , time , n)

p<- sample(1:n , N, replace = T)

p

kinnniku[1, ] <-p

sinnkei[1, ] <- tabulate(kinnniku[1, ])

for( j in 1: (time - 1) ){

p1<- sample(1:n, N, replace = T)

#一対一の戦いでなく、集合の戦いにする

p2<- sample(1:n, N, replace = T)

p3<- sample(1:n, N, replace = T)

    p <- sinnkei[ j , ]

       for( i in 1 : N ){

      a<-kinnniku[ j , i ]

      b<-p1[ i ]

      c<-p2[ i ]

      d<- p3[ i ]

          num.a <-  p[ a ] 

      num.b <-  p[ b ]

      num.c <-  p[ c ]

      num.d <-  p[ d ]

      prob.i <- c( num.a , num.b ,num.c ,num.d )/(num.a + num.b +num.c+num.d )

      new.i <- sample(c( a ,b ,c ,d ), (num.a + num.b + num.c +num.d ) , prob = prob.i , replace=T)

      kinnniku[ j+1, i ]<- new.i [1] 

}

sinnkei[ j+1, ] <-tabulate( kinnniku[j+1, ] , nbins=n)

}

 plot(sort(sinnkei[ time, ]))