search for: rmean

Hi, When I run the survfit function, I want to get the restricted mean value and the standard error also. I found out using the "print" function to do so, as shown below, print(km.fit,print.rmean=TRUE) Call: survfit(formula = Surv(diff, status) ~ 1, type = "kaplan-meier") records n.max n.start events *rmean *se(rmean) median 200.000 200.000 200.000 129.000 0.145 0.237 1.158 0.95LCL 0.95UCL 0.450 1.730 * re...

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 100.0 100.0 103.0 486.7 24.4 387.0 340.0 467.0 * restricted mean with upper limit = 819 Questions: 1. What is...

...rvival estimates) can compute mean restricted life. Frank Dinesh W wrote > I am using survfit to generate a survival curve. My population is such > that my x axis is in days and i have a starting population of say 10,000 > of which say only 2000 are left as of day 365. When I try to print rmean > it does not print and in any case I am not interested in that. I actually > want to sum up all the daily survival values starting at 1.0 (S_0) for day > 1 through 0.2 (S_365) through day 365. I then want to assume "r% retained > each day" and compute the remaining integral a...

Hi, I'm not sure if this is the proper way to ask questions, sorry if not. But here's my problem: I'm trying to do a bootstrap estimate of the mean for some survival data. Is there a way to specifically call upon the rmean value, in order to store it in an object? I've used print(...,print.rmean=T) to print the summary of survfit, but I'm not sure how to access only rmean because it does not show up under attributes for survfit. Thanks for any help in advance! Fayaaz Khatri [[alternative HTML version dele...

...#39;<0.0001' . > CMHtest(MSPatients[,,1]) Cochran-Mantel-Haenszel Statistics for New Orleans Neurologist by Winnipeg Neurologist AltHypothesis Chisq Df Prob cor Nonzero correlation 51.424 1 7.4426e-13 cmeans Col mean scores differ 55.393 3 5.6601e-12 rmeans Row mean scores differ 53.631 3 1.3450e-11 general General association 64.318 9 1.9580e-10 In the print.CMHtest() function below the lines before # # TODO give the output shown above. The lines below try to use print.anova(), but this gives something even worse: Cochran-Mantel-Haenszel S...

...Cat=VarCat),FUN=length) RLevCat <- aggregate(Rank,list(LevCat=VarCat),FUN=sum) MLevCat <- aggregate(Rank,list(LevCat=VarCat),FUN=mean) SampleSummary <- data.frame(LevCat,RLevCat[,2],NLevCat[,2],MLevCat[,2]) names(SampleSummary)<-c("Samples","RSum","N","RMean") SampleSummary <- SampleSummary[order(SampleSummary$RMean,decreasing=T),] NCat <- length(LevCat) N <- length(dat.multcomp$VarCat) Results <- data.frame(rep(NA,6),rep(NA,6),rep(NA,6),rep(NA,6)) names(Results) <- c("Comparison","Difference","SE",&quo...

Hello, I have estimated a survival model with six strata: >model.b <- survfit(Surv(time=start.tijd,time2=eind.tijd2,event=va)~strata(product.code) , data=wu.wide) I would like to save the output of >print(model.b,print.n="records",show.rmean=FALSE) in a dataframe so that I can export it later. How do I do this? Note that summary(model.b) gives an error: Error in as.matrix(x) : (subscript) logical subscript too long (R 1.9.1, windows xp, sp 2) Thanks for any help, Ruud

Dear r-help, I have a matrix, suppose, 10x10, and I need the matrix 5x5, having in each cell a mean value of the cells from the initial matrix. Please, point me to a function in R, which can help me doing that. Digging the documentation and mail archives didn't give me a result. Thank you. --- Best regards, Vladimir mailto:wl at eimb.ru

...1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 temp2 <- function(y, wt, x, parms, continuous) { # Center y n <- length(y) y <- y- sum(y*wt)/sum(wt) if (continuous) { # continuous x variable temp <- cumsum(y*wt)[-n] left.wt <- cumsum(wt)[-n] right.wt <- sum(wt) - left.wt lmean <- temp/left.wt rmean <- -temp/right.wt goodness <- (left.wt*lmean^2 + right.wt*rmean^2)/sum(wt*y^2) list(goodness= goodness, direction=sign(lmean)) } } Paolo Radaelli Dipartimento di Metodi Quantitativi per le Scienze Economiche ed Aziendali Facolt? di Economia Universit? degli Studi di Milano-Bicocca P.zza dell...

...quot;KM" shows the model It is in this very moment when it shows the mean value and the mean standard error. I'd like to know if there is any instruction that returns the mean value and its standard error, something like: "KM$rmean" Thank you very much in advanced. Sincerely yours, Alejandro G. Lobo

...In addition: Warning message: cannot open file `C:/PROGRA~1/R/rw1091/library/marray/swirldata/81' > fnames <- as.vector(swirlTargets@maInfo[,1]) > swirl.raw <- read.marrayRaw(fnames, path = datadir, + name.Gf = "Gmean", name.Gb = "morphG", + name.Rf = "Rmean", name.Rb = "morphR", + layout = swirl.layout, + gnames = swirl.gnames, + targets = swirlTargets) Error in file(con, "r") : unable to open connection In addition: Warning message: cannot open file `C:/PROGRA~1/R/rw1091/library/marray/swirldata/81' > array1.r...

Dear List, Does anybody no how to compare mean survival times for two (more) groups in R? What test statistics should I use? Thank you very much! Joe [[alternative HTML version deleted]]

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

Hi all, Using: R version 2.8.1 Patched (2009-03-07 r48068) on OSX (10.5.6) with survival version: Version: 2.35-3 Date: 2009-02-10 I get the following using the first example in ?summary.survfit: > summary( survfit( Surv(futime, fustat)~1, data=ovarian)) Call: survfit(formula = Surv(futime, fustat) ~ 1, data = ovarian) time n.risk n.event survival

Hello, I need some results from the survival analysis of my data that I do not know whether exist in Survival Package or how to obtain if they do: 1. The Mean survival time 2. The standard error of the mean 3. Point and 95% Lower & Upper Confidence Intervals estimates Any help will be greatly appreciated. Cem [[alternative HTML version