メモ Ruby対数正規分布の乱数

(2017/11/16)あああタイトル間違ってましたぁ・・更新前に見た方申し訳ないです・・

合っているかは自信ないです・・

ゲームの乱数くらいにお使いください。

require "random_bell"

# http://pmsl.planet.sci.kobe-u.ac.jp/~seto/?page_id=316
# をRubyに移植
def lognormal mu, sigma
    return Float::NAN if (mu <= 0)
    bell = RandomBell.new(
        mu: Math.log(mu*mu) - Math.log(mu*mu + sigma*sigma) / 2.0,
        sigma: Math.sqrt(Math.log(1 + (sigma / mu) * (sigma / mu))),
        range: -10..10)
    Math.exp(bell.rand)
end

分布確認用

require"./lognormal"
puts (1..10000).to_a.map{|n|lognormal(平均,分布).to_i}.inject({}){|r,n|r[n]||=0;r[n]+=1;r}.sort_by{|n,v|n}.map{|n,v|"#{n} #{"*"*v}"}

C++(標準ライブラリ)にもPython(NumPy)にもC#(Math.net)にもあるのに
Rubyにはない悲しみ

gem install random_bell random_bellを使用しています。

random_bell (0.1.1) ruby 2.3.3p222 (2016-11-21 revision 56859) [x64-mingw32]

Ruby 呼び出せないメソッド

define_method("!bar", ->{"foo"})

p method("!bar")

p method("!bar").call # 呼び出し方

!bar #< 呼び出せない

こんなことができるとは思ってなかった

ruby 2.3.3p222 (2016-11-21 revision 56859) [x64-mingw32]にて確認

Hamlの組み込み

Hamlをテキストから直接実行したいときのメモ

require "haml"
template = <<'EOS'
%p= a1
%p= a2*3
- a3.each do |i|
  %p= i**2
EOS

haml_engine = Haml::Engine.new(template)
puts haml_engine.render(Object.new, a1:123, a2:10, a3:[1, 2, 3, 5])

ヒアドキュメントはインデントが合わなくなるから嫌いです・・・

参考

File: REFERENCE — Haml Documentation

Ruby:実行中のソースコードのファイルがあるディレクトリ

いいページがあるものの、検索でなかなか引っかからないので自分でページを作ってみた。

d.hatena.ne.jp

これのものです。

File.expand_path(File.dirname(__FILE__))

100記事目！(この記事は101記事目)

先ほどの

soukouki.hatenablog.jp

が100記事目になります！

17ヶ月、少しずつためた記事です。

アクセスもだんだん増えていき、2847アクセスになりました。(こんな内容の記事でいいのか・・？

誰かの役に立てる記事が書けるようにしていきたいです。