search for: orchardsprays

...too long. The comparable line with shorter variable names is Browse[1]> x [[1]] [1] "factor(-(r <= 2) + (r >= 6), labels = c(\"1:2\", \"3:5\", \"6:8\"))" I am using rw2001 ## Code to duplicate the problem: options()$width library(lattice) data(OrchardSprays) OrchardSprays$r <- OrchardSprays$rowpos xyplot(decrease ~ treatment | factor(-(rowpos<=2)+(rowpos>=6), labels=c("1:2","3:5","6:8")), data=OrchardSprays) xyplot(decrease ~ treatment | factor(-(r<=2)+(r>=6), labels=c("1:2&quot...

I feel bad even asking, but: Rgames> data(OrchardSprays) Rgames> model<-lm(decrease~.,data=OrchardSprays) Rgames> model Call: lm(formula = decrease ~ ., data = OrchardSprays) Coefficients: (Intercept) rowpos colpos treatmentB treatmentC 22.705 -2.784 -1.234 3.000 20.625 treatmentD treatmentE...

Hi All, I'm trying to draw boxplots. I'm having a hard time to get "ticks labels" on multiple panels using 'alternating' option. # R Code: # May not be the best example, please just look into 'scales' option library(lattice) data(OrchardSprays) dta <- subset(OrchardSprays, OrchardSprays$rowpos %in% c(1,2,3)) # Original # This works fine, as you can see Y-AXIS labels (0,50,100) on both Panel-1 and Panel-3, and # X-AXIS lables (A,B,C,...) on both Panel-1 and Panel-2 bwplot(decrease ~ treatment | factor(rowpos), data=dta, pch="|&...

Full_Name: James Curran Version: 2.7.1 OS: Windows XP/Vista Submission from: (NULL) (130.216.51.81) stripchart ignores the main title option. For example, in the documentation stripchart(decrease ~ treatment, main = "stripchart(OrchardSprays)", vertical = TRUE, log = "y", data = OrchardSprays) does not yield a plot with "stripchart(OrchardSprays)" at the top of the plot. Looking at the source code of stripchart.default there is: .... title(xlab = xlab, ylab = ylab) .... I would have thought the ma...

Hi All, Can anyone please help me with getting a "single title" and "legend" for both the plots in the following R code. I'll eventually be using .wmf file. # R code: library(lattice) p1 <- xyplot(decrease ~ treatment, OrchardSprays, groups = rowpos, type = "a", main="Same title / legend", auto.key = list(space = "right", points = FALSE, lines = TRUE)) p2 <- xyplot(decrease ~ treatment, OrchardSprays, groups = rowpos, type = "a", main="Same title / legend&quo...

...point). However, I cannot add line legend and point legend together. The code looks like the following. (1) First, I tried the following code, however, it plot a line across the point legend("*"). boxplot(count ~ spray, data = InsectSprays) rb <- boxplot(decrease ~ treatment, data = OrchardSprays, col="lightgrey") title("Comparing boxplot()s and non-robust mean +/- SD") mn.t <- tapply(OrchardSprays$decrease, OrchardSprays$treatment, mean) xi <- seq(rb$n) points(xi, mn.t, col = "blue", pch = "*") h1 <- abline(h=110, col="red") h2 &...

...ize in lattice when changing resolution? ---- Example: ..... ###This give a totally squished graph, where the actual plotting area is minimal CairoPNG("example.png",width = 480, height = 480, dpi=600, pointsize = 12, bg = "white") bwplot(decrease ~ treatment, OrchardSprays, groups = rowpos, panel = "panel.superpose", panel.groups = "panel.linejoin", xlab = "treatment", key = list(lines = Rows(trellis.par.get("superpose.line"), c(1:7, 1)), text = list(lab = as.chara...

Hi, I am trying to plot a graph using stripchart() - my x-axis is categorical and my y-axis is numerical. I would like to reverse the order of my y-axis (i.e. with 0 at the top left and increasing numbers moving downward towards the bottom left corner). My data are subsetted into 3 categories that I've plotted on 1 graph using stripchart(add=TRUE) option. Thanks for any assistance, Regards,

...e right. Somehow the internal calculation which computes horizontal space for the legend is not flexible enough. ## Call a new graphics window {with default par()s !}: get(getOption("device"))() par(mfrow = c(2,2)) ## part of example(interaction.plot) _improved_ using with() : data(OrchardSprays) with(OrchardSprays, { interaction.plot(treatment, rowpos, decrease) interaction.plot(rowpos, treatment, decrease) }) data(esoph) with(esoph, { interaction.plot(agegp, alcgp, ncases/ncontrols) interaction.plot(agegp, tobgp, ncases/ncontrols, trace.label="tobacco",...

I am a new and unexperienced user of R and got so far as to know how to produce boxplots. I have no experience of messing with function code, so presently I do not know how to create a boxplot with group means instead of group medians. If somebody could help me either replace the median with the mean or superimpose the mean onto the existing boxplot, it would be appreciated.

