body(glm)
などとすれば関数の中身は調べられる.
しかしそれではうまくいかない場合もあり,
たとえば
library(MASS)
の AIC()
関数とかオブジェクト指向なやつはダメ.
その場合には,
methotds("AIC")
とすれば "AIC.default"
といったメソッドがあることがわかるので
(これは S3 オブジェクト指向限定の方法だけど),
getS3method("AIC", "default")
などとすればコーディングがわかる.
glm(..., family = "gaussian")
の dispersion parameter は
(weight
指定でなければ)
sum(fit$residuals^2) / fit$df.residual
と評価されていて,
これはつまり不偏分散のこと.