similar to: error bars in lattice barchart

Hello, I am attempting to write a script that adds error bars to a barchart. I basing my attempt heavily on the following thread: http://tolstoy.newcastle.edu.au/R/e2/help/06/10/2791.html I can't seem to get around the problem that was discussed in the thread. The following example should illustrate my problem. Sorry about the messy example but I am 1) trying to make it as close as possible

Hi all, I'm trying to plot error bars in a lattice plot generated with xyplot. Deepayan Sarkar has provided a very useful solution for simple circumstances (https://stat.ethz.ch/pipermail/r-help/2005-October/081571.html), yet I am having trouble getting it to work when the "groups" setting is enabled in xyplot (i.e. multiple lines). To illustrate this, consider the singer data

Hi, all, I'm trying to modify the code to a log scale for y-axis from the post http://tolstoy.newcastle.edu.au/R/help/06/06/28612.html However, the error bar did not change accordingly. The following is the code I used based on the singer.ucl data. Thanks in advance! prepanel.ci <- function(x, y, ly, uy, subscripts, ...) { y <- as.numeric(y) ly <-

All, I'm trying to adapt some code provided by Deepayan Sarkar from a previous thread (https://stat.ethz.ch/pipermail/r-help/2005-October/081571.html) on this topic. ## This code produces a graph w/o error bars: xyplot(Y ~ Hr, data, groups=DRUG, panel=function(x,y,...){ panel.xyplot(x,y,..., type=c("g", "l") ) panel.points(x,y,..., pch=16,

All, I have data across 5 time points that I am graphing via xyplot, along with error bars. For one of the variables I have missing data for two of the time points. The code below is okay but I can't seem to get the lines to connect across the missing time points. Does anyone now how to rectify this? Cheers, David Afshartous library(lattice) ## the data junk = data.frame( Visit =

Dear R users, I am trying to use xyplot to draw group mean and CI. The following is the sample code. But I want: 1. Use different colors and symbols to draw individual points, CI and the lines connect group means from different time points; 2. Add jitters to x axis to allow CIs not be overlapped Could anyone modify the attached code to achieve this? Thanks library(lattice)

Dear R users, About a year ago Deepayan offered a suggestion to incorporate error bars into a dotplot using the singer data as an example <<http://finzi.psych.upenn.edu/R/Rhelp02a/archive/63875.html>>. When I try to utilize this code with a grouping variable, I get an error stating that the subscripts argument is missing. I have tried to insert them in various ways, but cannot

Hi, [I'mm cc-ing r-help on this, since some others might be interested (Tony Rossini is, at least)] Yes, I saw your mail and was wondering about what would be the best approach to do it. It's definitely not trivial, mainly because Trellis is not designed to handle factors on the x-axis (chiefly for convenience, as factors typically have long names which are difficult to write on the

I would like to place the value for each bar in barchart (lattice) at the top of each bar. Something like the following code produces. library(lattice) mypanelfunc <- function(x, y, ...) { panel.barchart(x, y, ...) panel.text(x, y, labels=as.character(round(x,2)), ...) } myprepanelfunc <- function(x, y, ...) list(xlim=c(0, max(x)+.1)) mydata <- expand.grid(a=factor(1:5),

Hi. I've got a lattice plot with multiple panels and two groups superimposed on each panel. Each panel has an independently scaled y-axis (scales = list(relation = "free")). I've successfully put up 95%CI error bars using panel.arrows (and some help from the mailing list). My question is whether I can "unscale" the arrowheads so that they appear to have the same

Dear R-users, i want to plot gene regulation data in a lattice barchart. To illustrate the problem i encounter, the following code uses the "barley"dataset: library(lattice) barley[["SD"]] <- 5 PLOT <- barchart(data=barley, yield~variety|site, groups=year,origin=0, as.table=TRUE, scales=list(x=list(relation="same", rot=30), y=list(alternating=3,tck=-1)),

Hi all, I'm hoping someone could tell me what best practices are as far as keeping programs organized in R. In most languages, I like to keep things organized by writing small functions. So, suppose I want to write a function that would require helper functions or would just be too big to write in one piece. Below are three ways to do this: ################### Style 1 (C-style)

Hi, I'm trying to calculate a ratematrix for a RNA aligment (U instead of T) in order to use it as a ratematrix in Profidst (a phylogenetic program which takes into account both the primary sequence and the secondary structure of the RNA, in my case rRNA). The sequence-structure aligment has been made in 4SALE (a java app) and saved as one-letter encoded (using a 12 letters alphabet,

Hi, The problem is that I would like to do an interaction plot with intervals based on Tukey's honestly significant difference (HSD) procedure, but I do not know how to do it in R. I have 3 factors "A", "B" and "C" and a response variable "response". I would like to study a model where there are main effects and second order interaction effects.

Hi all, I already posted this to seqinr-forum@lists.r-forge.r-project.org but since (so far) no reply was given, I decided to try my luck in this list too. Dear list. I wish (for my thesis work), to import tRNA data into R and have it aligned. My questions are: 1) What resources can I use for the data. 2) What commands might help me with the import/alignment. So far, I found two nice

Hi, I have a CSV file with various biological reactions. Subscripts, superscripts, and italics are encoded in carats, and I was wondering if R can actually recognize those and print actual superscripts, etc. Here's an example: <i>S</i>-adenosyl-L-methionine + rRNA = <i>S</i>-adenosyl-L-homocysteine + rRNA containing

Hi, The functions prepanel.default.bwplot() and lpretty() are not running in Deepayan's barley example concerning vertical bars with the lattice function barchart(). Why, is there a restricted use for the package lattice-Internal? Urs Simmen mailto:usimmen at dtc.ch

Hi everybody, I'm having problems with the barchart command in the lattice package. I'm creating barcharts from matrices with with anything from 20 to 71 columns. When I leave the column names alone, that is they are set in the read.table command (and inherited by subsequent commands) the correct number barcharts is created by the "print(barchart(...))" command.

Dear list member, I struggle with the problem, why the probabilities of choosing one of three mutually exclusive alternatives don?t sum up to 1! Let?s assume we have three alternatives X, Y, and Z. Let?s further assume we know their respective utilities: uX, uY, uZ. I?m interested in calculating the probability of choosing X, Y, and Z. Since I assume that the alternatives are mutually

Dear all, How can I make an index plot with lattice, that is plotting a vector simply against its particular index in the vector, i.e. something similar to y <- rnorm(10) plot(y) I don't want to specify the x's manually, as this could become cumbersome when having multiple panels. I tried something like library(lattice) mp <- function(x, y, ...) { x <- 1:length(y)