similar to: differentiate groups on barplot

I would like to join repeated measures for patients across two visits using a line. The program below uses symbols to represent each patient. Basically, I would like to join each pair of symbols. library(lattice) patient <- c(1,2,3,4,5,6,7,8,9,1,2,3,4,5,6,7,8,9) var <- c(826,119,168,90,572,323,122,10,42,900,250,180,120,650,400,130,12,33) visit <- c(1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2)

I have an outlier that I would still like to display, but would prefer to shorten the axis. For example, display 0% - 40%, and 90% - 100%. Is this possible? I am using an xyplot. Thanks Murray -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

I would like to set the y axis limit of an xyplot using the object 'ylimit', but receive this error: [1] 990 Error in extend.limits(limitlist[[i]], axs = axs) : improper length of lim I get the same error if I use ylim. library(lattice) trellis.device(col = FALSE, theme = lattice.getOption("col.whitebg")) name <- "Variable name" symbols <-

Is there any way of adding text to a density plot? I have had a go using the text() function but I think the error is because this function doesn't work with densityplot(). Alternatively, I understand I can achieve pretty much the same result if I plot a density kernel estimate using plot() (which allows text()), but I do prefer densityplot(). Also, is it possible to specify the dimensions

I have checked out the help files, but cannot find details on how to use maths characters in ylab. Instead of m^2, I would like the 2 in superscript, if possible. I would also like to place more padding on the label so that the label is not obscured by the horizontal numbers. y <- 1:10 x <- rnorm(10,50000,2000) plot(x ~ y, ylab = 'Y Label (m^2)', las = 1, type =

I need an object created in a user defined function to be accessible to another user defined function. I am fairly certain the object is correctly created as it prints when I use the return() function, however the 2nd function seems unable to use it. I am guessing objects created in functions are not globally accessible, is this correct? What is an appropriate way to pass the object to the

Is it possible to add standard error bars to the means on interaction plots? Thanks Murray -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

Is there an implementation of the Cox-Snell residuals / Nelson-Aalen plot for goodness of fit? Or otherwise is there an appropriate Goodness of Fit diagnostic? Thanks Murray -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

How can I get the Log - Rank p value to be output? The chi square value can be output, so I was thinking if I can also have the degrees of freedom output I could generate the p value, but can't see how to find df either. > (survtest <- survdiff(Surv(time, cens) ~ group, data = surv,rho=0)) Call: survdiff(formula = Surv(time, cens) ~ group, data = surv, rho = 0) N Observed

Are there ways to load a csv file without row.names by read.csv? Thanks. [[alternative HTML version deleted]]

Is it possible to format the x axis, so that days are spaced proportionately? i.e. a larger space between 50 & 100, than 0 and 5? Have not had any luck with axis.POSIXct(). Many thanks Murray try <- structure(list(visit = structure(as.integer(c(1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5)), .Label = c("Screening", "0",

R help - I am attempting to write a script that has multiple subjects in 1 data file. Each subject has multiple rows with columns as variables. Here is my code, I am having problem executing it on each unique subject id (dat$Subject). getwd() setwd("/Users/edwardpatzelt/Desktop/Neuroimaging/MERIT/SRRT/merge") dat <- read.table("test2.txt", header = TRUE, na.strings

I am using the R code listed below to create 4 smoothed lines on a xyplot. I'm having trouble fine tuning it. First I think I may need a black and white plot so how do I get it to plot the lines with different characters, preferable the same characters used in the key (plus, X circle and triangle). I might also be interest in a version that draws four solid lines of different colors but

Dear all, I have a 2 by 2 matrix and I would like to do a barplot with it. (so 2 bars with each having 2 stacks.). I would like to have one colors per stack, so 4 different colors total. The problem is that R is only given me 2 colors (the same two for the bottom stack and the same two for the top stack). Any idea how I can do to have 4 colors? (without using ggplot2 preferably) Here is my

I would like to add a footnote to this graph but do not see a "footnote" command in the package:lattice documentation. I would like to note the "span=.8" as the footnote. postscript(file= ?C:/Documents and Settings/dsonneborn/My Documents/Slovak/output/pcb_tables/smooth_PCB_lines_four.ps?, bg=?transparent?, onefile=FALSE, pointsize=20,paper=?letter?, horizontal=TRUE,

I am using plot(survfit(Surv(time,status) ~...) and would like to add error bars rather than the confidence intervals. Am I able to do this at specified times? e.g. when time = 20 & 40. leukemia.surv <- survfit(Surv(time, status) ~ x, data = aml) plot(leukemia.surv, lty = 2:3,xlim = c(0,50)) #can i add error bars at times 20 & 40? legend(100, .9, c("Maintenance", "No

Dear R wizards, I am seeking advice on graphics in R. Specifically, how to manipulate the size and save a plot I have produced using the LDheatmap library. I confess I am relatively new to graphics in R, but I would greatly appreciate any suggestions you may have. LDheatmap produces a coloured triangular matrix of pairwise associations between 600 genetic markers in my dataset.

Hello there! I have a little trouble with using dovecot as lmtp and maillists. I use dovecot 2.0.13 with exim 4.72 on Debian 6.0. (All stuff at one host). info of maillists is store in mysql table, for example: +---------+-------------------------------------------------------------------------------------------------+ | name | recipients

Hi everybody, i've been using rsync for quite some time to sync files between my laptop and my desktop for example i use the following to backup the music library (from Mac to Linux) > rsync -rvzu --exclude=.DS_Store --delete --exclude 'Podcasts/*' > $LOCAL_MUSI $REMOTE_MUSIC >> $LOGFILE An odd thing i always noticed is that some files were always copied, even

Dear fellow R users, I would like to draw a "sunflowerplot" because I have data (decade by month) that plots multiple times on the same x-y co-ordinates. Further I would like to colour each of the points/sunflower leaves on the plot according to the group they belong to (i.e. which type of event each represents within that decade and month). I thought that this would be relatively