Использование функции rbinom в R в матрице - PullRequest
Новая
       6

Использование функции rbinom в R в матрице

0 голосов
/ 27 ноября 2018

Предисловие: я довольно новичок в использовании R, я использовал SAS всю свою взрослую жизнь и не привык работать с матрицами.

В настоящее время я работаю над проектом для класса эволюционной биологии.это требует запуска функции rbinom через вложенный цикл над матрицей.Вероятность в первой строке установлена ​​равной 0,1, но затем значение в последующих строках должно использовать вероятность из предыдущей строки.Я не могу понять, как ссылаться на значение в предыдущей строке.Мой код ниже, если кто-нибудь знает синтаксис этого, я был бы очень признателен!В настоящее время у меня установлено значение i-1, но я знаю, что это неправильно. 

#equation = rbinom(1,2*N,p) / (2*N)

p<-0.1
N<-10
T<-5 #number generations
L<-3 #number independent SNP's
alleles<-matrix(nrow=T,ncol=L) #initialize a matrix of allele frequencies 
each generation
alleles[1,]<-p #initialize first row to equal p

for (j in 1:ncol(alleles)) {
      for (i in 2:nrow(alleles)) {
        alleles[i,j]<-(rbinom(1,(2*N),(i-1))/(2*N))
      }
} 
alleles
...