similar to: bladder1 dataset in survival library

> On Apr 16, 2018, at 3:20 PM, David Winsemius <dwinsemius at comcast.net> wrote: > >> >> On Apr 16, 2018, at 2:58 PM, Therneau, Terry M., Ph.D. via R-devel <r-devel at r-project.org> wrote: >> >> A user asked me about this and I can't figure it out. >> >> tmt% R >> R Under development (unstable) (2018-04-09 r74565) --

A user asked me about this and I can't figure it out. tmt% R R Under development (unstable) (2018-04-09 r74565) -- "Unsuffered Consequences" Copyright (C) 2018 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) > library(survival) > data(cgd0) Warning message: In data(cgd0) : data set ?cgd0? not found ---- The data set is present and can be

Hello, I am interested in producing the expected number of events, in a recurring events setting. I am using the Andersen-Gill model, as fit by the function "coxph" in the package "survival." I need to produce expected numbers of events for a cohort, cumulatively, at several fixed times. My ultimate goal is: To fit an AG model to a reference sample, then use that fitted model

Dear users, In a study with recurrent events: My objective is to get estimates of survival (obtained through a Cox model) by rank of recurrence and by treatment group. With the following code (corresponding to a model with a global effect of the treatment=rx), I get no error and manage to obtain what I want : data<-(bladder)

Hello R-help. I am trying to make a stacked barplot where the color of the sections of each bar depend on another variable. > myData[1:11,] score percent marker cellType Malignant 1 0 100.00000 ESR1 (ER) Bladder.M(5) TRUE 2 0 80.00000 PAX8 Bladder.M(5) TRUE 3 1 20.00000 PAX8 Bladder.M(5) TRUE 4 0 100.00000 ESR1 (ER) Brain.N(3) FALSE 5 0

It would be nice if the 00check.log file also included this part of the output: Running ?bladder.R? Comparing ?bladder.Rout? to ?bladder.Rout.save? ... OK Running ?book1.R? Comparing ?book1.Rout? to ?book1.Rout.save? ... OK Running ?book2.R? Comparing ?book2.Rout? to ?book2.Rout.save? ... OK etc. The survival package has enough test scripts that it exceeds my terminal's scroll

The expected number of bladder cancer over next 20 years a tire industry is 1.8. Poission distribution is assumed to hold and 6 reported deaths are caused by bladder cancer among the employees. Trying to find how unusual this event is. > ppois(q=6, lambda=1.8, lower.tail = TRUE, log.p = FALSE) [1] 0.9974306 not sure if ppois is the right one to use and the parameters... thx much

Hi All, I am fitting a tree to censored survival data using the rpart package and wanted to better understand the results. I am trying to interpret the output from the tree. I am interested in understanding what "yval" is for a survival tree. I see in the output of summary, the phrase "estimated rate". The estimated rate is 1 for the entire tree, and more of less for each

Dear volunteer: I am a graduate student of medcine in china.And now,I am devoting myself to constructing a nomograms of bladder cancer.I want to do it with R-project.However, I do not know how to construct a nomograms with R-project.I want to get yours help,thank you! I wish you can tell me the operating procedure of the R-project. And I apologize for my english,it is poor,sorry!

Hi there, I perform a stratified cox and then I need the strata as a  numeric array "straft.ln" ft.ln <- coxph(Surv(times,deaths)~ages+chemos+chemos:f1+chemos:f2+horms+horms:f1+horms:f2+grades+grades:f1+grades:f2+positives+positives:f1+positives:f2+sizes+sizes:f1+sizes:f2+strata(stra),data=ddd) basehazzft.ln=basehaz(ft.ln,centered=FALSE) H0ft.ln=c(basehazzft.ln[,1])

Hello all: While I'm recovering from gall bladder surgery, I thought I'd do something useful (since sitting up is about the only comfortable position I can find)... like create a HOWTO Samba (-2.0.7-21ssl) on RH7 (w/ xinetd-2.1.8.9pre14). I have already sent off the request to the HOWTO coordinator at linuxdoc.org (has anyone noticed that the listed email address

dear professor: I am a doctor of urinary,and I am developing a nomogram of bladder tumor.Now I have a problem about this. I have got the result like this through analysing the dataset "exp11.sav" through multinominal logistic regression by SPSS 17.0.(the Sig. is high,that is good ,it is just aexperimental data ) Parameter Estimates Ya B Std. Error Wald df Sig. Exp(B) 95%

Dear Sir or Madam? I am a doctor of urology,and I am engaged in developing a nomogram of bladder cancer. May I ask for your help on below issue? I set up a dataset which include 317 cases. I got the Binary Logistic Regression model by SPSS.And then I try to reconstruct the model ?lrm(RECU~Complication+T.Num+T.Grade+Year+TS)? by R-Project,and try to internal validate the model through

Hi, I am currently analysis Raman spectroscopic data with the hyperSpec package. I consulted the documentation on this package and I found an example work-flow dedicated to Raman spectroscopy (see the address : http://hyperspec.r-forge.r-project.org/chondro.pdf) I am currently trying to remove outliers thanks to PCA just as they did in the documentation, but I get a message error I can't

Is there an easy way to override the compile options when installing packages from CRAN on Unix-style systems? I tried setting CFLAGS, CXXFLAGS, and FFLAGS but they are not picked up. For example (with a random package and adding -mtune only to try to get rid of the default) export CFLAGS="-march=native -mtune=native" export CXXFLAGS="-march=native -mtune=native" R CMD

Hi, I don't understand the following. When I create a small artificial set of date information in class POSIXct, I can calculate the mean and the median: a = as.POSIXct(Sys.time()) a = a + 60*0:10; a [1] "2009-03-10 11:30:16 EDT" "2009-03-10 11:31:16 EDT" "2009-03-10 11:32:16 EDT" [4] "2009-03-10 11:33:16 EDT" "2009-03-10 11:34:16

I don't understand what is happening! I have a (large) object sim1, an matrix list with dim c(101,101) where each element is an 3*3 matrix. I am subsetting that with a matrix coo, of dim c(100,2), of unique indices, but the resulting object has length 99, not 100 as expected. Code reproducing the problem follows: library(RandomFields) set.seed(123) sim0 <- GaussRF(x=seq(0, 100, by=1),

Hello, a. when I use hclust with the methods media, centroid, and mcquitty, and plot the results, the dendrograms have lines that are crossing each other. Is this ok? b. My next question refers to plot.dendrogram: How can I use parameters as "hang" or "cex" here? E.g. for st <- as.dendrogram(subtreeshc[[x]]) I would like to have something like this, where cex and hang

Hello! Can I use DBM databases (as they are used in Perl) for the import of large datasets in R? Thanks in advance for your help, Petra Steiner - --------------------------------------------------- Petra Steiner Arbeitsbereich Linguistik Universitaet Muenster Huefferstrasse 27 48149 Muenster

Hi all, I've read the R for windows FAQ and am a little confused re: memory.limit and memory.size to start using R 2.6.2 on WinXP, 2GB RAM, I have the command line "sdi --max-mem-size=2047M" Once the Rgui is open, memory.limit() returns 2047, memory.size() returns 11.315, and memory.size(max=T) returns 19.615 Shouldn't memory.size(max=T) return 2047? Upon running several