Top

個別メモ

Latest update on 2012年3月5日 (月) at 10:54:46.

【第2105回】 また雪（2012年2月2日）

• 今朝は地面が白かったが，それほど深くなかったので，あっさりと雪かきは終わった。微妙に降っているが郵便局で急いで投函しなくてはいけないものがあるので，自転車で出発し，中央郵便局経由で長野駅善光寺口側の夜間駐輪場を利用した。往路あさま510号に乗れて良かった。
• 高崎まで集中して今日の一木会発表準備ができた。
• 某質問紙調査の結果出力の方法を忘れていて，学生から質問されたのだが，思い出すのが大変だった。少し試行錯誤して，右ページは1枚の空白ページの後に1ページとして存在したのでA3横に割り付け印刷を必要部数しておき，それをプリンタの用紙カセットに再び入れてから左ページ（元はA4で連続している）を，用紙をA3横にして割り付けしないし拡大縮小も調整も無しで印刷することで成功した。
• デスクトップマシンで計算させている結果が今ひとつ思わしくない理由として，某所に正規分布族の分布を使っていることが悪い可能性に気づいた。欲しい形は頭に浮かんでいるのだが，どういう分布か名前がわからないので，柴田義貞『正規分布　特性と応用』東京大学出版会をパラパラとめくっていたら，エクスポ・ガンマ分布という分布が目に飛び込んできた。本に載っているPDFの式をRで次のように定義してみた。

  dexpgammab <- function(x, alpha=0.5, location=0) { 1/gamma(alpha)*exp(alpha*(x-location) - exp(x-location)) }

• curve(dexpgammab(x, 0.5, 2), -2, 5)などとして関数形を表示してみたら希望通りだった。欲しいのはCDFなのだが，前掲書にはCDFは載っていないようだった。もちろん積分すればいいのだが面倒だ。それに，たぶんCRANのパッケージのどれかには入っているだろう。しかしCRANを"exp gamma distribution"で検索しても何もヒットしない。仕方ないのでGoogleで検索したところ，MathematicaにはExpGammaDistributionという名前で入っていて，詳細説明を見ると，"generalized extreme distribution"としても知られていると書かれていた。なるほど，これは一般化極値分布だったのか。改めてCRANを検索するとVGAMパッケージに入っていて，dgev(x, location=2, scale=1, shape=-0.5)としたら，だいたいさっきのエクスポ・ガンマ分布の結果と近い形になった（shapeをマイナスにするという発想に至るまで少し試行錯誤が必要だった）。後は，pgevを，元々のコードのpnormだったところで使えばいい。これは，もしかしたら，大きな進歩かもしれない。
• 一木会プレゼンの英訳が終わったのは17:30を過ぎていた。それからLBP-1820の割り付け印刷機能を使って9枚／A4横として両面印刷を2枚という配付資料を20部作った。30分では上限くらいの量だと思っていたが，案の定そうであった。さすがにこのネタのバリエーションで喋るのは人類学会，地球研に続いて3度目なので，スムーズに説明できたし，質疑もまあまああったが，乗ろうと思っていたバスには間に合わなかった。時間が余ったのでリリカの中のスーパーで晩飯を探していたら，大間のマグロをふんだんに使った太巻きが賞味期限が近づいたために半額になっていたので購入。バス停で食べたが美味であった。復路あさま549号。夜間駐輪場は何段にも自転車が停められていて引っ張り出すのが大変だったが，確かにこちらの方が雪に埋もれないし鍵も凍らなくていいかもしれない。
• 帰宅後，息子の受験対策で小論文問題の添削をしてやるのに1時間半ほどかかった。問題は面白いのだが，赤本の解答はポイントを外している場合があるので要注意だ。たぶん高校生レベルではかなりの難問だと思うので，合格発表後に，解答例と解説も含めて大学が公表すべきではなかろうか。

▼前【2104】（もう2月（2012年2月1日） ）　▲次【2106】（民族衛生学会編集委員会（2012年2月3日） ）　●Top

△Read/Write COMMENTS

Notice to cite or link here | [TOP PAGE]