similar to: Schoenfeld Residuals with tied data

hi sorry if this has been discussed before, but I'm wondering why the scaled Schoenfeld residuals do not follow the defining formula for obtaining them from the ordinary Schoenfeld residuals, but are instead offset by the estimated parameter values. e.g. library(survival) attach(ovarian) sv<-Surv(futime,fustat) f1<-coxph(sv~age+ecog.ps) f1

Hi, I have a coxph model like coxph(Surv(start, stop, censor) ~ x + y, mydata) I would like to calculate the Schoenfeld residuals for the null, i.e the same model where the beta hat vector (in practical terms, the coeff vector spat out by summary()) is constrained to be all 0s --all lese stays the same. I could calculate it by hand, but I was wondering if there is a way of doing it with

Dear all, I use the Andersen plot to check for proportional hazards assumption for a factor (say x) in the Cox regression model and obtained a straight line that pass through the origin. However, the formal test done by the R-function cox.zph, which is based on the plot of Schonefeld residuals against time, indicates that proportional hazards assumption is violated. Further, a plot of the score

Dear All, Sorry to bother you again. I have a model: coxfita=coxph(Surv(rem.Remtime/365,rem.Rcens)~all.sex,data=nearma) and I'm trying to do a plot of Schoenfeld residuals using the code: plot(cox.zph(coxfita)) abline(h=0,lty=3) The error message I get is: Error in plot.window(...) : need finite 'ylim' values In addition: Warning messages: 1: In sqrt(x$var[i, i] * seval) : NaNs

I'm useing Samba (2.0.7) on an Compaq TRUE64 (Digital OSF) Server. Everything is working fine, but the dfree-script does not work. I think it's not even called. The corosponding entry in my smb.conf is: [global] server string = Samba %v Server, xxxx workgroup = xxxx keep alive = 30 dead time = 15 dfree command = /usr/local/samba/bin/dfree The script it self works fine

Hi there, slowly but sure I'm getting absolutely insane with the printer setup on a Samba PDC. I got Cups to work. It works like a charm if someone installs a printer via the cups interface. It also prints in a good quality. On Samba side I even have those "Printers and fax devices" in the visible shares. And if i install CUPS6 Drivers manual on the target system i can even print on

A time series text with a title that seems designed to hide its wide scope is: Forecasting with Exponential Smoothing The State Space Approach Hyndman, R.J., Koehler, A.B., Ord, J.K., Snyder, R.D. Springer 2009. This book is actually an excellent overview of time series theory, ARIMA as well as state space. It is of course, in part, a manual for the forecast and other packages in what has been

Hi My card is working, but when I reboot the machine, most of the times it is not working, I get "ztcfg: ZT_CHANCONFIG failed on channel 1: No such device or address (6)" To make it work again I have to shut down, remove the card, reboot so kudzu will remove the config. shut down again, put the card back in, reboot, now kudzu see it, I choose "Ignore" and then it's

Hello, I got error message when applying the plot function to the cox.zph object to create the Schoenfeld residual plots. > plot(zph.revasFit[1]) Error in plot.window(xlim, ylim, log, asp, ...) : need finite 'ylim' values In addition: Warning messages: 1: NaNs produced in: sqrt(x$var[i, i] * seval) 2: no non-missing arguments to min; returning Inf in: min(x) 3: no

To enhance my understanding, and that of my students, I have a question about cox.zph in the survival package. If I have correctly gleaned the high-level point from the 1994 Biometrika paper of Grambsch and Therneau, it looks to me like cox.zph provides a mechanism to test for a simple trend in plots of a function of time, g(t) versus the scaled schoenfeld residuals and it also provides some

Hi, I have just started to setup and configure my asterisk box and am having trouble with it. I have a dlink nic in the box as well as the digium card. When the nic is in there by itself it works, but when I put the tdm400p in, there seems to be some sort of conflict with the network card. The tdm appears as another type of network card, which I don't think it should, should it? It shows

I am trying to use R functions from my MATLAB code using a DCOM server. I am working under Windows XP. While I mange to run basic functions and run them using scripts, when I try a more complicated set of functions I fail. I am trying to perform survival analysis, using Surv Surv(recruTime,recur) Surv creates a somewhat complex output and I get the response Problem evaluating command

DeaR users. <framework> These days i'm working on fitting an extended Cox model with time-dependent covariables and possibly time-varying effects. My data are in counting process format as described in Therneau&Grambsh's `Modeling Survival Data', page 68. I'm trying to follow Harrell's `Regression Modeling Strategies' advices for the choice of my model. This

Hi all, I was wondering if anyone had any pointers on how to determine whether or not a module has gone wonky on the TDM400P? I have a 2 FXO (channels 3 and 4) and 2 FXS unit (channels 1 and 2). The bad (?) module in question is the FXO module on channel 3. I can't dial in to or out of that channel; dialing in gives a busy signal, dialing out just shows * hanging around after attempting a

I am running * in a development environment, adding functionality as I go. The * box has a X100P card in it which ztcfg enabled as channel 1 with fxsks signalling (fxsks=1). Everything worked fine and I was able to make inbound and outbound calls to/from the PSTN, the only issue being that some exchanges wouldn't handle the dtmf signalling, but I put that down to a peculiarity with some AU

Dear R users, I am doing Cox regression using coxph(Survival) or cph(Design). I have time varying effects (diagnosed with schoenfeld residuals Chi2 test and graph) so I first want to split time into 2 separate intervals : t<6months and t>=6months, to estimate one hazard ratio (hr) for each interval. I am analysing Overall survival according to 3 prognostic factors (age,deep,ldh).

Hi Folks, Hopefully this question has enough R and not too much stats to be appropriate for this list. Based on,* Hosmer et al. 2013. Logistic regression for matched case-control studies. Applied Logistic Regression *(eqtn. 7.8)*, *I am assessing GOF of conditional (or matched) logistic regression models with the *standardized Pearson residuals*. The authors define ?large? as delta chi-squared

Hi all, I have just purchased the DevKit from Digium and received a X100P and a TDM400P (it has one FXS module). The problem is that I can't get the kernel module (wcfxs) to load and run. I have searched the archives and can't find anything about this. Do the messages below ring any bells with anyone ? There is no interrupt clash (see pci list below). I am running on a AMD 1.2G

Hi Group, I have a problem with my FXO interface. It seems that asterisk cannot see any configured channels, only pseudo. Pls see the error and my config below. Rgrds Markus ----- May 2 15:09:00 NOTICE[17327]: app_dial.c:1010 dial_exec_full: Unable to create channel of type 'Zap' (cause 0 - Unknown) == Everyone is busy/congested at this time (1:0/0/1) == Auto fallthrough, channel

I am porting a program in matlab to R, The problem is that Matlab has a feature where symbols that aren't arguments are evaluated immediately. That is: Y=3 F=@(x) x*Y Will yield a function such that F(2)=6. If later say. Y=4 then F(2) will still equal 6. R on the other hand has lazy evaluation. F<-function(x){x*Y} Will do the following Y=3 F(2)=6 Y=4 F(2)=8. Does anyone know of away to