从 lme 摘要中提取 N答案

【问题标题】：extracting N from lme summary从 lme 摘要中提取 N
【发布时间】：2012-10-04 18:07:51
【问题描述】：

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

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

非常感谢。

【问题讨论】：

library(nlme); fm1 <- lme(distance ~ age, data = Orthodont); nobs(fm1); nobs(summary(fm1)) 对我来说很好......
我仍然遇到与 nobs 相同的错误：：找不到函数“nobs”。我正在运行 v 2.12.1。你认为这是一个版本的事情吗？
以前没问题的，呵呵。
@AWC 请注意，即使安装软件包也会变得困难，因为 R 的内部机制在 2.14 和 2.15 版本中发生了相当大的变化，更不用说由于编译了许多基本函数当前的 R 比您正在使用的恐龙快得多。

标签： r summary mixed-models

【解决方案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

【讨论】：