similar to: Passing arguments to panels in trellis plots

Dear all, is there a way to produce survfit plots in a trellis environment? I am trying this: print(Ecdf(~time | size*type, groups=alg,data=B,subscripts=TRUE, panel=function(x,groups,subscripts) { t <- survfit(Surv(time[subscripts],event[subscripts])~groups[subscripts],data=B,conf.type="none")

Dear list, I would like to plot the empirical cumulative distribution of the time needed by a treatment to attain a certain goal. A number of experiments is run with a strict time limit. In some experiments the goal is attained before the time limit, in other experiments time expires before the goal is attained. The situation is very similar to survivial analysis with censored data. I tryed

hi, I have to configure Xen server and DomU''s which should be able to talk to different vlans - depending on services put on DomU''s. I would like to enslave 2 NICs via bond interface. Both NICs are on trunk ports which are configured in the same way, the same native vlan. I don''t know, if I have to add vlans tagging to eth* interfaces then enslave tagged eth*

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)

What is the best practice in creating PVm domU''s - with or without swap? >From the admin point of view, especially in HA environment it is temptating not to have a seperate swap partition. Could it be performance issue if swapping to swap-file insted of swap-partition? What about pure swapless- no swapfile, no swap-partition? Any expirence out there? regards Longina -- --

Hi again, Bonding seems to work with mode =4. I have a strange problem now - can''t get rid of old virtual interfaces which were created on the first start up of xend. Even if I boot machine with xend disabled, I still have old stuff coming up: veth0-veth5, and vif0.0-vif0.3. I would like to configure all networkig manualy first - how can I remove those interfaces, where are they

Hi, I just installed Xen in Ubuntu Hardy. DomU is not accessible either from Dom0 or external network. I used default set up network-bridge/vif-bridge from xend-conf.sxp. Bridge has got default name eth1 which is my default ethernet interface. Before creating DomU: root@jaguar-0:/home/longina# xm list Name ID Mem VCPUs State Time(s) Domain-0

Hi all, My patch for supporting "Basic" HTTP Proxy Authentication has been updated for rsync 2.5.6 by Martin Ehmsen. The patch changes the interpretation of the RSYNC_PROXY environment variable so that the syntax user:pass@proxy.foo.com:port may be used instead of just proxy.foo.com:port (the old syntax is, of course, still supported). The updated patch is available at:

Hi all, I have put together a patch for supporting "Basic" HTTP Proxy Authentication. The patch changes the interpretation of the RSYNC_PROXY environment variable so that the syntax user:pass@proxy.foo:port may be used instead of just proxy.foo:port (the old syntax is of course still supported). The patch has only been tested lightly, but it should(TM) work. The patch (and

Dear useRs, What is the syntax to obtain survival curves for single strata on many subjects? I have a model based on Surv(time,response) object, so there is a single row per subject and no start,stop and no switching of strata. The newdata has many subjects and each subject has a strata and the survival based on the subject risk and the subject strata is needed. If I do newpred <-

This works ok: > cox = coxph(surv ~ bucket*(today + accor + both) + activity, data = data) > fit = survfit(cox, newdata=data[1:100,]) but using strata leads to problems: > cox.s = coxph(surv ~ bucket*(today + accor + both) + strata(activity), > data = data) > fit.s = survfit(cox.s, newdata=data[1:100,]) Error in model.frame.default(data = data[1:100, ], formula = ~bucket + :

Good afternoon, I am currently trying to do some work on survival analysis. - I hope to seek your advice re: 2 questions (1 general and 1 specific) (1) I'm trying to do a stratified Cox analysis and subsequently plot(survfit(object)). It seems to work for some strata, but not for others. I have tumor grade, which is a range of 1 - 4. When I divide this range of 1:4 into 2 groups, it

Dear all, I'm having difficulty getting access to data generated by survfit and print.survfit when they are using with a Cox model (survfit.coxph). I would like to programmatically access the median survival time for each strata together with the 95% confidence interval. I can get it on screen, but can't get to it algorithmically. I found myself examining the source of print.survfit to

Hello, I'm using plot.survfit to plot cumulative incidence of an event. Essentially, my code boils down to: cox <-coxph(Surv(EVINF,STATUS) ~ strata(TREAT) + covariates, data=dat) surv <- survfit(cox) plot(surv,mark.time=F,fun="event") Follow-up time extends to 54 weeks, but the last event occurs at week 30, and no more people are censored in between. Is there a

Dear users, In a study with recurrent events: My objective is to get estimates of survival (obtained through a Cox model) by rank of recurrence and by treatment group. With the following code (corresponding to a model with a global effect of the treatment=rx), I get no error and manage to obtain what I want : data<-(bladder)

hie when i use plot.survfit to plot more than one graph why I only see the last graph how do i see the other graphs.for example n=20 n1=n/2 n2=n/4 a11=4;a12=4 ;a21=4 ;a22=4 t1<-array(1,c(n1)) t2<-array(2,c(n1)) treatgrp=matrix(c(t1,t2))

I've a data set with 60000 rows of data representing 6000+ distinct loans. I did a coxph() regression on it (see call below), but a subsequent survfit() call on the coxph object is almost certainly wrong. It gives n=6 when it should be more like 6000+ (I think) > survfit(resultag) Call: survfit.coxph(object = resultag) n events median 0.95LCL 0.95UCL 6 489 Inf

Hi guys. DomU config disk = [ ''phy:/dev/xen00/debian5-base,sda1,w'', ''phy:/dev/xen00/debian5-swap,sda2,w'', ] ----------------------- mkswap /dev/xen00/debian5-swap But in domU in booting: Activating swap...failed. What wrong? -- Best Regards, alex.faq8@gmail.com

Hi! I came across R just a few days ago since I was looking for a toolbox for cox-regression. I?ve read "Cox Proportional-Hazards Regression for Survival Data Appendix to An R and S-PLUS Companion to Applied Regression" from John Fox. As described therein plotting survival-functions works well (plot(survfit(model))). But I?d like to do some manipulation with the survival-functions

I have a model with quant vars only and the error message does not make sense: (mod1 <- coxph(Surv(time=strt,time2=stp,event=(resp==1))~ +incpost+I(amt/1e5)+rate+strata(termfac), subset=dt<"2010-08-30", data=inc,method="efron")) Call: coxph(formula = Surv(time = strt, time2 = stp, event = (resp == 1)) ~ +incpost + I(amt/1e+05) + rate + strata(termfac),