similar to: Legend positioning in scaled survival plot

I posted several months about the problem with adding tick marks to curves using lines.survfit. This occurs when lines.survfit is used to add a curve to survival curves plotted with plot.survfit. The help for this function implies that mark.time=TRUE thus: plot(pfsfit,conf.int=FALSE,xscale=365.25,yscale=100,xlab="Years",ylab="% surviving",lty=2,mark=3)

I have attempted to follow posting guidelines but I have failed to find out what I am doing wrong here. I am trying to use lines.survfit to plot a second curve onto a survival curve produced by plot.survfit. In my case this is to be a progression free survival curve superimposed upon an overall survival curve, but I will illustrate my problem using the example given in the help for

I am sorry to trouble the list with this, but I'm hoping that someone will be able to help me track down this reference. I am looking for a good reference on repeated measures, and many of the relevant functions in R refer in their "help" to: Hand, D. J. and Taylor, C. C. (1987) _Multivariate Analysis of Variance and Repeated Measures._ Chapman and Hall. This sounds

Hi, there: It's my first time to post question in this forum, so thanks for your tolerance if my question is too naive. I am using a nonparametric smoothing procedure in sm package to generate smoothed survival curves for continuous covariate. I want to truncate the suvival curve and only display the part with covariate value between 0 and 7. The following is the code I wrote:

Deepayan: Sorry for the additional email, but I'm concerned that I wasn't being clear. In the end, I would like a graphic that portrays my y-values on a log scale in their native data range. That is to say, I would prefer to see the axis labels as 1000 rather than 3 (log10 of 1000). Thank you for your kind assistance, Hobie Perry St. Paul, MN Hobie Perry <hobie_perry@yahoo.com>

I am dying here with linphone (not sure if it is crap software or just me being an idiot) but out of the box debian installations of two linphones fail with a "Got SIP response 415 "Unsupported Media Type" back from 192.168.1.3" Can anybody recommend a particular SIP soft phone that broadly satisfies the following criteria? 1. Run on linux. 2. Simple to use and setup. 3. Is

Ok, I finally got the slideshow code to a state worth showing it off. The site is a very rough cut of a site I''m building for my wife''s photography, so ignore the unfinished design for now :) http://rachel.kathihill.com/ To see the ajax version, go to: http://rachel.kathihill.com/?ajax=1 To randomize the order the images show: http://rachel.kathihill.com/?random=1 To change

Dear R users, i try to calculate the probabilty to survive a given time by using the estimated survival curve by kaplan meier. What is the right way to do that? as far as is see i cannot use the predict-methods from the survival package? library(survival) set.seed(1) time <- cumsum(rexp(1000)/10) status <- rbinom(1000, 1, 0.5) ## kaplan meier estimates fit <- survfit(Surv(time,

Hi All, I want to take a simple random sample from a large dataset, gly, but I'm getting an error message. Any help? dim(gly) [1] 112371 37 > s1 <- sample(gly,100) Error in `[.data.frame`(x, .Internal(sample(length(x), size, replace, : cannot take a sample larger than the population when 'replace = FALSE' Thanks, Rachel [[alternative HTML version

Hi, I would like to view matrices I am working with in a clean, easy to read, separate window. A friend showed me how to do something like I want with edit(). I can view the matrix in the 'R Data Editor': For a sample matrix: > mat=matrix(1:15,ncol=3) > mat [,1] [,2] [,3] [1,] 1 6 11 [2,] 2 7 12 [3,] 3 8 13 [4,] 4 9 14 [5,] 5 10 15

Hello, I am trying to deal with a database table that contains a date in the form of the number of seconds since epoch. Data is inserted in this format by an existing script that I can''t change, so a schema change isn''t realistic. I would like to allow humans to view and edit this data in the form of an actual human-readable date. I figured a facade column would be the

Quick question please....A user logs into windowsXP and tries to create a folder/document and the ownership on the new file/folder defaults to nobody:nobody. I have the user set up in samba on the IRIX machine. All other users have no problem. Anyone have any suggestions? Thanks Rachel

I''m a rails newb and have been Googling about this, but I''m still stumped. Not showing everything here, but basically it should be a pretty common setup so I''m sure others know what I''m doing wrong. - A meter can belong to many meter_groups - A meter_group can have many meters. - A user can ''subscribe'' to viewing a meter_group (Subscription)

Hi All, I'm trying to automate a data summary using summary or describe from the HMisc package. I want to stratify my data set by patient_type. I was hoping to do something like: Describe(myDataFrame ~ patient_type) I can create data subsets and run the describe function one at a time, but there's got to be a better way. Any suggestions? Rachel [[alternative HTML

Hi, I am trying to create a function for a goodness-of-fit test for the Pareto Distribution for some loss data that I have. So far I have the following: function(X=OTOL) { n <- length(X)-1 #calculated the number of values (extra as 0 included) i <- 2:640 #values of i j <- 1:639 #values of i-1 Y <- (n-j+1)*((X[i])-(X[j])) #First part of GoF model Y } Where OTOL is the ordered loss

Hi, I''ve created a Rails application and am having a problem listing by a particular category in an associated table. I followed the instructions found at the O''Reilly tutorial (http://www.onlamp.com/pub/a/onlamp/2005/03/03/rails.html?page=3) but have not had success. Can anyone please outline a better set of instructions than in the "Showing Recipes in a

Hi there, I am currently working on a R programming project and got stuck. I am supposed to generate a set of possibilities of 1296 different combinations of 4 numbers, ie. 1111, 1234, 2361, (only contain 1 to 6) in a matrix form here is what I got which has not been working as it keeps coming out with the same number on the row.. The code: gl1<- gl(1296,1,length=1296, labels=1:1296,

Hi there Probably a very simple solution to this problem. I have a character vector eg c("name1","name2","name3") and I want to pull out these columns from a data.frame, converting each of these columns into factors also. Many thanks Rachel -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I am trying to plot survival curves using the following code as an example: >rs1799964.coxph<-(coxph(Surv(sassurvmonths,status)~age+stage+rs1799964_TNFA,method="efron")) >plot(rs1799964.coxph,lyt=c(1,3),xlab="Survival in Months",ylab="Proportion Surviving") I am gettingthe following error message: >Error in xy.coords(x, y, xlabel, ylabel, log) :

Hi! Whenever I save a graphic in ps/eps format generated by filled.contour or image, an undesired grid is added to it (not visible on the X11 screen). For example: postscript("volcano.eps") filled.contour(volcano,col=gray(seq(0,1,,50)),levels=seq(min(volcano),max(volcano),,50)) dev.off() Any ideia how to eliminate this grid? Thanks, Rachel [[alternative HTML version deleted]]