KuboLog 2006-05-17

苦情・お叱りは，たいへんお手数かけて恐縮ですが，久保 (kubo@ees.hokudai.ac.jp) までお知らせください．

2006 年 05 月 17 日 (水)

0710 起床．朝飯．コーヒー． 0900 自宅発．晴．昨日あたりからいきなり気温高い． 0915 研究室着．

本日もぢりぢり GLMM 解説原稿かき．

Perl-5.8.2 の中での日本語文字あつかいの問題に関して …… 霜月さんご指摘のように use Encode; してみる．

#!/usr/bin/perl -w
use strict;
use Encode;
my $file = shift @ARGV;
my $lines;
open IN, '<:encoding(euc-jp)', $file or die;
while (<IN>) {
	s/．/。/g;
	... (中略) ...
	$lines .= $_;
}
close IN;
open OUT, '>:encoding(euc-jp)', $file or die;
print OUT $lines;
close OUT;
__END__

さらに阿辺川さんの解説によると， use Encode; のかわりに (たとえば) use encoding 'euc-jp'; としておけば Perl スクリプトじたいが (この例の場合だと) EUC-JP で書かれていても動作する，とのこと …… 実験してみたらたしかにそうなった． Perl 内部のあつかいは UTF8 なので， EUC-JP で読みこんで書き出すならばこの場合でも入出力どちらの open にも encoding(euc-jp) うんぬんが必要である．

ということで，私などもよーやくにしてケガれ言語 Perl にて日本語文字がうまくあつかえるようになりました．めでたしめでたし．

日本語といえば …… R の最新ヴァージョン 2.3.0 は中間栄治さんのコードが取り込まれたとのことなんで日本語あつかいまわりはさらによくなっているはず．実験してみると png() device でも問題なく日本語がでる (ひょっとして以前から問題なかったっけ?)

昼飯．

夜まで GLMM 解説原稿かき．

「overdispersion の発見」について．世間では (residual deviance) / (degree of freedom) なる割り算値が 1 から逸脱だのどーのこーので overdispersion が見つかる，と信じられているわけだが …… まあ，なんというかこれはいかがわしいかんぢがする．

そもそもこれで何とかなる，と考えてる連中は「なぜそういう割り算値を使うのか?」を理解してるのだろうか? …… これは例によって近似の近似といった世界で，尤度比検定とおなじよーなかんぢで

-2 log[(full モデル尤度) / (いまのモデルの尤度)]

なる数量 (residual deviance) が標本数おおくなると自由度 (差) n の χ² 分布 (その平均値は n) に漸近するはず，といったハナシだ．これってホントかよ，というような．

で，まあ簡単に図にかける場合には (次元の高い場合は Table でも作りますかね?)，むしろ「推定されたモデルなら標本のばらつきはこれぐらいになるはず」といった範囲を示したほうがマシではないか …… といった説明にしようかと検討中．

[二項分布モデルから期待される]

標本のばらつきの範囲．新潟大会自由集会 GLMM 説明の解説用例題と対応している図．上下の破線にはさまれた範囲が 99% 区間．まあこういうふうに図で示せる場合なら， overdispersion あることは一目瞭然，ですよね ……

1910 研究室発． 1930 帰宅．運動．晩飯．

[今日の運動]

エアロバイク 50 分間．

[今日の食卓]

朝 (0810): 米麦 0.6 合．キャベツ・タマネギ・シイタケ・豆腐のカレー．
昼 (1310): 研究室お茶部屋．米麦 0.5 合．キャベツ・タマネギ・シイタケ・豆腐のカレー．
晩 (2140): 米麦 0.7 合．ニラ・マイタケ・卵炒飯．コマツナあえもの．

ぎょーむ日誌 2006-05-17

2006 年 05 月 17 日 (水)