search for: pung

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

...uot;6","7","8","9")))) # grid.lines(x = visit,y = var,draw = TRUE) ?? I am thinking I may need to use a function that joins coordinates (for example join (1,826) with (2,900)), but am hoping there may be a better way. Thanks for any help. Murray -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

...ot;,cex = 1.0), scales = list(relation = "free",x = list(tick.number = 1,at = c(1,3)),y = list(limits = ylimit)), auto.key = list(space = "right",cex = 1.1), par.settings = list(superpose.symbol = list(pch = symbols,cex = 1.1))) -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

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

...ould 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 = "n", bty = "n" ) Thanks Murray -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

...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 second function? (I have tried defining it as a parameter when calling the 2nd function from within the 1st). Cheers! -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

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

..., 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0)), .Names = c("patient", "score", "group"), class = "data.frame", row.names = 1:20) attach(test) colr <- c("gray","pink") barplot(score,beside = T,space = .4,col = colr[test$group]) -- 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]]

...es? 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 Maintenance"), lty = 2:3) Thanks Murray -- Murray Pung

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

...;)) densityplot(~try) normtest <- shapiro.test(try) normtest pvalue <- round(normtest$p.value,5) normtext <- paste(normtest$method,"p-value =",pvalue) normtext xcoord <- max(try)*0.6 text(xcoord,0.1,normtext) # alternative plot(density(try)) text(0,0.1,normtext) -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[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

...) ~ group, data = surv, rho = 0) N Observed Expected (O-E)^2/E (O-E)^2/V group=1 20 16 11.0 2.23 4.64 group=2 20 12 17.0 1.45 4.64 Chisq= 4.6 on 1 degrees of freedom, p= 0.0312 > survtest$chisq [1] 4.641028 > survtest$df NULL -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

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.

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

..."Mean of Variable", xtick = TRUE, xlab = "Day", x.leg = 1, y.leg = 25, cex.leg = 0.9, lty = c(1,5) ) detach(try) rm(try) -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

Dear R-Users, given you have been teaching R to students (grad level, mainly social science background, no previous programming experience, 80% know SPSS), what are your experiences concerning the style of teaching? Do you prefer to stand in front of the class like in "normal" lectures and you show them slides? Or do you you explain some concept (for example things like

I have specified error = "greenwood" in the survfit function: summary(survfit(Surv(time, status) ~ 1,error = "greenwood", data=aml)) The output does not appear to give variance, but standard error. Preferably I would like variance to be output, but is it okay to convert std errror to variance, and how can I do this correctly? Thanks Murray [[alternative HTML version