similar to: use xyplot to plot mean and CI by groups

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

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 =

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've read the emails of Dan, Deepayan and Sundar about adding error bars to the lattice plots ( https://stat.ethz.ch/pipermail/r-help/2006-October/114883.html), but I still have the problem when I want to adding error bars to barchart. I tried both the solution of Deepayan and Sundar but without luck. Here is my code (I changed prepanel.ci and panel.ci a little to plot bars

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,

Dear R help, I tried to add linear fit lines in intxplot by adding "index.cond = function(x,y) coef(lm(y ~ x))[1]" inside the intxplot().? It didn't help.? I would appreciate any help. My dataset and codes are pasted below. Thanks, A.K. datGreen<- read.table(text=" Time,??? DarkSt,??? LightSt 0,??? 0,??? 1 1,??? 0.03,??? 0.82 2,??? 0.06,??? 0.91 3,??? 0.12,??? 0.64 4,???

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 <-

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 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

Hi, How to unlist list column value and add column into data frame. Data frame ID ContractDe PassengersDe TrainnerDe 1 list(ConID=c("Zx","78yu"),ConRes = c("98","Tut")) list(PassID

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, 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.

In the example: x = rep(c(0,0,1,1),4) y = rep(c(0,1,0,1),4) z = c(1,0,1,0,0,0,1,1,0,1,0,0,1,1,1,1) f = as.factor(c(rep("a",4),rep("b",4),rep("c",4),rep("d",4))) levelplot(z~x+y|f,data.frame(x=x,y=y,z=z,f=f)) I noted that the last ("d") plot remains empty. I guess the reason for this is that the values are constant (1), but I consider it more

One of many scripts to produce 4 lattice plots on one page keeps throwing an error. I've tried manipulating the file to eliminate the error, but have not been able to do so. The error is: > source('bicarb.R') Error in source("bicarb.R") : bicarb.R:15:1: unexpected symbol 14: 15: hco33 ^ The 'h' is in column 0 so the caret would be column -1, but it's

hi... i want to use maildir format so like i see un some manuals i create this directory for my virtual users: # ls -la /export/home/vmail/prueba.uy/t* /export/home/vmail/prueba.uy/toto1: total 4 drwxr-xr-x 2 vmail vmail 512 Apr 10 15:09 . drwxr-xr-x 6 vmail vmail 512 Apr 9 12:11 .. /export/home/vmail/prueba.uy/toto3: total 10 drwxr-xr-x 5 vmail vmail 512 Mar

Hi there, I'm a newy at samba 4 and I'm trying to joing a samba 4 alpha 17 box to our domain as a DC so I can "drain" the domain info and use the linux box to test samba without disturbing the domain itself. The thing is that I get an error when trying to join the server : [root at vpdc samba]# bin/samba-tool domain join montecarlotv.com.uy DC -Uadministrador

First a bit of background: After doing a search for a flexible polar plot function and coming up empty, I have begun writing one myself. Since I am new to writing mid-level graphics routines, this has required some learning about lattice, grid and related things. I am to the point where I have a workable proof of concept, but still need to make some improvements. My goal is to have

Dear friends - please give me a hand. I have a dataset of 40 patients in two groups observed on three occasions. I only want to plot a line for each patient in the two groups. I use the ltext function to put the patient number but fail to make lattice understand the numbers as unique since apparently it starts all over with 1:20 for each panel instead of respecting my desires to have 1:20

Hello, R users, I have two factors (treat, section) anova design experiment where there are 3 replicates. The objective of the experiment is to test if there is significant difference of yield between top (section 9 to 11) and bottom (section 9 to 11) of the fruit tree under treatment. I found that there are interaction between two factors. I wonder if I can contrast means from levels of

Dear Thomas Lumley, and R-help list members, I have read your article "Network meta-analysis for indirect treatment comparisons" (Statist Med, 2002) with great interest. I found it very helpful that you included the R code to replicate your analysis; however, I have had a problem replicating your example and wondered if you are able to give me a hint. When I use the code from the