Rを使用して集団の遺伝子頻度をシミュレートするために、ランダムな頻度(つまり、頻度は1に等しい必要があります)を生成したいと思います。私の解決策は次のとおりです。
freq<- function(x,y) #random frequency generator, no defined distribution
{rn<-sample(1:y,x)
sum <- sum(rn)
for(i in 1:x){
RG[i]=rn[i]/sum}
return(RG)
}
除算前に合計を特定の値 (たとえば、乱数は 100 に等しくなければならない) に制限する方法に関する提案はありますか?