ぎょーむ日誌 2004-11-28

2004 年 11 月 28 日 (日)

• 0900 起床． 朝飯． コーヒー．
• 怠業 …… ばかりしてるわけにもいかんので， 他のことはさておき， とどこおっていたメイル返信． まずは間瀬さんにいろいろと教えていただいた御礼 …… で， とても trivial な内容なのに 1 時間以上を費す． うーむ， われながら遅い．
• これまたとどこおっていた， glmmPQL() に関する質問 (のひとつ) の回答を検討する …… 結果に Log-likelihood とか出てまぎらわしいな． しかも summary() とると， ご丁寧にも AIC & BIC 計算している． ふーむ． ともあれ返信．
• glmmPQL() は一般化線形モデルを quasi-likelihood (QL) ふうに 定式化して， lme() (in nlme) に渡してるだけなんで計算とかは確実っぽく見える (う? PQL の pennalized はどこにいった?)． そうか， なんでパラメーター推定値に関して t 検定やるのか， といった疑問については nlme 方面を調べないとわからないわけだな． いや， 最大 QL 推定値の分布ってわかってるのかしらん?
• ともあれ， 滞納なメイルは …… あ， アカマツ作図がまだ残ってた． とりあえずこれは明日の夜以降， ということで．
• 1205 自宅発北大構内走． 晴れたり雪ふったり． 気温はかなり低い． 1250 帰宅． 昼飯．
• 1350 自宅発． 晴． 1400 研究室着．
• 時間もなくなってきたんで， 明日の統計学授業の準備．
• 自然科学について何ごとかを修めるべく入学してきた ここ地環研大学院生は， 2 年間の修士課程の中で「データをどう数理的にとりあつかうか」 について 2 回分の講義しか受けない …… 他ならぬ私のやる これ だけ． しかも全員が受講するわけでもない．
• 受講生の大半は学部のときにも統計学とかやっていない (少なくとも「理解できる統計学授業」は受けてない) ので， 内容的には高校生に教えるつもりでやらねばならない． しかも， 高校数学を忘れてしまった高校生あいてに …… 本邦ではなぜこれほどまでに 「データ解析なんて簡単，やればできる」 という信仰が流布してしまったのか， なかなかナゾなところだ．
• で， 前回は統計モデルと推定の関係を理解しよう， と強調したんで， 今回もしつこくそのあたりのいろいろな例を挙げて， 「データの傾向とばらつきを図にする → 推定の試み」 という基本動作 を浸透させたいところなんだけど …… まあ， 「検定」の説明します， と予告したのでその準備をせんといかん．
• 「検定」などという難しいことは， (「理解できる統計学授業」を一回しか聞いたことない) 彼らにはまだまだ時期早尚というべきなんだけど …… その指導教官たちは 「データを統計学的に解析しろ」 とは言わずに 「検定にかけろ」 などと乱暴な 「指導」しやがることが多いので， 私としては先まわりにしてテをうたずにはいられない．
• で， その「検定」の説明であるが …… 昨年は連続確率分布から得られた標本に対する検定なるものを説明した． ありていに言って， 連続確率分布は説明するのもされるのも難しい． そこから得られる標本統計量の分布もまた連続分布なので， さらに難しい． このふたつについて院生が理解できるか， もしくは理解できたふりができなければ， 「検定」なんぞというしち面倒なしろものの説明は さらに難しくなる．
• ということで， 今年は離散な事例での説明を試みてみよう． 離散確率分布の密度関数は「度数グラフ」として簡単に説明できるんで， 確率なるけったいな概念と観測値たるカウント数の対応が断然つけやすくて， 説明するのもされるのもラクなのである． また R でプログラム作り …… 今ごろやってて間に合うのか?
• 1620． なんか進捗しないと思ったら， かとーオフィス内， けっこう寒くなっていた． ガスストーヴつける．
• すでに夜． そっかー， 離散確率分布の場合， 危険水準◯% の棄却域， というのを設定するときに 5% とかそういうのを 厳密には図示はできんわけだな． おもしろい．

[コイン投げ 20 回 × 一万回]

むろんこれはシミュレイションなんぞやらずとも， 単純な二項分布になるってのはわかるわけだが (R だと dbinom(0:20, 20, prob = 0.5)) …… 試行回数が 1000 ぐらいだと， ぜんぜん二項分布っぽくなくて面白い． で， 10000 とかにするんだったらむしろ barplot() すべきだったかも， と後になってから思いついたり．

• どうも EPS 形式で出力させるのはマズそうなので， jpeg() 出力と jpeg2ps 連携わざでいくことに． jpeg() では日本語が出力されんか． じゃ， 敵国語で， と調べてみる …… コインのオモテは obverse (俗語では heads)， ウラは reverse (tails) …… いや軸ラベル問題は LaTeX で片づけるか．
• ぢりぢりと進捗させるが前途遼遠． はらへったので撤退． 2040 研究室発． 寒い． 2050 帰宅． 晩飯．
• 晩飯後も準備はつづく． 明日はなんだかいろいろありそうなんで， あまり時間がない． ということで， 本日おいつめられてるわけ．
• 関係ないことにふと気づく． 先日， 提出した自由集会の要旨， 自由集会 web page の URL 書くの忘れてた …… ま， いいか．
• 準備が終わらん．
• 時刻 2732 (0332)， けっこうでかい地震． 震度 3 ぐらいかな． ラジオによると釧路あたりは震度 5 強だとか ……
• ときどき余震． お， 震度地図 でたか (tenki.jp で更新されてる)．
• さーて， 授業の準備のほうは …… 7 割ぐらいできたか? もう寝たいなぁ． 時刻すでに 2830 すぎ．
• [今日の運動]
  • 北大構内走 1205-1250． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0930): リンゴ．
  • 昼 (1320): 米麦 0.7 合． モヤシ・セロリ・ネギ・シイタケ・卵の炒めもの．
  • 晩 (2150): 米麦 1.0 合． セロリ・ニラ・エノキダケ・豆腐・卵の炒めもの．