library(gnlm)
があったんで，
これでいいんじゃないでしょうか，
と返信してしまった．
で，
今日になって気づいたんだけど，
これって CRAN pacakge じゃなくて，
Lindsey さん
の個人サイトにおいてあるやつだな．
library(gnm)
で
example(gnm)
したら問題なく終了したんだけど，
Warning message がでた
(この example の結果とは無関係)．
perl=TRUE は UTF-8 ロケールに対してのみ完全実装されています
……
うーむ R 内部の UTF-8 化ぷろぢぇくとが進行してますな．
help(gnlr)
だけではその挙動はよくわからず，
gnlr()
関数の定義をよく読む必要がある．
たとえば，
負の二項分布は以下のように定義されており，
尤度を評価している pnbinom()
の関数の使われかたに注意する必要がある．
`negative binomial` = {
fcn <- function(p) {
m <- mu1(p)
t <- sh1(p)
s <- exp(t)
-sum(wt * (cc * (lgamma(y[, 1] + s) - lgamma(s) +
s * t + y[, 1] * log(m) - (y[, 1] + s) * log(s +
m)) + log(lc - rc * pnbinom(y[, 1], s, 1/(1 + m/s))
)))
}
const <- wt * cc * lgamma(y[, 1] + 1)
}
メールしてんだけど。無視かい？本社セイコーマートの人とは違うのかい？そろそろ「不安」になってきたようだね． いーかげん北大にメイルだしてることに気づけよ．
（2009年度 正会員学会費：5,000円） 郵便振替： 00190-2-173844 応用統計学会 5,000円 もしくは、 みずほ銀行 大塚支店 普通口座 1281134 応用統計学会 5,000円会費はネット経由でぽいっとふりこむ．
make pdf
して，
rsync
による同期を命じて，
共立出版にメイルかいてみる．