從lme摘錄N摘要

Question

我正在尋找一種快速方法來提取在lme測試摘要中生成的觀察值數量（N），而不是複製和粘貼？我發現this在線，但該命令似乎不是內置的，並且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