我正在尝试对我正在使用的大量glm模型进行一些后处理,但我需要从glm对象中提取有关数据子集的信息。
作为一个玩具示例:
x <- rnorm(100)
y <- rnorm(100,x,0.5)
s<-sample(c(T,F),100,replace=T)
myGlm <- glm(y~x, subset= s)
由此,我需要知道从myGlm
中获取信息所使用的100个观察中的哪一个。我认为myGlm$data
会有子集数据,但它实际上有100个观察值。我查看str(myGlm)
无济于事。但是,很明显,在对象的某处,存储了有关子集s
的信息。
这似乎应该是完全无足轻重的!
提前致谢。
答案 0 :(得分:1)
as.numeric(rownames(myGlm$model))