MIBIM
laboratoria 5. - quasispecies

n=10;r=8;f1=1.;f2=0.2;μ=0.01;
MS=Table[1,{r}];
populacja=Table[1,{n},{r}]

If[populacja[[1]]==MS,cos1,cos2]

Part[populacja,4]
populacja[[4,3]]

RandomChoice[{f1,1-f1}->{0,1}]
RandomChoice[{f2,1-f2}->{0,1}]

randi=RandomInteger[n];

If[populacja[[randi]]==MS,RandomChoice[{f1,1-f1}->{populacja[[randi]]=Replik[RandomChoice[populacja],μ],}],RandomChoice[{f2,1-f2}->{populacja[[randi]]=Replik[RandomChoice[populacja],μ],}]]

populacja

Replik[l_,mu_]:=RandomChoice[{mu,1-mu}->{1-#,#}]&/@l

Replik[RandomChoice[populacja],μ]