dotplot question
I got it. I have to detach gdata, gplot, gtools and gmodels. Thank you very much!
On 3/27/2011 11:59 AM, Sebasti?n Daza wrote:
This is what I get: R version 2.12.2 (2011-02-25) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] grid splines stats graphics grDevices utils datasets [8] methods base other attached packages: [1] gregmisc_2.1.1 gplots_2.8.0 caTools_1.11 bitops_1.0-4.1 [5] gtools_2.6.2 gmodels_2.15.1 gdata_2.8.1 foreign_0.8-42 [9] Hmisc_3.8-3 survival_2.36-5 lme4_0.999375-39 Matrix_0.999375-48 [13] lattice_0.19-17 rj_0.5.2-1 loaded via a namespace (and not attached): [1] cluster_1.13.3 MASS_7.3-11 nlme_3.1-98 rJava_0.8-8 stats4_2.12.2 [6] tools_2.12.2 Thank you! On 3/27/2011 11:39 AM, Douglas Bates wrote:
2011/3/26 Sebasti?n Daza<sebastian.daza at gmail.com>:
Dear Douglas, I ran again the Dyestuff example, and I got a messy dotplot (see attachment). fm1<- lmer( Yield ~ 1 + (1| Batch ), Dyestuff ) dotplot(ranef(fm1, postVar = TRUE), strip =FALSE, scales=list(y=list(col="white")))
ranef(fm1)
$Batch (Intercept) A -17.6080025 B 0.3912889 C 28.5640930 D -23.0860478 E 56.7368971 F -44.9982287
Could you send us the output of sessionInfo() from your R session?
Thank you. Sebastian. On 3/21/2011 11:49 AM, Douglas Bates wrote:
2011/3/14 Sebasti?n Daza<sebastian.daza at gmail.com>:
Dear list,
I have a couple of questions related to dotplots and ranef:
1) Does anyone know how to order the distribution of a residuals in a dotplot:
dotplot(dotplot(ranef(model, postVar = TRUE)...
sometimes residuals are ordered but other times they are not.
Can you provide a reproducible example?
2) Is it possible to organize these plots in a matrix (as par function does)? I have more than a random term... so would be useful to order theses plots in a matrix...
We have good news and bad news. The good news is that you get more flexibility with lattice graphics than with traditional R graphics (which is what the par function controls). The bad news is that it requires a bit more effort to control the lattice graphics display in that you need to read the manual page for print.trellis carefully. An example of the dotplot for more than one grouping factor is provided on slide 31 of http://lme4.R-forge.R-project.org/slides/2011-03-16-Amsterdam/1Simple.pdf The code that generates that figure is qrr2<- dotplot(ranef(fm2, postVar = TRUE), strip = FALSE) print(qrr2[[1]], pos = c(0,0,1,0.75), more = TRUE) print(qrr2[[2]], pos = c(0,0.65,1,1))
3) Is it possible to change the strip name of these plots?
Yes. See the description of the strip.default function in the lattice package. In particular, the strip can be suppressed by setting strip=FALSE, as above.
Any clue about these questions will be appreciate. Thank you in advance! -- Sebasti?n Daza sebastian.daza at gmail.com
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
-- Sebasti?n Daza sebastian.daza at gmail.com
Sebasti?n Daza sebastian.daza at gmail.com