Re: [R] basehaz() in package 'Survival' and warnings() with coxph

2012-08-10 Thread hazbro
Okay the data sets dat1 and dat 2 are the same dat1 just has fewer covariates. David, I understand your concern with the number of events and number of variables I am using however. 611 is only the unique times at which the events occur where as there are 6987 events in my data of 77272 observatio

Re: [R] basehaz() in package survival and warnings with coxph

2012-08-10 Thread Terry Therneau
Since fit3.1 and fit2 are based on different data sets, why would I expect the same number of events? Also, when you have a large number of variables, are observations being deleted due to missing values? And to echo David W's comments -- it is hard for me to imagine a data set where this many

Re: [R] basehaz() in package 'Survival' and warnings() with coxph

2012-08-09 Thread David Winsemius
On Aug 9, 2012, at 5:53 PM, hazbro wrote: My sessionInfo is as follows: R version 2.15.1 (2012-06-22) Platform: x86_64-unknown-linux-gnu (64-bit) snip It will be difficult to reproduce an example here as the data set I am using in very large. I can give you an example: fit3.1<- coxph

Re: [R] basehaz() in package 'Survival' and warnings() with coxph

2012-08-09 Thread hazbro
My sessionInfo is as follows: R version 2.15.1 (2012-06-22) Platform: x86_64-unknown-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_GB.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_GB.UTF-8LC_COLLATE=en_GB.UTF-8 [5] LC_MONETARY=en_GB.UTF-8LC_MESSAGES=en_GB.UTF-8 [7] LC_PAPER=C

Re: [R] basehaz() in package survival and warnings with coxph

2012-08-09 Thread Terry Therneau
I've never seen this, and have no idea how to reproduce it. For resloution you are going to have to give me a working example of the failure. Also, per the posting guide, what is your sessionInfo()? Terry Therneau On 08/09/2012 04:11 AM, r-help-requ...@r-project.org wrote: I have a couple of

[R] basehaz() in package 'Survival' and warnings() with coxph

2012-08-08 Thread hazbro
Hello, I have a couple of questions with regards to fitting a coxph model to a data set in R: I have a very large dataset and wanted to get the baseline hazard using the basehaz() function in the package : 'survival'. If I use all the covariates then the output from basehaz(fit), where fit is a m