search for: gvhd

...plot(survfit()). The plot should, I believe, show two lines one for survival in each of two treatment (Drug) groups, however my plot shows only one line. What am I doing wrong? My code is reproduced below, my figure is attached to this EMail message. John > #Create simple survival object > GVHDdata<-list(Time=GVHD$Time,Time30=(GVHD$Time)/30, + Age=GVHD$Age,Drug=GVHD$Drug,Died=GVHD$Died, + AgeGrp=cut(GVHD$Age,breaks=c(0,15,25,45))) > > summary(GVHD$Drug) MTX MXT+CSP 32 32 > > > > fit0<-coxph(Surv(Time30,Died)~Drug,data=GVHDdata) > summary(fit0...

...in two types ('type') and the measurements were repeated a different time points ('day'). I read in the Lattice book that this can be done using the 'flowViz' package, so I tried to reproduce the example given in there and obtained an error: library('flowViz') data(GvHD, package='flowCore") densityplot(Visit~'FSC-H'|Patient, data=GvHD) Error: ' "FSC-H" ' is/are no valid parameter(s) in this frame I am not sure if I am missing something. I was hoping to apply something similar to my data, maybe like this? desityplot(id~pid|...

Dear R users: Analysis of the impact of a time-dependent covariate (GVHD or use of steroid after bone marrow transplantation) on two competing endpoints (invasive fungal infection and death) is frequently encountered in the setting of BMT data. Coxph package can be used as the following: for the analysis of GVHD: > gvhd -> coxph(Surv(start,stop,status = =1) ~...

Dear R users: Analysis of the impact of a time-dependent covariate (GVHD or use of steroid after bone marrow transplantation) on two competing endpoints (invasive fungal infection and death) is frequently encountered in the setting of BMT data. Coxph package can be used as the following: for the analysis of GVHD: > gvhd -> coxph(Surv(start,stop,status = =1) ~...

On Thu, 29 Jan 2009 20:55:19 +0100 r-help-request at r-project.org wrote: > Welcome to the R-help at r-project.org mailing list! > > To post to this list, send your email to: > > r-help at r-project.org > > General information about the mailing list is at: > > https://stat.ethz.ch/mailman/listinfo/r-help > > If you ever want to unsubscribe or change your

...yz argument. I bet there is a method to draw multiple gates through the panel-function, as curv1filter can also identify multiple peaks automatically and draw them into a densityplot... This script works for  xyplot but not for densitylot:     library(flowCore)     library(flowViz)     data(GvHD)     Filter1        <-    rectangleGate(filterId="Filter1", "FSC-H" = c(0, 200))     Filter2        <-    rectangleGate(filterId="Filter1", "FSC-H" = c(300, 400))     xyplot( `SSC-H` ~ `FSC-H` , data=GvHD[[1]],         panel = function(...) {        ...