search for: psline

Hi all, I have some questions about the predicted HR in coxph function including psline covariate. If just fitting covariate as linear form, by default, the reference value for each of predictions type (linear predictor, risk and terms) is the mean covariate within strata. If the psline is specified for the covariate, what's the reference value for the predictions? I did some te...

...ot explain. Framework: Ubuntu Linux, R2.11. For testing survival I have a separate directory and Makefile. I pull everything into the local .RData, no packages, library, or namespace. (It's easier to add test modifications to a routine in a chain of calls). A test of survreg + psline would fail because splines is not automatically loaded in this scenario. In my console window type > library(splines) > fit <- survreg(Surv(time, status) ~ ph.ecog + spline(age), data=lung) as the first 2 lines a a new session (the lung data is also saved) and up would pop an error...

...) above, but with the same model *convoluted with a gaussian* (whose height and width should be fitted). R looks really powerful for doing these kind of things, and I'd like to use it for analyzing X-ray spectra. After doing some research, I came accross a few packages: gss, stats, assist, psline; but I have not yet been able to figure out a way of performing the fits with splines described in (2) and (3) above. Any help would be welcome! EOL [new to R, and not a statistician--only a physicist] -- Dr. Eric-Olivier LE BIGOT (EOL) CNRS Associate Researcher ~~~o~o~~~~~~...

Hi there, I have a question on how to extract the linear term in the penalized spline in coxph. Here is a sample code: n=100 set.seed(1) x=runif(100) f1 = cos(2*pi*x) hazard = exp(f1) T = 0 for (i in 1:100) { T[i] = rexp(1,hazard[i]) } C = runif(n)*4 cen = T<=C y = T*(cen) + C*(1-cen) data.tr=cbind(y,cen,x) fit=coxph(Surv(data.tr[,1],

...n the search path; the wrong function could be picked up. Second, a model with survival::strata(inst) in the formula would not do what users expect. These have now been addressed. 1. For the coxph, concordance, survcheck, survfit, and survdiff modeling functions, any Surv(), strata(), cluster(), psline(), tt(), ridge(), frailty(), frailty.gamma(), frailty.gauss(), or frailty.t() model terms will use the versions from the survival namespace, automatically. A survival:: modifier is not required. 2. Any use of survival::strata(), survival::cluster(), survival::pspline(), ... (all but survival::Sur...