similar to: save print survfit object to data frame

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

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

One approach is to use the rms package's cph and Mean.cph functions. Mean.cph (cph calls coxph and can compute Kaplan-Meier and other survival 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

Hello all, Here goes one of my first functions. I want to make a nonparametric multiple sample comparison with unequal sample sizes (see Zar?s Biostatistical Analysis, 3rd. Ed., pg. 201 Example 10.11, pg. 288 Example 11.10). In the real world, I want to compare samples of fish length captured with different fishing gears. After using the Kruskal-Wallis test I want to check the differences

Dear all, I'm trying to manage with user defined split function in rpart (file rpart\tests\usersplits.R in http://cran.r-project.org/src/contrib/rpart_3.1-34.tar.gz - see bottom of the email). Suppose to have the following data.frame (note that x's values are already sorted) > D y x 1 7 0.428 2 3 0.876 3 1 1.467 4 6 1.492 5 3 1.703 6 4 2.406 7 8 2.628 8 6 2.879 9 5 3.025 10 3 3.494

I've found several similar issues with write.table/read.table with Dates on this list, but trying to follow this advice I still get an error. First, I read in data from several files, constructing several date/time columns using ISOdatetime > str(Tall$Begin) 'POSIXct', format: chr [1:40114] "2005-10-02 00:00:00" "2005-10-02 00:00:00" ... > length(Tall$Begin)

https://bugzilla.netfilter.org/cgi-bin/bugzilla/show_bug.cgi?id=95 ------- Additional Comments From niels.laukens@tijd.com 2004-07-04 23:29 ------- (sorry to disturb): will this (ever) fix? is it fixed with the attachement? or won't it fix? or is there an alternative? ------- You are receiving this mail because: ------- You are on the CC list for the bug, or are watching someone who

Hello, it is likely that I will have to analyze a rather sizeable dataset: 60000 records, 10 to 15 variables. I will have to make descriptive statistics, and estimate linear models, glm's and maybe Cox proportional hazard model with time varying covariates. In theory, this is possible in R, but I would like to get some feedback on the equipment I should get for this. At this moment, I have a

Dear Sirs, I want to estimate the survival mean of a few specific teams. I'm trying to calculate it through a Kaplan Meier estimator. For doing so, I load the "survival" package and run the following instructions: "options(survfit.print.mean=TRUE)" allows showing the mean and mean standard error "KM=survfit(Surv(Dias,Censura))"

Hello everyone, I have a problem with smbmount. After I mount a certain share, I can not view files with weird charaters in their name. This happends for example with the copyright char. After searching with Google I still didn't find an answer. I guess it has something to do with a 'charset' command, but I can't figure out wich one, and what character table I should use...

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

Hi, I installed samba4 (zentyal 3.4) and I followed the Samba4 AD DC Howto from the wiki to test the linux-side. There all seems OK. Unfortunately, when I try to register a Win7 PC, the domain is not found. So what can I do to test things on the Windows-side ? I did : C:\net view /domain:ace_domain returns : \\ZENTYAL1 which is my DC After manually "mounting" the test1-share : C:\net

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

Hi all, I am back for more punishment- Anyway I reinstalled my Linux using SuSe 8.2 pro. I did not install the defualt Samba , instead I got the latest tarball for 3.0b2 I also installed webmin. Now it comes up with error "The Samba server executable /usr/sbin/smbd was not found. Either Samba is not installed on your system or your module configuration is incorrect." Now when I

Hello, I am not very familiar with regular expressions and escaping. I need to replace the %-signs in a character vector with elements as "income 0%-33%# to be replaced by "income 0\%-33\%" (for later use in LaTeX). Using gsub("%","\\%","income 0%-33%") does not give the desired result. However, gsub("%","\\\\%","income

I was reading about the project drawer feature in Textmate, which is Mac only. Is there a similar feature in a Windows based text editor that works with R. This feature sounds really useful. Thanks, Graham [[alternative HTML version deleted]]

Hello, I am afraid I do not fully understand all intricacies of programming in lattice plots. In the code below I try to plot an ROC curve, following R-news 4(1). When I condition on the variable 'group' I get the error message below, when I plot the curve for all data (i.e., y ~ pred.prob), I get the plot I want. Can someone point out why conditioning gives that message? Thanks, Ruud

I have a print method for a set of statistical tests, vcdExtra::CMHtest, for which I'd like to have more sensible printing of pvalues, as in print.anova(). [Testing this requires the latest version of vcdExtra, from R-Forge **|install.packages("vcdExtra", repos="http://R-Forge.R-project.org")|** ] With my current print method, I get results like this, but all Prob values