1

次のようなデータセットの各グループの切片の信頼区間を推定しました。

d <- data.frame(
    g=sample(letters[1:5], 100, replace=TRUE),
    y=runif(100)
)

library(lme4)

fm <- lmList(y ~ 1 | g, data=d)
ci.fm <- confint(fm)
ci.fm
# An object of class "lmList.confint"
# , , (Intercept)
# 
#       2.5 %    97.5 %
# a 0.4253236 0.6668585
# b 0.4860047 0.6794625
# c 0.4071186 0.6822713
# d 0.2785848 0.5992378
# e 0.4110070 0.6890399

次に、各列に個別にインデックスを付ける方法が必要ですが、私にはわかりません。私は他の同様のものを見て、試しましたが、すべてうまくstr(ci.fm)いきませんでした。正しい方向へのポインタに感謝します。ありがとう。attributes(ci.fm)attr(ci.fm, dimnames[[2]]$"2.5 %")

4

2 に答える 2

2

これはS4オブジェクトです。で「スロット」のリストを作成し、 （ではなく）でslotNamesそれらにアクセスできます。この場合、唯一のスロットは、先頭にドットが付いているため、実際のデータが含まれています。@$.Dataci.fm@.Data

str(ci.fm)
# Formal class 'lmList.confint' [package "lme4"] with 1 slots
#   ..@ .Data: num [1:5, 1:2, 1] 0.35 0.348 0.315 0.143 0.388 ...
#   .. ..- attr(*, "dimnames")=List of 3
#   .. .. ..$ : chr [1:5] "a" "b" "c" "d" ...
#   .. .. ..$ : chr [1:2] "2.5 %" "97.5 %"
#   .. .. ..$ : chr "(Intercept)"

slotNames(ci.fm)
# [1] ".Data"

str( ci.fm@.Data )  # 3-dimensional array
#  num [1:5, 1:2, 1] 0.35 0.348 0.315 0.143 0.388 ...
#  - attr(*, "dimnames")=List of 3
#   ..$ : chr [1:5] "a" "b" "c" "d" ...
#   ..$ : chr [1:2] "2.5 %" "97.5 %"
#   ..$ : chr "(Intercept)"

ci.fm@.Data[,1,1]  # First column
#         a         b         c         d         e 
# 0.3496147 0.3482908 0.3147435 0.1434213 0.3876549
于 2012-06-28T15:27:12.300 に答える
2 
ci.fm[,"2.5 %",]
        a         b         c         d         e 
0.4463712 0.3297468 0.3622782 0.2901415 0.2847008
于 2012-06-28T15:28:19.383 に答える