たつをの ChangeLog : 2014-04-22

スタージェスの公式(Sturges' formula)とは、ヒストグラムを書くときに、観測値の数 n に対して、柱の数（ビンの数、階級数）k の最適値、とは言わないまでも、そこそこ有用な値を求める手段の一つ。

公式と計算例

スタージェスの公式：


Perl による計算例。

% perl -MPOSIX -le '$n=shift; $k=1+log($n)/log(2); print ceil$k' 10
5
% perl -MPOSIX -le '$n=shift; $k=1+log($n)/log(2); print ceil$k' 1000
11
% perl -MPOSIX -le '$n=shift; $k=1+log($n)/log(2); print ceil$k' 100000
18
% perl -MPOSIX -le '$n=shift; $k=1+log($n)/log(2); print ceil$k' 10000000
25

つまり、観測値が1000個ある場合は、11本の柱（ビン、階級）によるヒストグラムを作ればよいのである。

なお、柱の幅hは以下の式で求める。

ヒストグラム出力スクリプト

以下、スタージェスの公式を利用した、テキスト形式のヒストグラム出力スクリプト(Perl)。入力は一行一観測値。階級数はスタージェスで自動決定。

■ソース(hist.pl)：

#!/usr/bin/perl
use strict;
use warnings;
use POSIX 'ceil';
use List::Util qw(max min sum);

my $bar_width = 60;

### データ読み込み
my @vs;
while (<>) {
    next if not /^(-?[\d\.e-]+)/;
    push @vs, $1;
}

### 値補正：平均値が1以下だったら最小値を1以上にすべく
### すべての値を10^X倍する。
if (sum(@vs)/@vs <= 1) {
    my $l = abs(int(log(min(grep {$_} @vs))/log(10)))+1;
    @vs = map {$_*(10**$l)} @vs;
}

my ($max, $min) = (max(@vs), min(@vs));
my $k = ceil(1 + log(@vs)/log(2)); # 階級値
my $h = ceil(($max - $min) / $k); # 階級幅

### 各ビンの頻度
my @freq;
foreach my $t (sort {$a <=> $b} @vs) {
    my $a = ceil(($t-$min)/$h) || 1;
    $freq[$a]++;
}

### ビンの値の範囲（表示用）
my @kgr = ($min);
for (my $i = 1; $i <= $k; $i++) {
    $kgr[$i] = $kgr[$i-1] + $h;
}

my $max_f = max(map {$_||0} @freq);
my $keta_f = length($max_f);
my $keta_m = length(int($kgr[-1]));
for (my $i = 1; $i <= $k; $i++) {
    my $f = $freq[$i]||0;
    my $bar = "*" x ($max_f <= $bar_width ? $f : $f/($max_f/$bar_width));
    printf "%02d %${keta_f}d ( %${keta_m}d-%-${keta_m}d ) %s\n",
    $i, $f, $kgr[$i-1], $kgr[$i], $bar;
}

※平均値が1以下のときに、最小値が1以上10以下になるように全体の値を補正（10のN乗をかけている）している。

■実行例：

% cat hist.test
1
2
3
4
5
20
% head -5 hist.test | ./hist.pl 
01 2 ( 1-2 ) **
02 1 ( 2-3 ) *
03 1 ( 3-4 ) *
04 1 ( 4-5 ) *
% head -6 hist.test | ./hist.pl 
01 5 (  1-6  ) *****
02 0 (  6-11 ) 
03 0 ( 11-16 ) 
04 1 ( 16-21 ) *

ランダム：

% perl -le 'print int(rand(100)) for(1..10)' | ./hist.pl 
01 3 (  8-26 ) ***
02 1 ( 26-44 ) *
03 2 ( 44-62 ) **
04 2 ( 62-80 ) **
05 2 ( 80-98 ) **
% perl -le 'print int(rand(100)) for(1..100)' | ./hist.pl 
01  6 (   0-13  ) ******
02 20 (  13-26  ) ********************
03  7 (  26-39  ) *******
04 13 (  39-52  ) *************
05 20 (  52-65  ) ********************
06  8 (  65-78  ) ********
07 13 (  78-91  ) *************
08 13 (  91-104 ) *************
% perl -le 'print int(rand(100)) for(1..1000)' | ./hist.pl
01 100 (  0-9  ) ************************************************************
02  90 (  9-18 ) ******************************************************
03  98 ( 18-27 ) **********************************************************
04  91 ( 27-36 ) ******************************************************
05  84 ( 36-45 ) **************************************************
06  95 ( 45-54 ) *********************************************************
07  81 ( 54-63 ) ************************************************
08 100 ( 63-72 ) ************************************************************
09  87 ( 72-81 ) ****************************************************
10  78 ( 81-90 ) **********************************************
11  96 ( 90-99 ) *********************************************************

0以上1未満のとき（補正済み）：

% perl -le 'print int(rand(1)*100)/100 for(1..50)' | ./hist.pl
01  8 (  0-14 ) ********
02 10 ( 14-28 ) **********
03  5 ( 28-42 ) *****
04  8 ( 42-56 ) ********
05  6 ( 56-70 ) ******
06  8 ( 70-84 ) ********
07  5 ( 84-98 ) *****

とりあえず深く考えないで数字をつっこめばグラフを作ってくれるので結構便利かと思われ。

参考文献

- Wikipedia:ヒストグラム
- 東京大学教養学部統計学教室 (編集), "統計学入門", 東京大学出版会, 1991. (p.22)

買物メモ。

我が家のたためる携帯コップ[2012-04-17-3]に黄色が仲間入り。既存のオレンジは子供用、青は防災カバンに常備。

■マーナ たためる携帯うがいコップ イエロー


あと、100円ショップのたためるコップ[2013-10-19-2]もさらに2個買い増ししました。主に防災用。質の違いはあるけど、非常時には気にしないしね。

2014年4月22日のヲハニュースをお届けします。



▼GunosyでもSmartNewsでもなく、Pressoを愛用するたった１つの理由 - らっこのじゆうちょう
ソーシャルニュースアプリ「Presso」[2014-03-27-4]。私の愛用の理由は「その場ではてブできる」ってこと。重要。

▼How to Turn Your Phone Into a DIY Photo Projector for $1 | Photojojo
iPhone のバックライトを使った虫眼鏡式アナログプロジェクター。回転固定（逆さまに固定）で正しくうつる（左右も上下も反転するから）。そのうちやってみたい。

▼ダイソー風なブランド名とユニクロ風なロゴで無印良品風な商品を売る日本で大人気のショップ「MINISO 名創優品」がアツい : 市況かぶ全力２階建
おもしろーい！ 中国における日本のブランド力という観点からもおもしろい。

▼少年野球がトラウマになっている - 脱社畜ブログ

野球少年団では、「怒鳴る」「叩く」といった前時代的な指導方法があたりまえのようにまかり通っていた。僕もよく怒鳴られたし、叩かれた。叩かれたあとにお礼を言わなければならないシステムは、今思い返しても謎すぎる。あれは普通に人権侵害なのではないだろうか？

今考えると、大人たちの「怒鳴る」「叩く」には彼らなりの論理があったのだろうとは思う。もしかしたら、「お前たちのためを思って、あえて厳しくしているんだ」と彼らは思っていたのかもしれない。「厳しくしなければ、強くなれない」という信念があったのかもしれない。この手の論理は、大人になった今でも出会うことがある。会社の新人を散々詰った挙句、「厳しくしてるのは、お前たちに早く成長してほしいからだ」と平気な顔をして言う人は割とよく見かける。

ほんと、うんざりだよねー。

▼自己啓発書のリストを作って定期的にポストしてればプロブロガーになれる理論 - あざなえるなわのごとし
最近やってないなあ。またやるかな。

▼なんで、たわしが東大に!? - バッタもん日記
東大にタワシが。わたしではなくたわし。

▼『USB半田ごて2』の実力に迫る ネタ商品かと思ったらかなりの実用品でした！
良さげー。
どこでも電源とれる（海外でも）。

▼「口にお金を入れる」のはやめよう-紙幣に膨大なバクテリア - WSJ.com
言われなくてもやらないよー。

▼夫婦でカフェで過ごすということ（日々の日記） | 結城浩 | note
代官山のツタヤ・スタバが舞台。

▼Amazon.co.jp： Anker 40W 5ポート USB急速充電器 ACアダプタ PowerIQ搭載 iPhone5C/5S/5/4S/4/iPod/iPad/Xperia/GALAXY/ウォークマン等対応: 家電・カメラ

売れてるし、欲しくなる。今は充電器あまってるくらいだから不要だけど。

マイボトルドリンク「drop（ドロップ）」に初夏の新フレーバー「トロピカル烏龍茶」が登場。4月8日から発売されています。

- マイボトルドリンク「drop（ドロップ）」
http://www.suntory.co.jp/softdrink/drop/

で、drop の先行モニターをやっている関係で[2013-07-29-1]、この新しいのが入ったポーション12個セットを頂きました。毎回ありがとうございます。



「トロピカル烏龍茶」

トロピカルフルーツの甘い香りがリラックス気分へ誘う、香り立ち華やかな烏龍茶

トロピカルな香りの無糖烏龍茶。これまでと同様に[2013-07-29-1]、独断的な感想と評価。

■トロピカル烏龍茶
☆
夏っぽいトロピカルな香りが広がります。ホットで飲んだのですが、これは断然アイスが良さそうです。氷をいっぱいいれてキンキンに冷やして飲みたいですね。
 
というわけで、私の中のTOP3である「HOTでおいしいオレンジ」「あまとろ梅」「抹茶入り緑茶」（どれも☆☆☆）はまだまだゆるぎません。

12個セットの変遷メモ。今回は季節もののみ入れ替え。
- 2014年2月
-- OUT「香る烏龍茶」「芳香カシスコーヒー」
-- IN「さくら香る花茶」「甘い香りのマカダミアコーヒー」
- 2014年4月
-- OUT「さくら香る花茶」
-- IN「トロピカル烏龍茶」
 

関連記事

 
- 【お手軽ドリンク環境】サントリーとサーモスの新製品「マイボトルドリンク drop（ドリップ）」のブロガーイベントに行ってきました！[2013-07-29-1]
- マイボトルドリンク Drop（ドロップ）専用ボトルのボトムカバー！[2013-08-16-3]
- 「マイボトルドリンクdrop」に新しいポーション（さくら＆マカデミア）が登場！[2014-02-08-1]

先週のことになりますが、サントリーさんからお誘いいただいて、メーカーズマークのディッピング体験イベントに行ってきました。いつもありがとうございます。



メーカーズマークはコーンを主原料としたウイスキー、バーボンです。バーボンといえば、昨年のジムビームのイベント[2013-02-23-2]で「再開」したばかりですが、晩酌ではさりげなくローテーション入りしております。で、メーカーズマークはプレミアムバーボーンと呼ばれる高級カテゴリです。高いけどおいしい路線ですね。

甘さがポイント。バーボン自体、甘いのですが、メーカーズマークはしっかりと上品な甘さです。ジムビームも甘い系統でしたが、メーカーズマークの方が甘いですね。単に甘いだけでなくきっちり美味しいのが素人で飲み慣れない私にも分かります。ブラインドテストしても分かるよ。イベントでは他社のバーボンとの飲み比べも行いましたが、断然メーカーズマークが好みです。これは間違いない。あと、ソーダで割って飲むよりも、ロックやストレートが良いです、かなり。



§

で、イベントの本題。このメーカーズマークのビンは赤い蝋で密封されています。この「赤い封蝋（ふうろう）」がトレードマークみたいな感じですよね。一本一本手作業でやっているそうです。200度くらいで溶かしたロウにビンの口をさっとつけるという作業。なのでロウの垂れ具体が同じものは1本しかないとのこと。そりゃそうですね。

それを参加者全員が体験しました。自分のオリジナルのメーカーズマークを完成させるのです！ やってる様子、細かい作業内容等は他の参加者のブログ（末尾）をご覧下さい（自分では撮ってなかったので…）。

私がやったらこんなふうになりました。



漬けすぎ感がありますが、他の人も同じような状態。やはり、ロウをたくさんつけたくなるんですよね。モッタイナイ精神かも。でも、ロウがたくさんの方が素人風になってオリジナル度がアップでとっても素敵です。



§

会場はジムビームのイベントのときと同じ渋谷の NOS ORG でした。イベント後半、食事＆メーカーズマークを頂いたのですが、ソーダ割りでスタートしたはずがいつのまにやらロックになって、美味しく楽しく飲食できたのは良いのですが、結構酔っぱらってしまいました。翌日は軽く二日酔い。



- NOS ORG（ノス オルグ）
http://org.nos-tokyo.com/

§

オリジナル封蝋のメーカーズマークを一本頂いて帰りました。もったいなくて開けられない感じですが、心を鬼にして(!?)、一ヵ月以内にあけちゃう予定です。あけちゃう、って開封だけで、中身は大事に少しずつ飲みますよ。

そういうわけで、あらためて、ありがとうございました！ ごちそうさまでした！ ほんと楽しかったです。

§

■参考サイト（オフィシャル）：

- メーカーズマーク　世界に一本しかない、あなた宛てのバーボン。　サントリー
http://www.suntory.co.jp/whisky/makersmark/

- [SUNTORY公式ブログ] 【体験レポ】「メーカーズマーク」ディッピングイベントに参加してきました！ (サントリートピックス)
http://topics.blog.suntory.co.jp/006158.html

■いっしょに参加した方々のレポート（随時追加）：

- [N] 世界に1本だけのウイスキーを作る！ボトルを封蝋する「メーカーズマーク（Maker's Mark）」ディッピング体験に参加
http://netafull.net/whisky/046432.html

- それはエレガントな年上の女性のような　「メーカーズマーク（Maker's Mark）」ディッピング体験イベント」:裏[4k](ura_shike)落語・酒・デジタル　など
http://shike.blog.jp/archives/51938806.html

■メーカーズ マーク 45度 750ml ［バーボン］

ヤフーのリアルタイム検索関連アプリで、ツイッター等での話題（バズトピック）をいち早く手軽に閲覧できる「話題なう」ですが、このたびのアップデートで名前が「リアルタイム検索」になりました。基本的には今までの「話題なう」アプリと同じです。



- Twitter検索の決定版！Yahoo!リアルタイム検索公式アプリ 話題なう - Yahoo! JAPAN
http://promo.search.yahoo.co.jp/realtime/wadainow/

新しいキャッチフレーズ（？）は「ログイン不要！ Twitter検索の決定版！」とのことです。バズってる話題のチェックや今さっきのことの検索などに便利です。

新機能として、繰り返し検索するワードをお気に入り登録できるようになり、かつ、そのワードがバズったら通知してくれるようになったようです。定期的にウォッチしたいワードを登録しておくとネットウォッチがはかどるかもー！

関連記事

- 「話題なう」アプリのテレビ限定機能と通知機能[2014-02-07-4]
- Yahoo!リアルタイム検索に「感情分析(β版)」が登場[2013-09-12-3]
- Yahoo!リアルタイム検索の「注目のキーワード」が分かりやすく閲覧できるアプリ「話題なう」が iPhone にも登場[2012-09-05-2]
- 日本のヤフー(Yahoo! JAPAN)でツイッター検索ができるようになりました[2011-06-14-2]