syou6162 https://blog.hatena.ne.jp/syou6162/ yasuhisa's blog https://www.yasuhisay.info/ id:mrknさんがベキ分布を生成するRの関数を書いてくれたので、例のやつを書いてみようかと思いました。 > rzipf function (n, s=1, q=0, supp=c(1,n)) { x <- supp[1]:supp[2] p <- 1/(x + q)**s y <- sample(x, n, replace=T, prob=p) return(y) } > hist(rzipf(10000),nclass=50,plot=FALSE) hist(rzipf(10000),nclass=50,plot=FALSE) $breaks [1] 0 200 400 600 800 1… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fwww.yasuhisay.info%2Fentry%2F20070909%2F1189340158" title="ベキ分布で遊んでみる - yasuhisa&#39;s blog" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> http://farm2.static.flickr.com/1318/1350692964_f6df3e5699.jpg?v=0 Hatena Blog https://hatena.blog 2007-09-09 21:15:58 rich https://www.yasuhisay.info/entry/20070909/1189340158 1.0 100%