search for: datapharm

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

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

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

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

Are there ways to load a csv file without row.names by read.csv? Thanks. [[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.

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

...able", 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]]

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

...y) 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]]

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