Ah!
Sorry ... should have dug deeper into the examples section to notice that.
Thank you for the quick reply,
-steve
On Thu, Jun 2, 2016 at 8:59 AM, Frank Harrell <f.harr...@vanderbilt.edu> wrote:
> This happens when you have not strat variables in the model.
>
>
> ------------------------------
> Frank E Harrell Jr Professor and Chairman School of Medicine
>
> Department of *Biostatistics* *Vanderbilt University*
>
> On Thu, Jun 2, 2016 at 10:55 AM, Steve Lianoglou <lianoglou.st...@gene.com>
> wrote:
>
>> Hello foks,
>>
>> I'm trying to plot the number of patients at-risk by setting the
>> `n.risk` parameter to `TRUE` in the rms::survplot function, however it
>> looks as if the numbers presented in the rows for each category are
>> just summing up the total number of patients at risk in all groups for
>> each timepoint -- which is to say that the numbers are equal in each
>> category down the rows, and they don't seem to be the numbers specific
>> to each group.
>>
>> You can reproduce the observed behavior by simply running the code in
>> the Examples section of ?survplot, which I'll paste below for
>> convenience.
>>
>> Is the error between the chair and the keyboard, here, or is this perhaps
>> a bug?
>>
>> =========== code ===========
>> library(rms)
>> n <- 1000
>> set.seed(731)
>> age <- 50 + 12*rnorm(n)
>> label(age) <- "Age"
>> sex <- factor(sample(c('Male','Female'), n, rep=TRUE, prob=c(.6, .4)))
>> cens <- 15*runif(n)
>> h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))
>> dt <- -log(runif(n))/h
>> label(dt) <- 'Follow-up Time'
>> e <- ifelse(dt <= cens,1,0)
>> dt <- pmin(dt, cens)
>> units(dt) <- "Year"
>> dd <- datadist(age, sex)
>> options(datadist='dd')
>> S <- Surv(dt,e)
>>
>> f <- cph(S ~ rcs(age,4) + sex, x=TRUE, y=TRUE)
>> survplot(f, sex, n.risk=TRUE)
>> ===========
>>
>> I'm using the latest version of rms (4.5-0) running on R 3.3.0-patched.
>>
>> === Output o sessionInfo() ===
>> R version 3.3.0 Patched (2016-05-26 r70671)
>> Platform: x86_64-apple-darwin13.4.0 (64-bit)
>> Running under: OS X 10.11.4 (El Capitan)
>>
>> locale:
>> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
>>
>> attached base packages:
>> [1] stats graphics grDevices utils datasets methods base
>>
>> other attached packages:
>> [1] rms_4.5-0 SparseM_1.7 Hmisc_3.17-4 ggplot2_2.1.0
>> [5] Formula_1.2-1 survival_2.39-4 lattice_0.20-33
>>
>> loaded via a namespace (and not attached):
>> [1] Rcpp_0.12.5 cluster_2.0.4 MASS_7.3-45
>> [4] splines_3.3.0 munsell_0.4.3 colorspace_1.2-6
>> [7] multcomp_1.4-5 plyr_1.8.3 nnet_7.3-12
>> [10] grid_3.3.0 data.table_1.9.6 gtable_0.2.0
>> [13] nlme_3.1-128 quantreg_5.24 TH.data_1.0-7
>> [16] latticeExtra_0.6-28 MatrixModels_0.4-1 polspline_1.1.12
>> [19] Matrix_1.2-6 gridExtra_2.2.1 RColorBrewer_1.1-2
>> [22] codetools_0.2-14 acepack_1.3-3.3 rpart_4.1-10
>> [25] sandwich_2.3-4 scales_0.4.0 mvtnorm_1.0-5
>> [28] foreign_0.8-66 chron_2.3-47 zoo_1.7-13
>> ===========================
>>
>>
>> Thanks,
>> -steve
>>
>>
>> --
>> Steve Lianoglou
>> Computational Biologist
>> Genentech
>>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
--
Steve Lianoglou
Computational Biologist
Genentech
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
