2012-10-04 3 views
1

Je cherche un moyen rapide de tirer le nombre d'observations (N) généré dans le résumé d'un test lme, au lieu de copier et coller? J'ai trouvé this en ligne mais la commande ne semble pas être intégrée et le paquet texreg n'est pas téléchargé sur ce site.extraire N de lme résumé

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

Merci beaucoup.

+0

'bibliothèque (nlme); fm1 <- lme (distance ~ âge, données = orthodonte); nobs (fm1); nobs (summary (fm1)) 'fonctionne très bien pour moi ... –

+0

Je reçois toujours la même erreur avec nobs:: Impossible de trouver la fonction" nobs ". Je cours v 2.12.1. Pensez-vous que c'est une version? – AWC

+3

Une raison de ne pas mettre à jour? – Roland

Répondre

2
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 
+0

fit $ dims $ N- merci! Fonctionne avec brio. – AWC

+0

@AWC N'hésitez pas à accepter la réponse pour marquer la réponse à la question. – Roland

Questions connexes