similar to: survival5

Anyone using the penalised partial likelihood routines in survival5 should update their version. A bug has been fixed in the S package: in coxph() models with penalised likelihood and strata it was possible in some circumstances to get an infinite loop or perhaps an incorrect answer. The new version (2.3) is on cran.r-project.org and will percolate through CRAN in the next few days. -thomas

Anyone using the penalised partial likelihood routines in survival5 should update their version. A bug has been fixed in the S package: in coxph() models with penalised likelihood and strata it was possible in some circumstances to get an infinite loop or perhaps an incorrect answer. The new version (2.3) is on cran.r-project.org and will percolate through CRAN in the next few days. -thomas

I downloaded binary of survival5 for windows95 (also the splines, date). The function survsum shows up in content, has a help page, but the function itself is missing. Can we cut and paste the survsum from survival4 ? Mai Zhou P.S. I am not picking on R, not even complaining. I think R is great and I am try to contribute my (small) part to make it even better in way of bug report (or

SEE ALSO ORIGINAL POSTING IN PR#2933 On May 6, 2003 03:58 pm, Thomas Lumley wrote: > > Looking at a wider context in the code > > pfun <- function(coef, theta, ndeath) { > if (theta == 0) > list(recenter = 0, penalty = 0, flag = TRUE) > else { > recenter <- log(mean(exp(coef))) > coef <- coef - recenter

On Tue, 6 May 2003, Jerome Asselin wrote: > > I am confused on how the log-likelihood is calculated in a parametric > survival problem with frailty. I see a contradiction in the frailty() help > file vs. the source code of frailty.gamma(), frailty.gaussian() and > frailty.t(). > > The function frailty.gaussian() appears to calculate the penalty as the > negative

survexp in survival5 doesn't seem to work for me. see below: > library(survival5) Attaching Package "package:survival5": The following object(s) are masked from package:base : sort.list > library(chron) > data(ratetables) > survexp(~ratetable(year=julian(6,1,1991), + sex=1,age=35*365.24),times=(0:30)/6*365.24) Error in as.character(as.date(c(min(R[, 3]),

> data(package=survival5) Warning: no data listing found Warning: package `survival5' contains datasets but no index The data loads OK. Mai Zhou -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Prof Brian D Ripley wrote: > On 5 Nov 1999, Peter Dalgaard BSA wrote: > > > "Richard A. Bilonick" <rab at nauticom.net> writes: > > > > > Hi. I'm using 0.65.1 under Linux (Red Hat 6.0). I installed splines and then > > > survival5. I used library(survival) and was able to create obs using Surv. But I > > ... > > > -lg2c -lm

Dear R-Helpers: What other references are there on the capabilities of the survival5 package other than the help files and the chapter on survival analysis in every edition of Modern Applied Statistics with S? I'm thinking of something like "An Introduction to Survival Analysis in R" with worked examples that might complement or extend the chapter in MASS. Thanks, Spencer

I am confused on how the log-likelihood is calculated in a parametric survival problem with frailty. I see a contradiction in the frailty() help file vs. the source code of frailty.gamma(), frailty.gaussian() and frailty.t(). The function frailty.gaussian() appears to calculate the penalty as the negative log-density of independent Gaussian variables, as one would expect: >

** High Priority ** Hello R users My issue is both theorical and technical. I would like to run a penalised discriminant analysis with the fda() function, but I don''t know all the details of splines theory. I try on the example of the phonems from the article "Penalised Discriminant Analysis" of Hastie, Buja and Tibshirani 1994 : 5 groups and 256 variables. The 256

On a PMac 7200 (running LinuxPPC 3Q 99), trying to install "R-ctest-0.9.R16-1.ppc.rpm" and "R-splines-2.0.R3-3.ppc.rpm" rpms currently in the linuxppc binary download directory results in error messages saying that files in these rpms conflict with files from R-base 1.01-1. Of course, "R-rpart" and "R-survival4" then also can't be installed because the

Hi All, I am relatively new to R (took a class 2 yrs ago ...), and am hoping someone can point me in the right direction for a problem I'd like to solve in R : I would like to calculate the relative mortality of a particular impairment, relative to the standard population. I.e. I'm trying to find S_relative(t) in the eqn below : S_impar(t) = S_standard(t) * S_relative(t) where

A couple of clarifications for you. 1. I write mixed effects Cox models as exp(X beta + Z b), beta = fixed effects coefficients and b = random effects coefficients. I'm using notation that is common in linear mixed effects models (on purpose). About 2/3 of the papers use exp(X beta)* c, i.e., pull the random effects out of the exponent. Does it make a difference? Not much: b will be

I have installed R version 1 for windows and have installed the survival5 package in the library folder, but the base package only seems to partially acknowledge survival5's presence. How have I screwed up? The diagnostics are as follows... >library(survival5) ** no warning given.. seems OK ** >help(package=survival5) **this gives the contents of the package correctly,

Hi, I have been trying to fit an un-penalised gam in mgcv (in order to get more reliable p-values for hypothesis testing), but I am struggling to get the model to fit sucessfully when I add in a te() interaction. The model I am trying to fit is: gam(count~ s(x1, bs = "ts", k = 4, fx = TRUE) + s(x2, bs = "ts", k = 4, fx = TRUE) + te(x2, x3, bs =

Patrick Connolly wrote: > According to Richard A. Bilonick: > |> > |> I recently installed Red Hat Linux 7.0. I installed Version 1.2.0 (2000-12-15) of > |> R. I installed a number of packages including chron, date, ellipse, sgeostat, > |> survival5, and xtable. the "library()" command shows that they are installed. > |> > |> chron works but I get

Hi Everyone We would like to connect our Splicecom Maximiser PBX to our Asterisk box via H323 so that we can send our US calls via a low cost carrier (e.g. Broadvoice). Has anyone managed to do this in the past (I remember seeing some companies also worked with this system in the UK). The Maximiser only speaks H323 (not SIP) and can act as an H323 Gatekeeper, so in theory we should just be

The new code seems to have broken various things. Autoloading of libraries doesn't seem to work: > library(survival4) Autoloading required library: splines Error in pos.to.env(pos) : invalid "pos" argument > search() [1] ".GlobalEnv" "library:survival4" "library:date" [4] "library:base" The coxph function the

The new code seems to have broken various things. Autoloading of libraries doesn't seem to work: > library(survival4) Autoloading required library: splines Error in pos.to.env(pos) : invalid "pos" argument > search() [1] ".GlobalEnv" "library:survival4" "library:date" [4] "library:base" The coxph function the