ぎょーむ日誌 2012-11-(21-30)
2012 年 11 月 21 日 (水)
-
0650 起床.
朝飯の準備.
朝飯.
0740 自宅発.
晴.
気温ひくく,
水たまりなどは完全に凍結している.
0810 研究室着.
コーヒー.
-
統計モデリング入門
について「まちがいではないか?」とのご指摘をうけたので,
いろいろと検討してみる.
勉強になる.
-
返信かいて午前終了.
北 7 条郵便局に用事があるので,
ついでに昼飯も外で食べることにする.
天気予報はずれて晴.
ここらへんから見える最高峰である手稲山,
山頂まわりが冠雪しているな.
-
で,
郵便局での用事を済ませて (このあたりについては後述),
郵便局の横にある蕎麦屋へ.
北大のすぐ南,
北 7 条郵便局のとなりにある蕎麦屋.
以前に一度はいったことがあるはずなのだが
……
ぎょーむ日誌に記録なし.
前世紀のできごとだったのか?
ともあれ,
手打ち蕎麦 + 温かい鶏つゆをたのんでみた.
この店は機械で切る平打ちと手打ちがあり,
後者はかくのごとき田舎蕎麦.
新蕎麦のしーづんなのだが
……
あまり香りがするといったかんじではない.
「蕎麦たべた」という満足感はあります.
[北大生協会館店・書籍部]
この書籍部では何となく
「
統計モデリング入門」
は冷遇されていると思っていたんだけど
(いつ行っても置いてなかったので)
……
今日は分厚い「確率論ハンドブック」と星野さんの「調査観察データの統計科学」
に挟まれて「統計モデリング入門」も置いてありました!
-
なぜかその前にソーバー「科学と証拠」もあり.
ぱらぱらと立ち読み.
本文はなんとなくびみょーかなあと思ったんだけど,
松王さんの解説はおもしろそうな気もするので
……
うーん,
どうしようかなあ.
-
で,
研究室にもどってカネ勘定をすこしばかり.
[岩波書店に 33 万円を払込]
「統計モデリング入門」
割引販売
の本代・送料総額から全印税分をさしひいた差額.
これについては後述
……
それよりおどろいたのは,
ゆうちょ銀行 ATM に
スキャナがついてたこと.
郵便局窓口で「支払いはこのゆうちょ銀行通帳から……」と依頼すると,
「それなら ATM でやってみてください」とのことなので「?」
といぶかしみつつ試してみた.
よく見ると ATM に
払込用紙を投入するスロット
があり,
岩波書店から送られてきた用紙をつっこむと ATM
が「口座番号はこれで OK?」とか「このスキャン画像みて払込金額を入力して」
などと interactive に作業がすすみ,
最後にゆうちょ銀行通帳をいれてパスコードいれたら払込が終了した!
利用明細票はかくのごとし.
……「枯れた技術」のよせあつめとはいえ,
ちゃんと動作してすごいとおもった.
-
あとからいただいたご指摘:
「えー,ずっと前からやってるよ
……
知らなかったの」
「だいたいゆうちょ銀行どうしのお金のやりとりなら,
郵便局に行かなくても web browser からぽいぽいとやればいいし,
それだったら手数料 290 円はゼロでしょ」
……
みなさん,
つよくいきましょう
……
-
さてさて,
上記の金額だが
……
てきとーに数字を列挙してみよう:
-
統計モデリング入門
は現時点まで約 3700 部ほど印刷されてるので,
これが売れようが売れまいが印税は入る.
印税率 8%,
すなわち 3700 部 × 3800 円 × 0.08 で 112 万円ほど.
-
ただし岩波書店は
割引販売
の著者割本の代金 (久保が買い取って割引購入者に売ってるということ)
とその送料を立てかえていて,
それが上記印税で相殺されているのだが
……
現時点で 33.5 万円ほど足りない (ので上のごとく,
久保から岩波書店に払込をした
…… 岩波としては「貸金」を年内に回収したいそーで)
-
いっぽうで割引販売
(定価税込 3990 円を 3000 円)
の代金は久保が受け取っていて,
いままで合計 420 部ほど売れたので
420 部 × 3000 円 (割引価格) で 126 万円ほど.
-
これから今日の払込金を引くと 93 万円ほどになり,
印税と比較すると,
あたかも割引販売によって 112 - 93 ということで 19 万円ほど
ソンをしているように見える
しかしながら,
印税というのは「何部すられたか」で決まるおカネであり,
割引販売で 420 冊も売れてなければ,
第 4 刷 (9/25) の 700 部などはまだ印刷されなかったかもしれないので,
まあ別にソンしてるというほどでもない
……
そもそも採算だの効率だのがどうこうというハナシでもなく,
まあ,
割引販売で誰も不幸にはなっていない,
ということで良いのではないでしょうか.
-
午前の検討その他を反映して,
「統計モデリング入門」まちがいなど
ペイジを更新.
-
Gamma 分布の GLM
に dispersion parameter
phi
についての補足.
これの推定量は,
sum((object$weights * object$residuals^2)[object$weights > 0])/df.r
といったものである.
-
総称関数
summary()
の glm()
版である summary.glm()
のソースコードは
このようなもの
であり,
その中に上のごとく記述されている.
-
これは,
いわゆる「赤 GLM 本」こと Faraway ``Extending the Linear Model with R''
の p.137 に紹介されている推定量
Χ2 / (n - p)
である.
蛇足ながら,
私は
Χ2
を「ああカイ二乗ってことは deviance ってことかな」
などとかんちがい
していたのだが,
じつは「大 Pearson の
Χ2」
であり,
上の R コードのごとき式となる (weight を考慮した場合).
-
さらに補足
(どうせ自分でもすぐにどわすれするので)
……
summary.glm()
のコードを取得するためには
body()
関数を使って body(summary.glm)
などとする.
関数の引数を知りたい時には
args()
……
ながねん R を使っていながら,
こんなことも本日まで知らなかった.
-
1730 研究室発.
買いもの.
1800 帰宅.
晩飯.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0710):
米 0.7 合.
ダイコン・ニンジン・シイタケ・鶏だんごの味噌汁.
- 昼 (1240):
盤渓そばの手打ちあつもり 790 円.
- 晩 (1940):
米 0.7 合.
ダイコン・ニンジン・シイタケ・鶏だんごの味噌汁.
サンマのバルサミコ煮.
コマツナ・油揚げの煮びたし.
2012 年 11 月 22 日 (木)
-
0650 起床.
ラジオ体操.
朝飯の準備.
朝飯.
0730 自宅発.
曇天.
0800 研究室着.
コーヒー.
-
ちょっと脱力ぎみ
……
-
xargs -P が激しく便利なのでメモ
は興味ぶかい記事だな.
xargs -P 8
などと指定すると,
8 個のプロセスに分割して処理される.
-
並列処理はいいねえ
……
R の場合
library(parallel)
を使うのが
かんたん.
-
では,
ついでに Perl では
……
なとど
Parallel::ForkManager
にテをだしたのがうんのつきで,
午前はこれでつぶれてしまって,
しかもかんばしい成果を得られなかった.
-
昼飯.
-
季節モノというべきなのか何なのか
……
データ解析こんさるメイルかき.
-
またわきみちにそれてしまった
……
Ubuntu のコマンドである
usb-creator-gtk
を使って USB メモリスティック boot な Ubuntu づくり.
このコマンドは Ubuntu スティック
(容量は 2GB で OK)
しか作れない.
-
とりあえずふつーにいれると上記のごとく,
GNOME based な Unity desktop になるのだが
……
どうもこいつの動作はもっさりしているね.
-
USB な Ubuntu 作ってる理由は「緊急時の起動用」なので,
もっとさくさくうごいてほしい
……
ということで Xubuntu な ISO ファイルをダウンロードして
……
おお,
usb-creator-gtk
が使えた.
-
上のようなかんじの UI となる
……
ふーむ,
これが標準の Xubuntu (Xfce4 desktop)
か.
ふだんから Xubuntu 使っているんだけど,
desktop な pacakge は
sudo apt-get remove
しまくったので,
ひじょーに殺風景な user interface なんだよね.
とても気に入っているんだけど.
-
またこんさるメイルかき.
-
てなかんじでオワってしまった
……
1730 研究室発.
なぜかしら札幌駅ビルの大丸地下食品売り場に突入して,
買いもの.
すごいこんでいて,
しかも性比がすごくかたよっている
……
ほうほうのていで脱出して,
1810 帰宅.
晩飯の準備.
晩飯.
-
[今日の運動]
-
とらちゃんえくささいず
……
真夜中すぎまでつづく.
背筋・肩が「板」みたいなかんぢ
……
-
[今日の食卓]
- 朝 (0700):
米 0.7 合.
シジミの味噌汁.
サケの塩焼き.
- 昼 (1200):
研究室.
米 0.7 合.
サンマのバルサミコ煮.
レタス.
ゴボウサラダ.
リンゴ.
- 晩 (1920):
米 0.7 合.
カブ・ジャガイモ・ニンジンのスープ.
焼き鳥あれこれ.
コンニャクの煮物.
2012 年 11 月 23 日 (金)
-
0730 起床.
とらちゃんを
Ergobaby
でだっこしつつ,
朝飯の準備.
Ergo 効果で寝てくれた
……
のでテイブルの横でそっと寝かせて朝飯.
[エルゴ]
と通称されるだっこ支援装備.
育児なヒトたちのあいだでは「最強」とのウワサも.
[らぐびーな服]
朝からまた格納容器からの漏れだしをともなうばくはつぶつ事件があったので,
きがえ.
何やらこのような服をいただきました
……
お,
けっこうにあっている.
-
昨日は Ubuntu の
usb-creator-gtk
なるコマンドを使って USB bootable な Linux
(このコマンドだと *Ubuntu 限定だが)
を作ったわけだが
……
もっと汎用な
usb-imagewriter
なるツールがあるようだ.
この
How to install Linux Mint via USB
はなかなkよい.
-
昼飯の準備.
昼飯.
-
午後に何をやっていたのか記憶がないのだが
……
とらちゃんと遊んでた?
雪の中のコープ札幌まで買いものにいってみた.
これはひとりで.
-
晩飯の準備.
晩飯.
-
またねてくれないとらー
[真夜中のとらちゃんだんす]
ねてくれないー
-
[今日の運動]
-
[今日の食卓]
- 朝 (0820):
大丸地下で買ってきたパン.
生ハム.
チーズ.
- 昼 (1300):
スパゲッティーだったはずだが
……
自分で作っておきながら,
内容を思い出せない.
- 晩 (1900):
米 0.7 合.
カブ・ニンジン・ジャガイモのスープ.
コマツナ・マッシュルームの炒めもの.
ゆでブロッコリー,タマネギ・ひき肉のソース.
2012 年 11 月 24 日 (土)
-
0800 起床.
朝飯の準備.
朝飯.
コーヒー.
-
今日はとらちゃんを新千歳空港に送りとどけ,
その後ただちに撤退
……
とゆーナゾのミッションをこなさねればならぬ.
[ということで,初 JR]
いきなり新千歳空港行きの快速エアポートにのせられてしまったとら坊.
その移動速度におどろいているもよう.
[新千歳空港]
広さにおどろいているもよう
……
このあと,
私は授乳室なる不思議な (そして何ともいづらい)
空間で,
ばくはつぶつ処理のおむつがえ作業にとりくまねばならぬことに.
[帰りの快速エアポート]
ぐったりつかれておりました.
帰りの快速エアポートは札幌駅までいって L 特急になる
列車だったので,
ふつーの快速エアポートより座席が快適,
往路よりも乗客密度が低くくてのんびりできました.
-
AEON でちょっと買いものして,
帰宅.
哺乳瓶作業でけっこう時間とられる.
おそい昼飯.
-
ぐったり.
少しだけそうじ.
-
また哺乳瓶 &
ねかしつけ.
-
とらちゃん,
夜になって,
また不定愁訴的な泣き
……
対応するこちらもたいへん.
-
晩飯時刻もどんどんおそくなる.
今日はかんたん晩飯.
-
そしてまた真夜中すぎのとらちゃんだんす.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0830):
米 0.7 合.
ネギ・ナメコ・油揚げの味噌汁.
カブシュート・油揚げの煮びたし.
- 昼 (1420):
AEON 桑園の「パン工房」のバゲットとか.
ネギオムレツ.
- 晩 (2120):
スパゲッティーペペロンチーノ.
2012 年 11 月 25 日 (日)
-
0800 起床
……
昨日はばてばてだったが
……
まあ,
ある程度は回復?
朝飯の準備.
朝飯.
コーヒー.
-
午前中はもっぱら怠業.
天気は晴.
-
今日のミッションは,
とらちゃんを Pluto なる人力四輪車 (六輪車だったか?)
に搭乗させて,
札幌駅方面を強行偵察してみること.
-
駅ビル南部である ESTA B1F にあるフードコート
……
時刻はすでに 14 時ぐらいなんだが,
かなりこんでる!
なんとか場所をみつけて昼飯.
-
札幌駅周辺をふらふら.
ヨドバシとか.
-
自宅方面に戻る.
コープさっぽろ植物園店で買いもの.
晩飯.
[Linux 練習中]
ThinkPad をたたくとら先生.
-
[今日の運動]
-
とらちゃんえくささいず
……
今日は Pluto による移動がおおかったのであまりだっこしなかったわけだが
……
それは「休日のわりには」ということで,
真夜中になるころにはまた背中がけっこうやられました
……
ともあれ,
「三連休」はおわった
……
-
[今日の食卓]
- 朝 (0700):
パン工房のバゲットとか.
生ハム.
- 昼 (1430):
ESTA フードコートで石焼ビビンバ.
- 晩 (1830):
米 0.7 合.
ビーフストロガノフ.
トマト・レタス・生ハムのサラダ.
2012 年 11 月 26 日 (月)
-
0710 起床.
朝飯の準備.
朝飯.
コーヒー.
0810 自宅発.
曇天.
気温はやや高め,
過去数日にくらべれば.
0840 研究室着.
-
なんだかメイルなあれこれ
……
-
……
で,
午前終了か?
昼飯.
ちょっと外をうろうろ.
-
来年の 8 月で統計モデリングの「なつやすみ講習」
みたいな授業の一部を分担させていただくことになりました.
-
せとぎわな院生のデータ解析こんさる
……
でオワってしまったかんじ.
-
1800 研究室発.
買いもの.
1830 帰宅.
晩飯の準備のてつだい.
晩飯.
[本が売れる季節?]
「R グラフィックス」
の
割引販売
でひさしぶりに一冊うれました.
「
統計モデリング入門」
も順調なようで,
岩波書店から第 5 刷を来月だす予定という連絡いただきました.
合計印刷数 4700.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
米 0.7 合.
ネギ・油揚げの味噌汁.
- 昼 (1240):
なぜか正門前の CoCo 壱番でソーセイジカレーを食べてしまった.
- 晩 (1910):
イサキのアクアパッツァ.
リゾット.
2012 年 11 月 27 日 (火)
-
0650 起床.
朝飯の準備.
朝飯.
コーヒー.
0750 自宅発.
ものすごい風雪.
0820 研究室着.
[いきなりまっしろ]
気温はそんなに下がってないんだけど,
まあこんなかんじで.
-
EP の中で,
あやしげな server がこわれていることに気づいた.
記録のため書いとくと,
plotnet なデータを置いてた
eco1.ees.hokudai.ac.jp
.
[暗い区画におしこめといたのだが]
HDD のファイルシステム ext3 がこわれたもよう.
先日つくった USB な ubuntu 12.10 を boot して
fsck -f -y -C0 /dev/sda2
したんだけど
……
なおらん!
もうハードウェアが壊れたのか?
かなり古い機器だし
……
ここで投げだし,
甲山さんに報告しといた.
あとは野となれ山となれ.
-
1030 より研究室セミナー,
今日は北條さんで石狩浜 (この地名をみるとアタマがぐるぐるする今日このごろ)
銭函の海岸カシワ林の樹高 130 cm 以下の稚樹あれこれ
……
こいつら小さいくせに萌芽更新するので,
データ構造がすごくややこしくなる,
と判明した.
-
昼飯.
-
で,
その後はカシワな検討に
……
-
1730 研究室発.
買いもの.
1800 帰宅.
晩飯の準備のてつだい.
晩飯.
-
夜になってとら坊がわーわー (というかもっと瀕死感のある調子で)
泣くのは
たそがれ泣き
という普遍的な乳児行動みたいだね.
たそがれというより,
夜になってから泣くのは遅寝遅おきなせいなのかな?
-
[今日の運動]
-
たそがれ泣きのとらちゃんとすろー & 長時間なだんす
-
[今日の食卓]
- 朝 (0730):
米 0.7 合.
ネギ・アサリの味噌汁.
卵.
- 昼 (1200):
研究室.
米 0.7 合.
ヒジキの煮物.
リンゴ.
- 晩 (1900):
米 0.7 合.
サツマイモのポタージュ.
鶏もも肉グリル.
ニンジンのサラダ.
2012 年 11 月 28 日 (水)
-
12/17-19 に静岡・清水の
国際水研
で統計モデルせみなーみたいなことをやらないといけない
……
何も準備ができてないのだけど,
よーやく日程の「わく」みたいなものはできてきたので,
いろいろメモをかいてみる.
-
またこんさるメイル
……
尤度比検定が何をやっているのか,
いっしゅんわからなくなった.
ばてぎみなのかな
……
-
北大くらうどに置いてる web server 機に ssh
でアクセスできなくなった
……
httpd も止まっていたので,
(web なくらうど UI で)
とりあえず
sudo sevice apache2 restart
したところ,
なぜかしら ssh もつながるようになった.
メモリりーくとか?
-
Cardamine 作図やりなおしの苦闘.
-
1800 研究室発.
1830 帰宅.
晩飯の準備.
晩飯.
-
とらちゃん,
また連続ばくはつぶつ事件.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0700):
米 0.7 合.
サツマイモスープ.
ニンジンサラダ.
- 昼 (1200):
研究室.
米 0.7 合.
ニンジンサラダ.
リンゴ.
- 晩 (1940):
ダイコン・ネギ・シュンギク・ぷちがんも・豆腐・鶏肉の水炊き.
うどん.
2012 年 11 月 29 日 (木)
-
0740 起床.
不調.
あけがた 5 時ごろ,
とーとつに吐瀉してしまった.
まあ,トイレに行って吐くぐらいの余裕はあったけど.
起床後も吐く.
うーむ,
胃潰瘍でもなければ風邪でもないみたいなので,
食あたりか何かだと思うんだけど,
原因が特定できづ.
朝飯は食べても吐くだけ,
というかんじがする.
コーヒーとかも無理.
0810 自宅発.
曇天.
ちょっと気温あがったけど,
路面はあいかわらず凍結・つるつる.
0840 研究室着.
歩くとかなり「ヨワっている」のが実感される
……
何か昨日,
ヘンなもの食ったかな?
-
液体ならだいじょうぶだろうとリンゴジュースのんでいたのだが
……
これまたぜんぶ吐いてしまった.
胃はすっきりしたけれど.
-
Wikipedia の
嘔吐
の項目を読んで,
つめたい飲みものはやめて,
スポーツドリンクのお湯わりをのむ
……
お,
これはなかなか胃にいいな.
-
午前中はかくのごとくぢたばたしつつ,
データ解析こんさるメイルかき,
作図プログラミングのやりなおしなど.
昼飯
……
は食う気しないな.
スープ春雨かったので,
これだけ食べてみようかな.
-
食べてみたらやっぱりちょっとキモチわるい
……
とはいえ,
回復しつつあるよーで,
吐かずにすんだ.
-
「統計モデリング入門」
第 5 刷の
修正箇所
のとりまとめと連絡とか.
-
寝不足なので,
ねむい.
-
ちょっと早めに帰宅.
晩飯,
すこしだけ.
-
体温 37.8°C.
明日は休む.
-
[今日の運動]
-
[今日の食卓]
- 朝:
食ってない.
- 昼 (1300):
研究室.
春雨スープ.
- 晩 (1900):
米 0.4 合.
鶏スープ.
2012 年 11 月 30 日 (金)
-
昨晩は
37.7°C
だった体温,
今朝は
37.1°C
に.
8 時ごろ起床かな.
-
まだ食欲はないので,
ミルクコーヒーとマリービスケット 3 枚.
-
熱は下がり胃腸まわりも復旧しつつあるが,
まだまだへろへろなので,
とらちゃんをだっこして動きまわれない
……
ここでまた試行錯誤と創意工夫が要求される.
-
とりあえずの結果:
-
部屋をくらくして,
いっしょにふとんに入る
-
とら坊,両手をひろげてえーんと泣く
-
広げた両手それぞれを,
軽くにぎってやってあれこれあやす
(単調な語りかけとか)
-
とらちゃん,
5-10 分ぐらいで寝る
今日はこれで 3 回ほど寝かしつけに成功.
もちろん,
空腹だったりすると,
ぜんぜん無視されてしまって,
泣きつづける.
-
夜のたそがれ泣きにも通用しなかった.
Beco によるだっこで泣きやんだ.
-
昼飯も食う気がおきない.
かなり回復してると思うんだけど
……
しかし,
ビスケット食ってもかなり胃がもたれたので,
胃がヨワっているんでしょう.
-
午後もとらちゃんのあいてを.
-
[今日の運動]
-
とらちゃんえくささいず,
少しだけ.
しかしながら,
かなりばてた.
-
[今日の食卓]
- 朝 (0900):
ビスケット 3 枚.
ミカン.
- 昼 (1300):
ビスケット 3 枚.
ミカン.
- 晩 (1900):
ビスケット 3 枚.
ミカン.