similar to: global object in user defined function

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

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 =

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

I would like to add an extension to the current name of a variable to create a new variable that is its sqrt transform. Each piece of the equation below works independently, but the left side definition fails on run. I also tried creating the variable name first, but ended up with an object that toString() did not fix. Better ideas? Example 1

Hi everyone, I'm a new R user (if this is a really basic question, please do excuse me...) and I'm having some questions regarding a deciles problem. I have a variable which I need to categorize according to its deciles (X). However, this categorization should be made into another variable (call it NewVar). Ex. for the quartiles case (just for the sake of exposition, since I need

Hi, there I'm currently trying to figure out how to keep my "factor" levels for a variable when moving it from one data frame or matrix to another. Example below: vec1<-(rep("10",5)) vec2<-(rep("30",5)) vec3<-(rep("80",5)) vecs<-c(vec1, vec2, vec3) resp<-rnorm(2,15) dat<-as.data.frame(cbind(resp, vecs))

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

Hello, If I have two variables that are factors or characters and I want to create a new variable that is the combination of both what function can I use to accomplish this? Ex. Var1 Var2 SA100055113 19851113 And I want NewVar SA10005511319851113 Thanks in advance. Cameron Guenther, Ph.D. Associate Research Scientist FWC/FWRI, Marine Fisheries Research 100 8th Avenue S.E. St.