昨日の続き。

R 統計

でも！あまり好まれていないglmmPQLではできた！idを50にしても。

> fglmmPQL<-glmmPQL(y2~1,random=~1|id,family=gaussian,data=dat2)
iteration 1
iteration 2
> fglmmPQL
Linear mixed-effects model fit by maximum likelihood
  Data: dat2 
  Log-likelihood: NA
  Fixed: y2 ~ 1 
(Intercept) 
  0.9502553 

Random effects:
 Formula: ~1 | id
        (Intercept) Residual
StdDev:    3.735389 1.400771

Variance function:
 Structure: fixed weights
 Formula: ~invwt 
Number of Observations: 50
Number of Groups: 50

推定結果は微妙ではあるが。
このできるできないの違いは何なんだろう。
いやはや，統計ソフトにぶち込めばイイということがまかり通るようになってしまったがタメのツケです。

電子レンジばっか使っていると料理の腕がなまる。素材を生かし，おいしく，より美しくみせるには電子レンジを使ってもいいけど電子レンジの使いどころを見極め，料理人ならではの工夫が必要です。

電子レンジ分解するのなかなか大変ですよ。