similar to: retrieving stats from bwplot

I?m entirely stumped on this particular issue and really hoping someone has some advice for me. I am running a covariant model in lm I would like to give the standard errors or the confidence intervals for the fitted lines. I?ve been using the dataset OrangeSprays where I want lines for each level of treatment over the covariant ?colpos?. I?ve been able to calculate intercepts and slopes for

Hi, I wonder if it is possible with the graphics package to overlay one or several plots (barplots, for example) over a map. Data for the map is in a data frame with the latitude and longitude coordinates, and then: > plot(map$long, map$lat, type ="l") produces the map. I want to put each barplot in specific locations on the map, namely at the center of "statistical

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 treatmentF treatmentG treatmentH

Hi, I cannot seem to get lme4 to work. I have installed the lme4 and Matrix package with apt-get. and both can be found in /usr/lib/R/site-library. When I tried an example for lmer, R could not find the function lmer(), > library(lme4) Attaching package: 'lme4' The following object(s) are masked from package:nlme : getCovariateFormula getResponseFormula

Hello all, I wrote a small function to add labels for outliers in a boxplot. This function will only work on a simple boxplot/formula command (e.g: something like boxplot(y~x)). Code + example follows in this e-mail. I'd be happy for any suggestions on how to improve this code, for example: - Handle boxplot.matrix (which shouldn't be too hard to do) - Handle cases of complex

Dear list, I would like to plot arrows with different colors according to arrow length, and also (if possible) with head size proportional to arrow length. The idea is to make a quiver-like plot of matlab with wind speed data. So far, I´ve been able to use different colors, but I need to find a more efficient way to recode arrow length intervals into colors. On the contrary, I can't define

Dear R users I use bwplot to plot some figures. There are two troubles: 1. How to change the dot of the mean to a line, like the style in boxplot 2. How to hide some tick marks. For example, I have seven tick: A, B, C, D, E, F, G, but I want show four marks: A, C, E, G on the x-axis. Thanks! -- Kind Regards, Guanghong [[alternative HTML version deleted]]

Hi, I'm sure this one is very easy.... I am trying to write a function where one of its arguments has two posible (strings) values, defaulting to one of them if none is specified. My problem is that when evaluating the function the following warning is produced: "the condition has length > 1 and only the first element will be used" I've read the help page of if() but I

Junto con saludar, Quer?a saber si los que han usado boxplot.stats, saben como ubicar los resultados que arroja como "out" (outliers) directamente en la base de datos, ya que, muestra los resultados en valor de cada dato at?pico. Por ejemplo, estoy usando una base de datos con 300 datos y 10 columnas, en este caso estoy buscando los datos at?picos de la columna del precio de x

Hi, is there an equivalent to the 'range' option of the boxplot function to be found in the bwplot function of the trellis package? regards Vincent [[alternative HTML version deleted]]

Searched for this and found some help, but I still can't figure it out. I have trying to enlarge the x and y labels on my box plot. I understand that you can do this using "cex.lab", but it does not seem to be working for me. I must be adding it in the wrong spot. Any help would be greatly appreciated. Here is my code: bwplot(hr~Herd, data=telemetry, notch=T, ylab="Home

Hi, I'm having difficulty to replicate in R a nested (hierarchical) anova example found in p. 308 of Zar, J.H. 1996. Bostatistical Analysis. Prentice Hall. 3rd ed. The example (15.1) is as follows: The variable is blood cholesterol concentration in women (in mg/100 ml of plasma). This variable was measured after the administration of one of three different drugs, each drug having been

*This message was transferred with a trial version of CommuniGate(tm) Pro* A simple question, perhaps to simple, but I am new to R: For a data frame with k variables and n1+n2 observations (two groups) how can i substract the group means for every variable, so both group variables become centered to mean zero? Thanks in advance for your help H?ctor

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

Hello, In doing a series of non-linear estimations of a function which is a sum of a varying number of sinusoids, I would like to "autogenerate" the arguments needed by nls() depending on that number. For example, when there are two sinusoids: > nls( y ~ mu + A1 * cos(2*pi*f1*x - P1) + A2 * cos(2*pi*f2*x - P2), data = some.xy.data, start = list( mu=some.value0,

hello everybody, I want to print two plots in one png file, I tried several options 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

Hi everybody, Excuse me for this silly question, but after searching the help archives I'm still unable to find my way to read the source code of the "biplot" function in R. I have installed the mvbutils package, and tried: > fixr(biplot) which only gives me: function (x, ...) UseMethod("biplot") with no further details.... I'd like to read the code to see if

Hi, how can I replace NA value with 0: 1991 217 119 103 109 137 202 283 240 146 NA 1992 270 174 149 144 166 239 278 237 275 NA 1993 146 111 104 89 98 131 153 148 175 NA 1994 177 123 146 124 121 200 266 191 240 106 1995 145 98 95 89 95 130 183 161 164 129 1996 145 98 89 90 93 138 158 131 161 161 1991 217 119 103 109 137 202 283 240 146 0 1992 270 174 149 144 166 239 278 237

Hello, I am trying to use "bwplot" to display whisker boxes according to some conditioning factors ("age" has two values 1/2). I get the following messages: > library(trellis) > bwplot(dvk95~age|site*season,panel=function(x,y){panel.bwplot(x,y)}) Error in tapply(1:0, structure(list(INDICES = numeric(0)), .Names = "INDICES"), : arguments must have

I want to turn off selected axes in bwplot. I would like to only have the bottom axis drawn, with the others off. I have a series of bwplots that I want on one device, like this: p1<-bwplot(x1,box.ratio=.1) p2<-bwplot(x2,box.ratio=.1) ... print(p1,position=c(0,.8,1,1),more=T) print(p2,position=c(0,.6,1,.8),more=T) ... I know about the panel functions panel.bwplot() and panel.axis(), but