パッケージlme4の関数lmer()などの結果のオブジェクト

2009.1.25最終修正
（2009.1.22より　）

　関数lmer()などの結果オブジェクトのクラスはmer、そのsummary()の結果オブジェクトのクラスはsummary.merである。いずれもS4オブジェクトである。

　例として使ったのは、こちらのもので、lmer()の結果オブジェクトがresm13.lmer、sumresm13.lmerはsummary(resm13.lmer)である。

　まず、summary() の結果の方から。。オブジェクト名を入力すると

> sumresm13.lmer
Linear mixed model fit by REML
Formula: yy1 ~ 1 + (1 | xf11)
Data: mydata13
AIC BIC logLik deviance REMLdev
-12.55 -10.24 9.277 -20.74 -18.55
Random effects:
Groups Name Variance Std.Dev.
xf11 (Intercept) 0.0723802 0.269036
Residual 0.0066021 0.081253
Number of obs: 16, groups: xf11, 4

Fixed effects:
Estimate Std. Error t value
(Intercept) 3.166 0.136 23.27
と表示される。

　slotNames(）でスロットを確認する。以下のものがある、
> slotNames(sumresm13.lmer)
[1] "methTitle" "logLik" "ngrps" "sigma" "coefs" "vcov"
[7] "REmat" "AICtab" "env" "nlmodel" "frame" "call"
[13] "flist" "X" "Xst" "Zt" "pWt" "offset"
[19] "y" "Gp" "dims" "ST" "V" "A"
[25] "Cm" "Cx" "L" "deviance" "fixef" "ranef"
[31] "u" "eta" "mu" "muEta" "var" "resid"
[37] "sqrtXWt" "sqrtrWt" "RZX" "RX" "ghx" "ghw"

いくつかのスロットを見てみる。

　logLikは対数尤度、
> sumresm13.lmer@logLik
'log Lik.' 9.277319 (df=3)

　AICtabはAICなど、
> sumresm13.lmer@AICtab
AIC BIC logLik deviance REMLdev
-12.55464 -10.23687 9.277319 -20.73898 -18.55464

　devianceはデビアンスなど、
> sumresm13.lmer@deviance
ML REML ldL2 ldRX2 sigmaML
-20.73898394 -18.55463751 15.21336858 -1.03075342 0.07867308
sigmaREML pwrss disc usqr wrss
0.08125321 0.09903125 0.07966662 0.01936556 0.07966569
dev llik NULLdev
NA NA NA

関数str()で見ると、42個のスロットがある
> str(sumresm13.lmer)
Formal class 'summary.mer' [package "lme4"] with 42 slots
..@ methTitle: chr "Linear mixed model fit by REML"
..@ logLik :Class 'logLik' : 9.3 (df=3)
..@ ngrps : Named int 4
.. ..- attr(*, "names")= chr "xf11"
..@ sigma : Named num 0.0813
.. ..- attr(*, "names")= chr "sigmaREML"
..@ coefs : num [1, 1:3] 3.166 0.136 23.27
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr "(Intercept)"
.. .. ..$ : chr [1:3] "Estimate" "Std. Error" "t value"
..@ vcov :Formal class 'dpoMatrix' [package "Matrix"] with 5 slots
.. .. ..@ x : num 0.0185
.. .. ..@ Dim : int [1:2] 1 1
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ uplo : chr "U"
.. .. ..@ factors :List of 2
.. .. .. ..$ Cholesky :Formal class 'Cholesky' [package "Matrix"] with 5 slots
.. .. .. .. .. ..@ x : num 0.136
.. .. .. .. .. ..@ Dim : int [1:2] 1 1
.. .. .. .. .. ..@ Dimnames:List of 2
.. .. .. .. .. .. ..$ : NULL
.. .. .. .. .. .. ..$ : NULL
.. .. .. .. .. ..@ uplo : chr "U"
.. .. .. .. .. ..@ diag : chr "N"
.. .. .. ..$ correlation:Formal class 'corMatrix' [package "Matrix"] with 6 slots
.. .. .. .. .. ..@ sd : num 0.136
.. .. .. .. .. ..@ x : num 1
.. .. .. .. .. ..@ Dim : int [1:2] 1 1
.. .. .. .. .. ..@ Dimnames:List of 2
.. .. .. .. .. .. ..$ : NULL
.. .. .. .. .. .. ..$ : NULL
.. .. .. .. .. ..@ uplo : chr "U"
.. .. .. .. .. ..@ factors :List of 1
.. .. .. .. .. .. ..$ Cholesky:Formal class 'Cholesky' [package "Matrix"] with 5 slots
.. .. .. .. .. .. .. .. ..@ x : num 1
.. .. .. .. .. .. .. .. ..@ Dim : int [1:2] 1 1
.. .. .. .. .. .. .. .. ..@ Dimnames:List of 2
.. .. .. .. .. .. .. .. .. ..$ : NULL
.. .. .. .. .. .. .. .. .. ..$ : NULL
.. .. .. .. .. .. .. .. ..@ uplo : chr "U"
.. .. .. .. .. .. .. .. ..@ diag : chr "N"
..@ REmat : chr [1:2, 1:4] "xf11" "Residual" "(Intercept)" "" ...
.. ..- attr(*, "dimnames")=List of 2
.. .. ..$ : chr [1:2] "" ""
.. .. ..$ : chr [1:4] "Groups" "Name" "Variance" "Std.Dev."
..@ AICtab :'data.frame': 1 obs. of 5 variables:
.. ..$ AIC : num -12.6
.. ..$ BIC : num -10.2
.. ..$ logLik : num 9.28
.. ..$ deviance: num -20.7
.. ..$ REMLdev : num -18.6
..@ env :<environment: 0x0ba1c010>
..@ nlmodel : language I(x)
..@ frame :'data.frame': 16 obs. of 2 variables:
.. ..$ yy1 : num [1:16] 3.28 3.09 3.03 3.03 3.52 3.48 3.38 3.38 2.88 2.8 ...
.. ..$ xf11: Factor w/ 4 levels "1","2","3","4": 1 1 1 1 2 2 2 2 3 3 ...
.. ..- attr(*, "terms")=Classes 'terms', 'formula' length 3 yy1 ~ 1
.. .. .. ..- attr(*, "variables")= language list(yy1)
.. .. .. ..- attr(*, "factors")= int(0)
.. .. .. ..- attr(*, "term.labels")= chr(0)
.. .. .. ..- attr(*, "order")= int(0)
.. .. .. ..- attr(*, "intercept")= int 1
.. .. .. ..- attr(*, "response")= int 1
.. .. .. ..- attr(*, ".Environment")=<environment: R_GlobalEnv>
.. .. .. ..- attr(*, "predvars")= language list(yy1)
.. .. .. ..- attr(*, "dataClasses")= Named chr "numeric"
.. .. .. .. ..- attr(*, "names")= chr "yy1"
..@ call : language lmer(formula = yy1 ~ 1 + (1 | xf11), data = mydata13)
..@ flist :'data.frame': 16 obs. of 1 variable:
.. ..$ xf11: Factor w/ 4 levels "1","2","3","4": 1 1 1 1 2 2 2 2 3 3 ...
.. ..- attr(*, "assign")= int 1
..@ X : num [1:16, 1] 1 1 1 1 1 1 1 1 1 1 ...
.. ..- attr(*, "assign")= int 0
..@ Xst :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int(0)
.. .. ..@ p : int 0
.. .. ..@ Dim : int [1:2] 0 0
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num(0)
.. .. ..@ factors : list()
..@ Zt :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int [1:16] 0 0 0 0 1 1 1 1 2 2 ...
.. .. ..@ p : int [1:17] 0 1 2 3 4 5 6 7 8 9 ...
.. .. ..@ Dim : int [1:2] 4 16
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num [1:16] 1 1 1 1 1 1 1 1 1 1 ...
.. .. ..@ factors : list()
..@ pWt : num(0)
..@ offset : num(0)
..@ y : num [1:16] 3.28 3.09 3.03 3.03 3.52 3.48 3.38 3.38 2.88 2.8 ...
..@ Gp : int [1:2] 0 4
..@ dims : Named int [1:18] 1 16 1 4 1 1 0 1 2 5 ...
.. ..- attr(*, "names")= chr [1:18] "nt" "n" "p" "q" ...
..@ ST :List of 1
.. ..$ : num [1, 1] 3.31
.. .. ..- attr(*, "dimnames")=List of 2
.. .. .. ..$ : chr "(Intercept)"
.. .. .. ..$ : chr "(Intercept)"
..@ V : num[0 , 0 ]
..@ A :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int [1:16] 0 0 0 0 1 1 1 1 2 2 ...
.. .. ..@ p : int [1:17] 0 1 2 3 4 5 6 7 8 9 ...
.. .. ..@ Dim : int [1:2] 4 16
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : chr [1:4] "1" "2" "3" "4"
.. .. .. ..$ : NULL
.. .. ..@ x : num [1:16] 3.31 3.31 3.31 3.31 3.31 ...
.. .. ..@ factors : list()
..@ Cm :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int(0)
.. .. ..@ p : int 0
.. .. ..@ Dim : int [1:2] 0 0
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num(0)
.. .. ..@ factors : list()
..@ Cx : num(0)
..@ L :Formal class 'dCHMsimpl' [package "Matrix"] with 10 slots
.. .. ..@ x : num [1:4] 6.7 6.7 6.7 6.7
.. .. ..@ p : int [1:5] 0 1 2 3 4
.. .. ..@ i : int [1:4] 0 1 2 3
.. .. ..@ nz : int [1:4] 1 1 1 1
.. .. ..@ nxt : int [1:6] 1 2 3 4 -1 0
.. .. ..@ prv : int [1:6] 5 0 1 2 3 -1
.. .. ..@ colcount: int [1:4] 1 1 1 1
.. .. ..@ perm : int [1:4] 0 1 2 3
.. .. ..@ type : int [1:4] 2 1 0 1
.. .. ..@ Dim : int [1:2] 4 4
..@ deviance : Named num [1:13] -20.739 -18.5546 15.2134 -1.0308 0.0787 ...
.. ..- attr(*, "names")= chr [1:13] "ML" "REML" "ldL2" "ldRX2" ...
..@ fixef : Named num 3.17
.. ..- attr(*, "names")= chr "(Intercept)"
..@ ranef : num [1:4] -0.0568 0.2683 -0.3453 0.1338
..@ u : num [1:4] -0.0172 0.081 -0.1043 0.0404
..@ eta : num [1:16] 3.11 3.11 3.11 3.11 3.43 ...
..@ mu : num [1:16] 3.11 3.11 3.11 3.11 3.43 ...
..@ muEta : num(0)
..@ var : num(0)
..@ resid : num [1:16] 0.1712 -0.0188 -0.0788 -0.0788 0.0861 ...
..@ sqrtXWt : num[0 , 1]
..@ sqrtrWt : num(0)
..@ RZX : num [1:4, 1] 1.98 1.98 1.98 1.98
..@ RX : num [1, 1] 0.597
..@ ghx : num(0)
..@ ghw : num(0)

　次にlmerの結果オブジェクト（ここの例ではresm13.lmer）。

> resm13.lmer
Linear mixed model fit by REML
Formula: yy1 ~ 1 + (1 | xf11)
Data: mydata13
AIC BIC logLik deviance REMLdev
-12.55 -10.24 9.277 -20.74 -18.55
Random effects:
Groups Name Variance Std.Dev.
xf11 (Intercept) 0.0723802 0.269036
Residual 0.0066021 0.081253
Number of obs: 16, groups: xf11, 4

Fixed effects:
Estimate Std. Error t value
(Intercept) 3.166 0.136 23.27
上のsummary()の結果と比べるとそっくりである。

　slotNames(）でスロットを確認する。

> slotNames(resm13.lmer)
[1] "env" "nlmodel" "frame" "call" "flist" "X"
[7] "Xst" "Zt" "pWt" "offset" "y" "Gp"
[13] "dims" "ST" "V" "A" "Cm" "Cx"
[19] "L" "deviance" "fixef" "ranef" "u" "eta"
[25] "mu" "muEta" "var" "resid" "sqrtXWt" "sqrtrWt"
[31] "RZX" "RX" "ghx" "ghw"
> str(resm13.lmer)
Formal class 'mer' [package "lme4"] with 34 slots
..@ env :<environment: 0x0ba1c010>
..@ nlmodel : language I(x)
..@ frame :'data.frame': 16 obs. of 2 variables:
.. ..$ yy1 : num [1:16] 3.28 3.09 3.03 3.03 3.52 3.48 3.38 3.38 2.88 2.8 ...
.. ..$ xf11: Factor w/ 4 levels "1","2","3","4": 1 1 1 1 2 2 2 2 3 3 ...
.. ..- attr(*, "terms")=Classes 'terms', 'formula' length 3 yy1 ~ 1
.. .. .. ..- attr(*, "variables")= language list(yy1)
.. .. .. ..- attr(*, "factors")= int(0)
.. .. .. ..- attr(*, "term.labels")= chr(0)
.. .. .. ..- attr(*, "order")= int(0)
.. .. .. ..- attr(*, "intercept")= int 1
.. .. .. ..- attr(*, "response")= int 1
.. .. .. ..- attr(*, ".Environment")=<environment: R_GlobalEnv>
.. .. .. ..- attr(*, "predvars")= language list(yy1)
.. .. .. ..- attr(*, "dataClasses")= Named chr "numeric"
.. .. .. .. ..- attr(*, "names")= chr "yy1"
..@ call : language lmer(formula = yy1 ~ 1 + (1 | xf11), data = mydata13)
..@ flist :'data.frame': 16 obs. of 1 variable:
.. ..$ xf11: Factor w/ 4 levels "1","2","3","4": 1 1 1 1 2 2 2 2 3 3 ...
.. ..- attr(*, "assign")= int 1
..@ X : num [1:16, 1] 1 1 1 1 1 1 1 1 1 1 ...
.. ..- attr(*, "assign")= int 0
..@ Xst :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int(0)
.. .. ..@ p : int 0
.. .. ..@ Dim : int [1:2] 0 0
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num(0)
.. .. ..@ factors : list()
..@ Zt :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int [1:16] 0 0 0 0 1 1 1 1 2 2 ...
.. .. ..@ p : int [1:17] 0 1 2 3 4 5 6 7 8 9 ...
.. .. ..@ Dim : int [1:2] 4 16
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num [1:16] 1 1 1 1 1 1 1 1 1 1 ...
.. .. ..@ factors : list()
..@ pWt : num(0)
..@ offset : num(0)
..@ y : num [1:16] 3.28 3.09 3.03 3.03 3.52 3.48 3.38 3.38 2.88 2.8 ...
..@ Gp : int [1:2] 0 4
..@ dims : Named int [1:18] 1 16 1 4 1 1 0 1 2 5 ...
.. ..- attr(*, "names")= chr [1:18] "nt" "n" "p" "q" ...
..@ ST :List of 1
.. ..$ : num [1, 1] 3.31
.. .. ..- attr(*, "dimnames")=List of 2
.. .. .. ..$ : chr "(Intercept)"
.. .. .. ..$ : chr "(Intercept)"
..@ V : num[0 , 0 ]
..@ A :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int [1:16] 0 0 0 0 1 1 1 1 2 2 ...
.. .. ..@ p : int [1:17] 0 1 2 3 4 5 6 7 8 9 ...
.. .. ..@ Dim : int [1:2] 4 16
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : chr [1:4] "1" "2" "3" "4"
.. .. .. ..$ : NULL
.. .. ..@ x : num [1:16] 3.31 3.31 3.31 3.31 3.31 ...
.. .. ..@ factors : list()
..@ Cm :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int(0)
.. .. ..@ p : int 0
.. .. ..@ Dim : int [1:2] 0 0
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num(0)
.. .. ..@ factors : list()
..@ Cx : num(0)
..@ L :Formal class 'dCHMsimpl' [package "Matrix"] with 10 slots
.. .. ..@ x : num [1:4] 6.7 6.7 6.7 6.7
.. .. ..@ p : int [1:5] 0 1 2 3 4
.. .. ..@ i : int [1:4] 0 1 2 3
.. .. ..@ nz : int [1:4] 1 1 1 1
.. .. ..@ nxt : int [1:6] 1 2 3 4 -1 0
.. .. ..@ prv : int [1:6] 5 0 1 2 3 -1
.. .. ..@ colcount: int [1:4] 1 1 1 1
.. .. ..@ perm : int [1:4] 0 1 2 3
.. .. ..@ type : int [1:4] 2 1 0 1
.. .. ..@ Dim : int [1:2] 4 4
..@ deviance: Named num [1:13] -20.739 -18.5546 15.2134 -1.0308 0.0787 ...
.. ..- attr(*, "names")= chr [1:13] "ML" "REML" "ldL2" "ldRX2" ...
..@ fixef : Named num 3.17
.. ..- attr(*, "names")= chr "(Intercept)"
..@ ranef : num [1:4] -0.0568 0.2683 -0.3453 0.1338
..@ u : num [1:4] -0.0172 0.081 -0.1043 0.0404
..@ eta : num [1:16] 3.11 3.11 3.11 3.11 3.43 ...
..@ mu : num [1:16] 3.11 3.11 3.11 3.11 3.43 ...
..@ muEta : num(0)
..@ var : num(0)
..@ resid : num [1:16] 0.1712 -0.0188 -0.0788 -0.0788 0.0861 ...
..@ sqrtXWt : num[0 , 1]
..@ sqrtrWt : num(0)
..@ RZX : num [1:4, 1] 1.98 1.98 1.98 1.98
..@ RX : num [1, 1] 0.597
..@ ghx : num(0)
..@ ghw : num(0)

　str()で中味のつくりを見る、

> str(resm13.lmer)
Formal class 'mer' [package "lme4"] with 34 slots
..@ env :<environment: 0x0ba1c010>
..@ nlmodel : language I(x)
..@ frame :'data.frame': 16 obs. of 2 variables:
.. ..$ yy1 : num [1:16] 3.28 3.09 3.03 3.03 3.52 3.48 3.38 3.38 2.88 2.8 ...
.. ..$ xf11: Factor w/ 4 levels "1","2","3","4": 1 1 1 1 2 2 2 2 3 3 ...
.. ..- attr(*, "terms")=Classes 'terms', 'formula' length 3 yy1 ~ 1
.. .. .. ..- attr(*, "variables")= language list(yy1)
.. .. .. ..- attr(*, "factors")= int(0)
.. .. .. ..- attr(*, "term.labels")= chr(0)
.. .. .. ..- attr(*, "order")= int(0)
.. .. .. ..- attr(*, "intercept")= int 1
.. .. .. ..- attr(*, "response")= int 1
.. .. .. ..- attr(*, ".Environment")=<environment: R_GlobalEnv>
.. .. .. ..- attr(*, "predvars")= language list(yy1)
.. .. .. ..- attr(*, "dataClasses")= Named chr "numeric"
.. .. .. .. ..- attr(*, "names")= chr "yy1"
..@ call : language lmer(formula = yy1 ~ 1 + (1 | xf11), data = mydata13)
..@ flist :'data.frame': 16 obs. of 1 variable:
.. ..$ xf11: Factor w/ 4 levels "1","2","3","4": 1 1 1 1 2 2 2 2 3 3 ...
.. ..- attr(*, "assign")= int 1
..@ X : num [1:16, 1] 1 1 1 1 1 1 1 1 1 1 ...
.. ..- attr(*, "assign")= int 0
..@ Xst :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int(0)
.. .. ..@ p : int 0
.. .. ..@ Dim : int [1:2] 0 0
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num(0)
.. .. ..@ factors : list()
..@ Zt :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int [1:16] 0 0 0 0 1 1 1 1 2 2 ...
.. .. ..@ p : int [1:17] 0 1 2 3 4 5 6 7 8 9 ...
.. .. ..@ Dim : int [1:2] 4 16
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num [1:16] 1 1 1 1 1 1 1 1 1 1 ...
.. .. ..@ factors : list()
..@ pWt : num(0)
..@ offset : num(0)
..@ y : num [1:16] 3.28 3.09 3.03 3.03 3.52 3.48 3.38 3.38 2.88 2.8 ...
..@ Gp : int [1:2] 0 4
..@ dims : Named int [1:18] 1 16 1 4 1 1 0 1 2 5 ...
.. ..- attr(*, "names")= chr [1:18] "nt" "n" "p" "q" ...
..@ ST :List of 1
.. ..$ : num [1, 1] 3.31
.. .. ..- attr(*, "dimnames")=List of 2
.. .. .. ..$ : chr "(Intercept)"
.. .. .. ..$ : chr "(Intercept)"
..@ V : num[0 , 0 ]
..@ A :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int [1:16] 0 0 0 0 1 1 1 1 2 2 ...
.. .. ..@ p : int [1:17] 0 1 2 3 4 5 6 7 8 9 ...
.. .. ..@ Dim : int [1:2] 4 16
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : chr [1:4] "1" "2" "3" "4"
.. .. .. ..$ : NULL
.. .. ..@ x : num [1:16] 3.31 3.31 3.31 3.31 3.31 ...
.. .. ..@ factors : list()
..@ Cm :Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
.. .. ..@ i : int(0)
.. .. ..@ p : int 0
.. .. ..@ Dim : int [1:2] 0 0
.. .. ..@ Dimnames:List of 2
.. .. .. ..$ : NULL
.. .. .. ..$ : NULL
.. .. ..@ x : num(0)
.. .. ..@ factors : list()
..@ Cx : num(0)
..@ L :Formal class 'dCHMsimpl' [package "Matrix"] with 10 slots
.. .. ..@ x : num [1:4] 6.7 6.7 6.7 6.7
.. .. ..@ p : int [1:5] 0 1 2 3 4
.. .. ..@ i : int [1:4] 0 1 2 3
.. .. ..@ nz : int [1:4] 1 1 1 1
.. .. ..@ nxt : int [1:6] 1 2 3 4 -1 0
.. .. ..@ prv : int [1:6] 5 0 1 2 3 -1
.. .. ..@ colcount: int [1:4] 1 1 1 1
.. .. ..@ perm : int [1:4] 0 1 2 3
.. .. ..@ type : int [1:4] 2 1 0 1
.. .. ..@ Dim : int [1:2] 4 4
..@ deviance: Named num [1:13] -20.739 -18.5546 15.2134 -1.0308 0.0787 ...
.. ..- attr(*, "names")= chr [1:13] "ML" "REML" "ldL2" "ldRX2" ...
..@ fixef : Named num 3.17
.. ..- attr(*, "names")= chr "(Intercept)"
..@ ranef : num [1:4] -0.0568 0.2683 -0.3453 0.1338
..@ u : num [1:4] -0.0172 0.081 -0.1043 0.0404
..@ eta : num [1:16] 3.11 3.11 3.11 3.11 3.43 ...
..@ mu : num [1:16] 3.11 3.11 3.11 3.11 3.43 ...
..@ muEta : num(0)
..@ var : num(0)
..@ resid : num [1:16] 0.1712 -0.0188 -0.0788 -0.0788 0.0861 ...
..@ sqrtXWt : num[0 , 1]
..@ sqrtrWt : num(0)
..@ RZX : num [1:4, 1] 1.98 1.98 1.98 1.98
..@ RX : num [1, 1] 0.597
..@ ghx : num(0)
..@ ghw : num(0)
となる。

@ envと @ nlmodeは、この例では使われていないーこの２つのスロットはnlmer()の結果のときに使われる。l