similar to: Getting estimates from survfit.coxph

Hi List, A friend of mine recently asked the same question as Heinz T?chler. Since I've already written the code I'd like to share with the list. # x is an object returned by "survfit"; # "smed" returns a matrix of 5 columns of # n, events, median, 0.95LCL, 0.95UCL. # The matrix returned has rownames as the # group labels (eg., treatment arms) if any. smed <-

Hey guys. I''ve just found some odd behaviour within RSpec 1.1.12 , and would like to know whether this is a bug, or I''m doing something wrong. When I give multiple return values to a stub, like this: SubtitleFile.stub!(:new).and_return @sf1, @sf2 RSpec complains: Mock ''SubtitleFile_1001'' received unexpected message :size with (no args) However, if I

Terry, I recently noticed the censor argument of survfit. For some analyses it greatly reduces the size of the resulting object, which is a nice feature. However, when combined with the id argument, only 1 prediction is made. Predictions can be made individually but I'd prefer to do them all at once if that change can be made. Chris ##################################### # CODE # create

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

I have setup a Zentyal 3.2 Samba4 AD. I am using RSAT to administer the AD server. Windows7 and Linux (CentOS) clients are able to join the AD; uid/gid of directories/files show up as per the definitions in Samba4. With Mac OS X clients, there is a problem. The OS X nodes are able to join the domain (quite a few blogs on this subject). The AD connection shows green button. The AD user is

Hello, I am trying to compute the mdeian of the survival time from the function survfit: > fit <- survfit(Surv(time, status) ~ 1) > fit Call: survfit(formula = Surv(time, status) ~ 1) records n.max n.start events median 0.95LCL 0.95UCL 111 111 111 20 NA NA NA The results is NA? the fit$surv gives values between 1 and 0.749! Am I doing this correct?

Hi all! I am trying to extract output information from the survfit function in order to generate a matrix of select output for multiple factors. Specifically, I am interested in extracting the number of events (in the output below: 106, 2, 3). The variable names represented in my function (ee) are shown below, but none of those variables correspond to the column of events as shown in the output.

Hello, I create a plot from a coxph object called fit.ads4: plot(survfit(fit.ads4)) plot is located at: https://www.dropbox.com/s/9jswrzid7mp1u62/survfit%20plot.png I also create the following survfit statistics: > print(survfit(fit.ads4),print.rmean=T) Call: survfit(formula = fit.ads4) records n.max n.start events *rmean *se(rmean) median 0.95LCL 0.95UCL 203.0

Hi, I'm sure this should be simple but I can't figure it out! I want to get the median survival calculated by the survfit function and use the value rather than just be able to print it. Something like this: library(survival) data(lung) lung.byPS = survfit(Surv (time, status) ~ ph.ecog, data=lung) # lung.byPS Call: survfit(formula = Surv(time, status) ~ ph.ecog, data = lung) 1

Hi List I have data in the following form: Gene    TFBS NUDC     PPARA(1) HNF4(20) HNF4(96) AHRARNT(104) CACBINDINGPROTEIN(149) T3R(167) HLF(191) RPA2     STAT4(57) HEB(251) TAF12     PAX3(53) YY1(92) BRCA(99) GLI(101) EIF3I     NERF(10) P300(10) TRAPPC3     HIC1(3) PAX5(17) PAX5(110) NRF1(119) HIC1(122) TRAPPC3     EGR(26) ZNF219(27) SP3(32) EGR(32) NFKAPPAB65(89) NFKAPPAB(89) RFX(121)

https://bugzilla.samba.org/show_bug.cgi?id=12527 Bug ID: 12527 Summary: Sender waits for timeout when fuzzy basis file found Product: rsync Version: 3.1.2 Hardware: All OS: All Status: NEW Severity: critical Priority: P5 Component: core Assignee: wayned at samba.org

--begin inclusion --- I am trying to extract output information from the survfit function in order to generate a matrix of select output for multiple factors. Specifically, I am interested in extracting the number of events (in the output below: 106, 2, 3). The variable names represented in my function (ee) are shown below, but none of those variables correspond to the column of events as shown

Hello R users Would somebody know how to estimate survival from a frailty cox model, using the function survfit and the argument newdata ? (or from any other way that could provide individual expected survival with standard error); Is the problem related to how the random term is included in newdata ? kfitm1 <- coxph(Surv(time,status) ~ age + sex + disease + frailty(id,

Hello: I would like to obtain probability of an event for one single patient as a function of time (from survfit.coxph) object, as I want to find what is the probability of an event say at 1 month and what is the probability of an event at 80 months and compare. So I tried the following but it fails miserably. I looked at some old posts but could not figure out the solution. Here's what I did

Hello R helpers: *My first message didn't pass trough filter so here it's again* I would like to obtain probability of an event for one single patient as a function of time (from survfit.coxph) object, as I want to find what is the probability of an event say at 1 month and what is the probability of an event at 80 months and compare. So I tried the following but it fails miserably. I

Hello everyone, Once again I'm trying to set up our Samba server to operate as a PDC, now that I've verified that I have a clean install of the Samba 3 software. However there are some difficulties with this. On my client machine when I try to set the domain name in the system control panel, I get an error message saying that the domain could not be contacted. I've set up my

Hi all, I found a discrepancy between the sum() in R and either a sum done in C or Fortran for vector of just 5 elements. The difference is very small, but this is a very small part of a much larger numerical problem in which first and second derivatives are computed numerically. This is part of a numerical method course I am teaching in which I want to compare speeds of R versus Fortran (We

Hi, Thanks in advance for any advice you can give me, I am very stumped on this problem... I use R every day and consider myself a confident user, but this seems to be an elementary problem.. Outline of problem: I am analysing the results of a study on protein expression in cancer tissues. I have raw intensities from 2 different types of cancer and normal tissue, which can be taken from several

I am unable to get the basehaz function to apply a proportional hazards model to a new data frame. I replicated my specific situation with the example for coxph in the help, where I changed the x value of the first record from 0 to 1. Is there something incorrect in the syntax that I am using? Thanks in advance! test1 <- list(time= c(4, 3,1,1,2,2,3), status=c(1,NA,1,0,1,1,0),

Dear all, Does R or S-plus or any of their packages provide Non-parametric "Run test" (which tests whether a sequence of numbers might be random or not)? If yes, i'd like a numerical illustration of this test. Any response / help / comment / suggestion will be greatly appreciated. Thanks in advance. ------------------------------- Mohammad Ehsanul Karim <wildscop at