NetaTaneMenu >>>> [[モンティ・ホールのジレンマ]] >>>> [[いくらもらえる？　−　確率を解釈する]]+ [[ビリヤードのパズル]] >>>> [[仮説検定]] >>>> コーヒー豆ブレンド検定
*珈琲豆のブレンドの状態はどの程度同じなのん？
----
[[仮説検定]]の話のついでに提供した珈琲豆の話。
ドンガバチョの「製品ホゲーX」のシムは挫折？したので、せめてこちらの方を。
では問題
----
COLOR(#fe891c){問題は}毎回同じブレンドの珈琲を飲んでいると言えるか？COLOR(#fe891c){です。あるいは，}毎回同じブレンドの珈琲を飲もうと思えば、どの程度の数の豆を一度に挽けば良いのか？COLOR(#fe891c){です。}
----
そこで「十進ベーシック」ソースはこんなのです。

 !ブラジル５０％、メキシコ２５％、キリマンジャロ２５％
 LET  n=10000 !豆10000個
 LET  a=50    !豆1の割合(%)
 LET  b=25    !豆2の割合(%)
 LET  c=25    !豆3の割合(%)

この配合は，どこかの珈琲屋さんのブレンドです。
って豆ってこんなに大雑把なのでしょうか？それはさておき，

 LET  bean$=""
 !豆を測る。
 LET  ma=INT(n*a/100)
 LET  bean$=REPEAT$("a",ma)
 LET  mb=INT(n*b/100)
 LET  bean$=bean$&REPEAT$("b",mb)
 LET  mc=n-(ma+mb)
 LET  bean$=bean$&REPEAT$("c",mc) 

豆は文字列で。最初は配列にしようかとも思ったのですが，
シャッフル（ブレンド）が面倒だと思い文字列にしました。
この時点で，Rubyの方が楽なのにと思いましたが，
今年の職場のコンピュータにはRubyをまだ入れてません。致し方なし。

 !PRINT bean$
 !ブレンドする。
 RANDOMIZE

乱数の初期設定です。

 FOR i=1 TO n
   LET  l=INT(RND*n)+1
   LET  r=INT(RND*n)+1
   IF l>r THEN
      LET  dum=l
      LET  l=r
      LET  r=dum
   END IF
   LET  bean$=bean$(l:r)&bean$(1:l-1)&bean$(r+1:n)
 NEXT i
 PRINT bean$

トランプのように文字列を繰ります。
この回数も問題だとは思いますが...

 !テスト
 LET  tn=100     ! 資料とする豆粒の個数 100
 LET  tt=20      ! 試行の回数

で，混ぜ終わった豆から一掴み取り出して，
調べます。
これを20回繰り返します。

 FOR j=1 TO 10
   LET  sd=0
   FOR ti=1 TO tt
      LET  s=INT((n-tn)*RND)+1
      LET  tbean$=bean$(s:s+tn-1)
      !      PRINT tbean$
      LET  cnta=0
      LET  cntb=0
      LET  cntc=0
      FOR i=1 TO tn
         SELECT CASE tbean$(i:i)
         CASE "a"
            LET  cnta=cnta+1
         CASE "b"
            LET  cntb=cntb+1
         CASE "c"
            LET  cntc=cntc+1
         END SELECT
      NEXT i
      LET  d=SQR((cnta-INT(tn*a/100))^2+(cntb-INT(tn*b/100))^2+(cntc-INT(tn*c/100))^2)
      LET  sd=sd+d
      !      PRINT  cnta;cntb;cntc;d
     NEXT ti
   PRINT sd/tn
 NEXT j
 END

というセットを10回繰り返します。
取り合えず，評価は期待される豆の個数との距離の平均を取ってみました。
では，結果です。
----

 1.25391339658267 
 1.3575955160437 
 1.29658418011372 
 1.33927642901613 
 1.48064839386292 
 1.08039141138613 
 1.39465573097473 
 1.40410594821044 
 1.58761416119748 
 1.45843188521922 

----
評価の仕方もですが，ブレンドの割合，豆の種数，ブレンドする量，取り出す量，
少なくともこの４つぐらいは数値を動かして色々試して見ましょう。（って他人事か？）
COLOR(#fe891c){個人的には前にも話題になったような気がしますが，混ぜ方を変えるとどうかってことが気になります。つまり，トランプのシャッフルのような混ぜ方が妥当かどうか？ですね。うーん。謎が謎を呼ぶ。}
トップ   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS