similar to: Survfit plots in trellis graphics

Dear all, I am trying to produce survfit plots in a trellis environment and I would like the plots to be logarithmic. 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)

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

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

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)

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

I need to use a red-black tree, which package provides that data structure? -- Best regards Rune Schjellerup Philosof Ph.d-stipendiat, Forskningsenheden for Biostatistik Telefon: 6550 3607 E-mail: rphilosof at health.sdu.dk Adresse: J.B. Winsl?wsvej 9, 5000 Odense C SYDDANSK UNIVERSITET _______________________________________________________________ * Campusvej 55 * 5230 * Odense M * 6550

How come the k is 3 in all of this output? I expected it to be equal to r. > tmp3 <- lapply(1:3, function(k) function(r) print(paste(r, "<- r | k ->", k))) > for (i in 1:3) { tmp3[[i]](i) } [1] "1 <- r | k -> 3" [1] "2 <- r | k -> 3" [1] "3 <- r | k -> 3" -- Med venlig hilsen Rune Schjellerup Philosof Ph.d-stipendiat,

I often find myself making lists similar to this list(var1=var1, var2=var2) It doesn't seem list has an option, to make it use the name of the variable as name in the list. Is there another function that does this? -- Med venlig hilsen Rune Schjellerup Philosof Ph.d-stipendiat, Forskningsenheden for Biostatistik Telefon: 6550 3607 E-mail: rphilosof at health.sdu.dk Adresse: J.B.

Firstly, I am running beta2 as AD member. Auth through Winbind. Upgrading from beta1 to beta2 using RPM on RH 9 I noticed that Winbind will not start and smb will not stop using Init scripts. Seems to me RH9 /etc/init.d/winbind need the "winbind uid" parameter ? Here is the /var/log/samba/log.winbind output. [2003/07/04 13:54:52, 1] nsswitch/winbindd.c:main(811) winbindd version

Hello R-help list,   I am looking for way to toggle the case of the characters like a flip-flop; that is from ''Hello'' to "hELLO" or vice versa.   I know that there are a number of functions like casefold, tolower, toupper, etc. but these functions change the case in an uniform way.   Thanks in advance,   Antonio Rivero Ostoic         Antonio Rivero Ostoic PhD Student,

Browse[1]> j <- c(1,2,NA) Browse[1]> j[j==1][-1] [1] NA Browse[1]> j[j==1][-2] [1] 1 Browse[1]> j[j==1][-2] <- 2 Error during wrapup: NAs are not allowed in subscripted assignments As far as I can see, I have no NA in the lhs (not after the second subscript anyway). Besides, I have a single value on the rhs, so it should be allowed to have NAs in the lhs, according to

When I execute mclapply it creates the needed processes, but these processes never begin computing anything, they just wait indefinitely. I recently upgraded to version 2.9.1, which might have caused the problem. -- Med venlig hilsen Rune Schjellerup Philosof Ph.d.-studerende, Statistik, IST, SDU Telefon: 6550 3607 E-mail: rphilosof at health.sdu.dk Adresse: J.B. Winsl?wsvej 9, 5000 Odense

Hi I want to use a function (update) that in its body uses eval(call, parent.frame()) I would like to use this function in a function that does not contain the variables referred to in 'call'. Those variables are instead in the parent.frame() of my function (named 'second' below) Like this: a <- 2 evalu <- function(obj) { call <- obj$call eval(call, parent.frame()) }

Dear R users, I have spent most of this day figuring out how to read STATA data into R (which eventually worked) and to run a simple OLS regression. It seems that the manuals are written in the most general and abstract way which does not really make it easy to understand what's going on. I'd be glad if somebody could save my day by: telling me how to define a matrix from the

I am using the coxph, survfit and summary.survfit functions to calculate an estimate of predicted survival with confidence interval for future patients based on the survival distribution of an existing cohort of subjects. I am trying to understand the calculation and interpretation of the std.err and confidence intervals printed by the summary.survfit function. Using the default confidence

I am confused when trying the function survfit. my question is: what does the survival curve given by plot.survfit mean? is it the survival curve with different covariates at different points? or just the baseline survival curve? for example, I run the following code and get the survival curve #### library(survival) fit<-coxph(Surv(futime,fustat)~resid.ds+rx+ecog.ps,data=ovarian)