個別鵯記

【第1393回】 今日も東奔西走しながら【ご恵贈御礼】『Wonderful R 1: Rで楽しむ統計』第一印象など（2016年9月6日）

• 7:00起床。素麺を茹でて朝食。まずは名谷キャンパスに出勤。
• 昨日tweetはしたが，ここでも御礼を書いておく。奥村晴彦『Wonderful R 1: Rで楽しむ統計』共立出版，ISBN 978-4-320-11241-4（Amazon | honto | e-hon）をご恵贈いただいた。普通の統計のテキストとは一味違った構成で面白そう。例えば，「Chapter 4　事件の起こる確率」では富の分布で2項分布をポアソン分布で近似できることを示した後，偶然の再分配を導入すると指数分布に近づくという図を示しておいて「貧富の差のシミュレーションは練習問題として残しておく」という心憎いアクティブラーニング誘導。さらに続いて，地震の確率をポアソン過程と考えた場合の説明からBPTで周期的な発生確率の振動を考えた場合の計算を示し，次にポケットガイガーの測定値を使って「ランダムに事象が起きる」ことのシミュレーションとの類似を見せておいて，『ショパンの「雨だれ」（図4.6）のように一定の時間間隔で起きるのではなく……』と楽譜を提示するセンス。果てはポアソン分布の信頼区間について簡単な求め方と問題点を示し，素粒子物理で使われるようになった新しい方法であるFeldman-Cousinsの信頼区間の話にまで展開する。教科書にはしにくいかもしれないが，Rを自習したい人にはうってつけと思うし，もう10年以上Rを使ってきたぼくでも新しく知ったことも多いので（binom.test()は良く使うが，poisson.test()は知らなかったし，もちろんFeldman-Cousinsの信頼区間は知らなかった），少しでもRやデータをいじることに関心がある人なら読んで考えて手を動かしてみて損はないと思う。なお，本書では付値に<-ではなくて=が使われている。p.3に書かれている通り，かなり前から=も使えるのだが若干問題が生じる場合もあり，ぼくは未だに<-を使ってしまうが，入力の手間はもちろん=の方が少ないわけで，これからRを使う人は=で付値としても良いのかもしれない。また，「Chapter 5　分割表の解析」の中で，オッズ比とp値，信頼区間についてさまざまな求め方が説明されていて，拙作fmsbのoddsratio()も取り上げていただいているが，去年4月や一昨年12月に書いた通り，この関数はp値を求めるときと信頼区間を求めるときに違う考え方を使っているために，時々見かけ上変なことが起こる（という話を他ならぬ奥村先生とかつてtwitterでやりとりしたことがある。本書ではp.79のfisher.test()での非心超幾何分布を使った最尤推定の中身の説明――これ自体，こんなにわかりやすく説明されている本は実はあまりない――に続いて，fisher.test()の場合に同じ問題が起こる話とexact2x2パッケージのfisher.exact()関数を使った解決策について書かれていた）ので，今度のfmsbアップデート（健康余命を求める関数を入れるために近々予定している）には，確実に対処したい。
• 昼から六甲台に移動し，入試の面接の担当分が終わってからは19:00まで研究室で待機（しながら採点業務の続き）。19:10から20分ほど会議があって，この夏の修士課程入試業務は終わり。名谷キャンパスに移動して1時間弱，院生の研究指導（投稿した論文に帰ってきた査読結果に答えて直すところなので，ほぼテクニカルなものだけ）をして，これから帰途に就く。湊川公園廻りの終バス1本前か。
• 今日のドラゴンズは力負けか。高橋周平選手のFCはぎりぎりのプレーだったようだし，仮にファーストに投げて１つアウトを増やしていたとしても1点を失ったには違いないし，そうしていた場合は次打者は併殺ではなく１つのアウトになったはずだから，FCの影響はゼロだったと思う。攻撃も選手はいろいろ工夫していたがカープのセカンドベースマン菊池選手が上手すぎ。
• 都賀川の大腸菌計測のためのサンプリングが悪天候のため2度中止になっている。川の水量や流速は重要な情報だけれども測れないから仕方ないという話になっていたのだが，工学部の藤田一郎教授が開発したドローンから空撮した画像を解析して水量と流速を推定するシステムの存在を知った。これは協力をお願いするべきかもしれない。そもそものシステム開発のきっかけが，ゲリラ豪雨によって都賀川で起こった水害事故だったそうだし。参考情報としてこれとかこれ。とりあえずメールで提案してみようか。

Tweet

△Read/Write COMMENTS

▼前【1392】（東奔西走＋『罪の声』『水族館の殺人』『モンスターマザー』感想メモ（2016年9月5日） ）　▲次【1394】（会議と成績入力と（2016年9月7日） ）　●Top

Notice to cite or link here | [TOP PAGE]