similar to: survreg (survival) reports erroneous results for left-censored (PR#2293)

On Wed, 13 Nov 2002, Jan de Leeuw wrote: > > No problemo. And, in fact, I get the same results in > the R-1.6.0 Carbon version. I don't. Could there be a G3/G4 issue? -thomas > --- Jan > > On Wednesday, November 13, 2002, at 02:05 PM, tim@timcohn.com wrote: > > > Full_Name: Tim Cohn > > Version: 1.6.1 > > OS: Macintosh OS X > > Submission

Full_Name: Tim Cohn Version: 1.6.1 OS: Macintosh OS X Submission from: (NULL) (130.11.34.250) The Mac version of survreg does not handle left-censored data correctly (at least the results are not what I get doing it other ways, and they are not the same as I get running R 1.6.1 in Windows 98se; the Windows 98 results are correct). On the windows version of R 1.6.1. >

I have build R-1.6.0 on Jaguar (actually on 10.2.2), using gcc 3.1 from Apple, g77 3.1 from fink, ATLAS from fink. It passes all checks OK. I also re-compiled the 233 packages listed below with the same setup. I have NOT tested all these packages in any detail. A tar.gz file with the whole /usr/local/lib/R will be on ftp://gifi.stat.ucla.edu on Tuesday. It will NOT include anything from fink

This is probably only interesting to zero or one or two persons. gcc-3.3 from Apple's CVS does not handle the -bundle flag properly. I had to hack configure.ac, which now has shlib_ldflags="-nostdlib /usr/lib/bundle1.o -Wl,-bundle -Wl,-bundle_loader,/usr/local/lib/R/bin/R.bin -lgcc" ;; LIBR_LDFLAGS="-dynamiclib -install_name \$(Rexeclibdir)/libR.dylib -framework

I don't think it is necessary to put binary versions of the Darwin/X11 version on CRAN any more. The current sources compile without problems on both Mac OS X 10.1 and 10.2, using either gcc/g77 2.95.2 or gcc/g77 3.1, or f2c. It seems to me that if you are using R from the command line, using X11, that you are sophisticated enough to install the tools below. And, if you want to install and

Can anybody give me a neat example of interval censored data analysis codes in R? Given that suvreg(Surv(c(1,1,NA,3),c(2,NA,2,3),type="interval2")~1) works why does survreg(Surv(data[,1],data[,2],type="interval2")~1) not work where data is : T.1 T.2 Status 1 0.0000000 0.62873036 1 2 0.0000000 2.07039068 1 3 0.0000000

Hello, I'm working with some left censored survival data using accelerated failure time models. I am interested in fitting different distributions to the data but seem to be getting the same results from the model fit using survreg regardless of the assumed distribution. These two codes seem to provide the same results: aft.gaussian <-

I am trying to fit a lognormal distribution on interval-censored data. Some of my intervals have a lower bound of zero. Unfortunately, it seems like survreg() cannot deal with lower bounds of zero, despite the fact that plnorm(0)==0 and pnorm(-Inf)==0 are well defined. Below is a short example to reproduce the problem. Does anyone know why survreg() must behave that way? Is there an alternate

Can survreg() handle interval-censored data like the documentation says? I ask because the command: survreg(Surv(start, stop, event) ~ 1, data = heart) fails with the error message Invalid survival type yet the documentation for Surv() states: "Presently, the only methods allowing interval censored data are the parametric models computed by 'survreg'"

Byron I meant there is no native (Aqua) version of GTK yet (there is wxWindows and Qt and Tk, so possibly GTK will follow). The X11 version of RGtk works for me too. I agree with your diagnosis of the GUI problem. Best -- Jan On Monday, August 19, 2002, at 12:08 AM, Byron Ellis wrote: > > On Monday, August 19, 2002, at 12:09 AM, Jan de Leeuw wrote: >> > > I accidentally

homals-0.9.0 is on CRAN -- by Jan de Leeuw and Patrick Mair This package implements the methods discussed in Gifi, Nonlinear Multivariate Analysis, Wiley, 1990. In the Gifi terminology it covers homals, princals, canals, morals, criminals, and overals. The R implementation fills several gaps in Gifi, adding multiple ordinal, numerical, and polynomial data transformations. Differences with

homals-0.9.0 is on CRAN -- by Jan de Leeuw and Patrick Mair This package implements the methods discussed in Gifi, Nonlinear Multivariate Analysis, Wiley, 1990. In the Gifi terminology it covers homals, princals, canals, morals, criminals, and overals. The R implementation fills several gaps in Gifi, adding multiple ordinal, numerical, and polynomial data transformations. Differences with

-------------- next part -------------- An embedded message was scrubbed... From: Terry Therneau <therneau at mayo.edu> Subject: Re: Interval censored Data in survreg() with zero values! Date: Tue, 30 Dec 2008 16:46:37 -0600 (CST) Size: 4268 URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090102/abe75d02/attachment-0002.eml>

The Journal of Statistical Software was founded by Jan de Leeuw in 1996. Currently, there are 26 volumes with 260 aricles and the journal publishes about 4 new volumes per year. The web address is www.jstatsoft.org. All articles are pdf files, since 2005 all produced from JSS LaTeX templates. JSS is edited by Jan de Leeuw, since 2005 jointly with Achim Zeileis of the Wirtschaftsuniversität Wien.

Pretty soon, hopefully, the RAqua version will make the Darwin/X11 version unnecessary. In the meantime, for convergence, I'll modify the Gifi version in various ways. -- It's 1.7.1 -- It no longer supports gnome -- It still uses Tcl/Tk for X11 (using 8.5 from CVS) !! It no longer uses anything from fink (readline and dlcompat as on Stefano's site, jpeg and png and teTeX from

Hello, I have interval censored data, censored between (0, 100). I used the tobit function in the AER package which in turn backs on survreg. Actually I'm struggling with the distribution. Data is asymmetrically distributed, so first choice would be a Weibull distribution. Unfortunately the Weibull doesn't allow for zero values in time data, as it requires x > 0. So I tried the

With the current version from rsync (8/16, 19:00) I get make[3]: *** No rule to make target `helpsearch.c', needed by `Makefile'. Stop. make[2]: *** [R] Error 1 make[1]: *** [R] Error 1 make: *** [R] Error 1 === Jan de Leeuw; Professor and Chair, UCLA Department of Statistics; Editor: Journal of Multivariate Analysis, Journal of Statistical Software US mail: 9432 Boelter Hall, Box

RAqua now builds fine, but for some reason tcl/tk support cannot be compiled in. I get checking for /Library/Frameworks/Tcl.framework/tclConfig.sh... /Library/Frameworks/Tcl.framework/tclConfig.sh checking for /Library/Frameworks/Tk.framework/tkConfig.sh... /Library/Frameworks/Tk.framework/tkConfig.sh checking whether compiling/linking Tcl/Tk code works... no This is with configure

R-devel has the --with-lapack flag for configure. If you build R with --with-blas="-framework vecLib" --with-lapack="--framework vecLib" then the build goes through, using the native optimized BLAS and Lapack in /System/Library/Frameworks/vecLib.framework. Also, it works, in the sense that it does eigenvalue problems correctly. I don't have any timings yet, because I

1. R-devel on OS X aborts when the code below is run after spitting out the message Loading required package: MASS Parameter 5 to routine DSTEIN was incorrect Mac OS BLAS parameter error in DSTEIN, parameter #0, (unavailable), is 0 2. I compiled R-patched using the included BLAS and Lapack, although Brian says this is impossible on OS X, because they are in FORTRAN (??). This version of R