「統計学関連なんでもあり」の過去ログ--- 041

No.05536　glmmLMについて　　【Rはじめました】　2008/01/25(Fri) 18:50

いつもお世話になっております。

A <- c(19,4,9,2,32,2,3,0,20,22,4,2,8,7,2,3,10,5,4,4,7,2,5,2)
B <- c(rep(1,12),rep(2,12))
C <- rep(seq(1,4),6)
D <- seq(1,24)
res <- data.frame(A,B,C,D)

というデータに対して，BとCがDに与える影響の有無を知りたいときに，

library(glmmML)
ex <- glmmML(D~B+C,cluster=A,family=poisson,data=res)
step(ex)

とすると，B，C両方の効果を含むモデルと，B,Cの効果いずれか一方を考慮しない場合のAICが求まりますが，両方とも考慮しない場合のAICは出てきません。
BもCもDに作用していない場合のモデルとの比較はどのようにすればよいのでしょうか。
例えば

ex1 <- glmmML(D~1,cluster=A,family=poisson,data=res)

のようにしてもエラーになります。

No.05546　Re: glmmLMについて　　【ごう】　2008/01/26(Sat) 18:56
ex.sigma <- glmmML(D~1,cluster=A,family=poisson,data=res,start.sigma=10)
summary(ex.sigma)
http://www.google.com/search?hl=ja&ie=UTF-8&domains=hosho.ees.hokudai.ac.jp&sitesearch=hosho.ees.hokudai.ac.jp&q=start.sigma

参考になりますかどうか。

No.05550　Re: glmmLMについて　　【Rはじめました】　2008/01/26(Sat) 22:14

ごう様，的確なアドバイスをいただき，誠にありがとうございました。
大変参考になり，よくわかりました。