
Benilton, Usando o R com R --vanilla consegui usar a lme4, veja
require(lme4) Carregando pacotes exigidos: lme4 Carregando pacotes exigidos: Matrix Carregando pacotes exigidos: lattice
Anexando pacote: 'Matrix' The following object(s) are masked from 'package:base': det Anexando pacote: 'lme4' The following object(s) are masked from 'package:stats': AIC, BIC
sigma <- 1 # desvio-padrão do efeito de lote i <- 10 # níveis de lote j <- 4 # número de repetições por lote m <- 50 # número de sementes por ensaio b <- rnorm(10, 0, sigma) # efeitos dos lotes z <- gl(i, j) # vetor de níveis de lote y <- rbinom(i*j, size=m, prob=1/(1+exp(-rep(b, e=j)))) # germinação observada plot(y, col=z) mm0 <- glmer(cbind(y, 50-y)~(1|z), family=binomial) summary(mm0) Generalized linear mixed model fit by the Laplace approximation Formula: cbind(y, 50 - y) ~ (1 | z) AIC BIC logLik deviance 77.98 81.36 -36.99 73.98 Random effects: Groups Name Variance Std.Dev. z (Intercept) 0.66017 0.81251 Number of obs: 40, groups: z, 10
Fixed effects: Estimate Std. Error z value Pr(>|z|) (Intercept) -0.1926 0.2615 -0.736 0.462
O que esse "vanilla" faz e como eu posso usar o R assim por meio do Emacs, o meu editor de código? Grato. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================