...ptions but i didn't succeed the first plot (bwplot) print to the defined position, but the second (ggplot) doesn't Any idea? Thanks a lot Christophe # Example: #------------------------------------- library(ggplot2) library(lattice) library(grid) one <- bwplot(decrease ~ treatment, OrchardSprays, groups = rowpos, panel = "panel.superpose", panel.groups = "panel.linejoin", xlab = "treatment", key = list(lines = Rows(trellis.par.get("superpose.line"), c(1:7, 1)), text = list(lab = as.chara...

...I?ve run the summary (model) function, but only get the differences in standard error between lines. Also by running confint(model), I have the same problem. How can I get standard error per or the CI-interval for each fitted line? I?ve attached the coding underneath. Thank you, Sigrid > data(OrchardSprays) > model<-lm(decrease~rowpos+colpos+treatment+treatment:colpos) > summary(model) -- View this message in context: http://r.789695.n4.nabble.com/se-s-and-CI-s-for-fitted-lines-in-multivariate-regression-analysis-tp4645703.html Sent from the R help mailing list archive at Nabble.com.

...el <- ourlier_df[ss,label_name] temp_x <- boxdata_outlier_df[i,"x"] temp_y <- boxdata_outlier_df[i,"y"] text(temp_x, temp_y, current_label,pos=4) } list(boxdata_outlier_df = boxdata_outlier_df, ourlier_df=ourlier_df) } # example: boxplot(decrease ~ treatment, data = OrchardSprays, log = "y", col = "bisque") boxplot.add.outlier.text(OrchardSprays, "treatment", "decrease", "colpos") ----------------Contact Details:------------------------------------------------------- Contact me: Tal.Galili@gmail.com | 972-52-7275845 Rea...

I have four boxplots, stacked on top of one another (vertically). There are five classes that the continuous variable is set against. The problem I have is the data for the second boxplot only contains values for four of the five classes. I would like to include the space in this graph so there is a continuity for all of the boxplots, in terms of number of classes. I would include a set of

...#39;panel' function and I'm aware of the option 'stats' which defaults to 'boxplot.stats' but I didn't understand it well and therefore I am unable to get what I need. Thanks in advance. Héctor Villalobos Windows XP, R 2.5.0 > bwplot(decrease ~ colpos | treatment, OrchardSprays, horizontal = FALSE,cex = 1, panel = function(...) {panel.bwplot(decrease ~ colpos | treatment, stats=boxplot.stats)}) Error in all(is.na(x) | is.na(y)) : argument "y" is missing, with no default In addition: Warning message: is.na() applied to non-(list or vector) in: is.na(x) [[alte...

...nlme : getCovariateFormula getResponseFormula groupedData Error in autoloader(name = "confint", package = "MASS") : autoloader did not find 'confint' in 'MASS' > (fm1 <- lmer(decrease ~ treatment + (1|rowpos) + (1|colpos), + OrchardSprays)) Error: couldn't find function "lmer" > Is this a bug with the lme4 package for Debian (r-cran-lme4)? Peter ######################### R : Copyright 2005, The R Foundation for Statistical Computing Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0 R is free so...

...m to support multiple print commands at once (please see below). Is there an alternative to get plots file which can be used in powerpoint/word? # R code: # this does not work; but same thing works with pdf() library(lattice) win.metafile("test.wmf") p1 <- xyplot(decrease ~ treatment, OrchardSprays, groups = rowpos, type = "a", auto.key = list(space = "right", points = FALSE, lines = TRUE)) p2 <- xyplot(sunspot.year ~ 1700:1988, xlab = "", type = "l", scales = list(x = list(alternating = 2)), main =...

Greetings I have a very simple question that I have not been able to solve by reading the manual. When I produce a stripchart with two straight columns of dots representing individual observations, one representing one group of subjects and the other representing another, the columns wind up at the far left and far right sides of the plot, and I'd like them to be closer together, about as far

Dear Rhelp, I would like to get alternate tick labels for the xyplot: library(lattice) library(grid) xy <-xyplot(decrease ~ treatment, OrchardSprays, ??? main= "Some plot", ??????? groups = rowpos, type = "a", ??? ??????????? page???? = function(n){ ???????????????????????? grid.text(LETTERS[j], ???????????????????????? y = 0.95, ???????????????????????? x = 0.15, ???????????????????????? default.units = "npc", ???...

...y part, but nevertheless, its a question. Is it possible to add titles to colorkey legends? As far as I can tell, there is a command to do it for normal "key" legends, but not for "colorkeys". eg it works for a normal key, created through auto.key xyplot(decrease ~ treatment, OrchardSprays, groups = rowpos, type = "a", auto.key = list(space = "right", points = FALSE, lines = TRUE,title="Key title")) but there is no comparable command for a colorkey x <- seq(pi/4, 5 * pi, length = 100) y <- seq(pi/4, 5 * pi, length = 100) r <- as....