r - 从 lme 摘要中提取 N

Question

我正在寻找一种快速的方法来提取测试N摘要中生成的观察数 () lme，而不是复制和粘贴？我在网上找到了这个，但该命令似乎不是内置的，并且texreg包没有上传到这个站点。

n <- nobs(model) #extract number of observations

非常感谢。

Answer 1

library(nlme)

x <- 1:80
g <- sort(rep(c(1,2,3,4),20))
set.seed(1)
y <- 3*x+g*0.01+rnorm(80,sd=0.01)

plot(y~x)
fit <- lme(y~x,random=~1|g)

nobs(fit)
[1] 80 #works just fine

fit$dims$N
[1] 80 #alternative