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