similar to: maths characters in labels & ylab padding

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

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

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

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

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

To differentiate between groups on the barplot, I guessed that col = colr[test$group] would have worked. How can I do this? Many Thanks Murray test <- structure(list(patient = 1:20, score = c(100, 95, 80, 75, 64, 43, 42, 40, 37, 35, 30, 29, 27, 26, 23, 22, 19, 18, 17, 16), group = c(1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0)), .Names = c("patient",

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

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.

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

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

I'd like to paste a superstring with a number in an object. Thanks for any help. Murray mycor <- cor(1:10,1:10) plot(1:10,1:10) text(8,2,paste(expression(R^2)," = ",mycor)) [[alternative HTML version deleted]]

Dear List I am new to R and to the list and will try best as I can be clear and concise. My apologies if anything I write contravenes the posting code on this list. I would also like to say I have run through most of the material on the R website before writing this email however, I am stuck. Here is what I want to do and what I have done 1. Read a comma seperated text file into R I have used

Full_Name: Murray Pung Version: 1.9.0 OS: OSX Mac Submission from: (NULL) (134.148.20.33) Date/Time: 2004-06-05 12:32:30 +1000 OS Version: 10.3.4 (Build 7H63) Report Version: 2 Command: R.bin Path: /Library/Frameworks/R.framework/Resources/bin/R.bin Version: 1.9.0 (R 1.9.0) PID: 358 Thread: 0 Exception: EXC_BAD_ACCESS (0x0001) Codes: KERN_PROTECTION_FAILURE (0x0002) at

hi all, in R, what's the best way to have a rotated ylab in a graph plotted with either xyplot or xYplot? I tried this but it didn't work. xyplot(y~x,data.frame(x=1:10,y=runif(10)),ylab=list(srt=90,crt=90,rot=90,label="my label")) more generally, how do you output a text at an angle in a lattice graph? what would be a good reference for R lattice graphics? I need more

Hi! How thy font size of xlab & ylab can be set? Eryk

Hi all, Has anyone tried to plot math symbols with vis.gam() - a function that uses (I think) plot.gam() internally? This is from the package mgcv - I am new to gam(). For example, with the normal plot(), the expression() works fine for Math symbols: plot(xx, ylab=expression(paste(Delta,"H",sep=""))) # works fine But, with vis.gam(), the same

Hi All, I'm trying to write "ug/m3" as y-label, with greek letter "mu" replacing "u" AND "3" going as a power. These commands works in general: plot.new() text(0.5, 0.5, expression(symbol("m"))) But, I'm sure about how to do it using barchart() from "Lattice". Can anyone help please? Thanks, Peng Cai [[alternative HTML