Hi, this is my first time using R. I want to simulate the following process:
"in a population of size N, there are i individuals bearing genotype A, the
number of those bearing A is j in the next generation, which following a
binominal distribution (choose j from 2*N, the p is i/2*N), to plot the
probability of the next generations, my script is as follows. It cannot run
successfully, declaring that the "ylim should be limited. " I wonder where
the bug is. Thanks very much!
generation<-function(i,N)
{
m<-1;gen<-numeric();
for(m in 1:50)
{
testp<-runif(1,0,1);
j<-0; sump<-0;
while(sump < testp)
{ sump<-sump+dbinom(j,2*N,i/(2*N));
j<-j+1;
}
i<-j; gen[m]<-j/(2*N); m<-m+1;
}
plot(m, gen[m]);
}
--
View this message in context:
http://www.nabble.com/-help--simulation-of-a-simple-Marcov-Stochastic-process-for-population-genetics-tp19085705p19085705.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
