similar to: p-value=0 running coxph

Hi, In my analysis of impacts of insecticide-treated bednets on malaria, I look at the relationship between malaria incidence and mosquito behaviors. The condensed data set is copied here. Ordinary regression (lm) shows that Incidence was negatively related to Mortality. This makes sense because the latter reflected the strength of killing mosquitoes by insecticide-treated nets. Since the

hi all, I have a simple problem that i am not able to solve. I've a list called datalist with the following structure: [...] [[10]] [[10]]$a -1 0 1 -1 31 5 2 0 6 7 5 1 1 7 36 [[10]]$b -1 0 1 -1 31 5 2 0 6 7 5 1 1 7 36 [[10]]$c [1] 0.855 [[10]]$d [1] 0.855 [...] with [[1]] ... [[100]]. How can i get the mean value of datalist[[x]]$d, where x

Hi, When I run the survfit function, I want to get the restricted mean value and the standard error also. I found out using the "print" function to do so, as shown below, print(km.fit,print.rmean=TRUE) Call: survfit(formula = Surv(diff, status) ~ 1, type = "kaplan-meier") records n.max n.start events *rmean *se(rmean) median 200.000

great!! thanks very much, mean(unlist(lapply(listdata, function(z) z$c))) works well. and what about getting the average table $a (displaying the average elements across all 1000 matrix)? could you please help me? I am struggling with this... thanks in advance Roberto mean(unlist(lapply(x, function(z) z$d))) should do the trick On Sun, 12 Feb 2006 20:06:12 +0000, statistical.model at

Dear R-tisans, I am trying to calculate the 12th root of a transition (square) matrix, but can't seem to obtain an accurate result. I realize that this post is laced with intimations of quantitative finance, but the question is both R-related and broadly mathematical. That said, I'm happy to post this to R-SIG-Finance if I've erred in posting this to the general list. I've

Hello, I am hoping for some advice regarding warning/error messages I received when running a Cox regression # message 1 - obtained while creating a plot of residuals > plot (NV.zph, main = "groupNUSM - UNFIT", var= 'groupNUSM') Warning messages: 1: In approx(xx, xtime, seq(min(xx), max(xx), length.out = 17)[2 * : collapsing to unique 'x' values 2: In

Hello dear all, I am unable to understand why when I run the following three lines: set.seed(4254) > a <- data.frame(y = rnorm(40), x=ordered(sample(1:5, 40, T))) > summary(lm(y ~ x, a)) The output I get includes factor levels which are not relevant to what I am actually using: Call: > lm(formula = y ~ x, data = a) > Residuals: > Min 1Q Median 3Q Max >

Hi, I''ve used old ZPH patch under squid 2.4 Stable4 and it works great ! Now I want to patch squid 2.4 stable 5, with new patch, on http://www.it-academy.bg/zph/ I''ve patched and installed squid 2.5 stable 5 succefully, but I can''t get ZPH works. I''m trying with ... $TC class add dev $LANDEV parent 1: classid 1:7 htb rate 1Mbit $TC filter add dev $LANDEV

I've been Googling, but have not come up with anything. It worked before the update, bit now samba won't start. I see the following in the logs. -- Unit smb.service has begun starting up. Apr 22 18:51:47 001cognos01.ccnva.local systemd[1]: smb.service start operation timed out. Terminating. Apr 22 18:53:18 001cognos01.ccnva.local systemd[1]: smb.service stop-final-sigterm timed out.

You report an error message: > plot(zph.revasFit[1]) Error in plot.window(xlim, ylim, log, asp, ...) : need finite 'ylim' values I have never seen this error before, and I cannot guess what causes it. You need to provide more information, and likely a small data set that produces the problem. Perhaps you have an x variable that is a constant? Terry Therneau

Dear all, Using the data set and code below, I am interested in modelling how egg temperature (egg.temp) is related to energy expenditure (kjday) and clutch size (treat) in incubating birds using the lme-function. I wish to generate the F-distribution for my model, and have tried to do so using the anova()-function. However, in the resulting anova-table, the parameter kjday has gone from being

The following works, you need to include x=TRUE in the call to coxph. Passing the time and status variables as additional arguments is a matter of personal preference. f.coxph.zph<-function(x, timeVar, statusVar) { cox.fit <- coxph(Surv(timeVar, statusVar) ~ x, na.action = na.exclude, method = "breslow", x=TRUE) fit.zph<-cox.zph(cox.fit) fit.zph$table[,3] } time.cox <-

Hej hej, is there a way to extend the SpatialPointsDataFrame data slot?This is the structure of an object of it:> str(coord) Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots ..@ data :'data.frame': 214 obs. of 2 variables: .. ..$ location.long: num [1:214] -79.8 -79.8 -79.8 -79.8 -79.8 ... .. ..$ location.lat : num [1:214] 9.16 9.16 9.16

Hi, I wrote a function that worked well in R, but not in S-Plus, can anyone suggest a solution? > f.coxph.zph<-function(x) { cox.fit <- coxph(Surv(time.cox, status.cox) ~ x, na.action = na.exclude, method = "breslow") fit.zph<-cox.zph(cox.fit,transform='log') fit.zph$table[,3] } yyy is my data frame that contains survial time, censor status and predictor

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

Hello, Excuse me for a more methodological than technical question. I'm developing a Cox model with 10 covariates. One of them is age (named "eta"). I've checked proportionality with cox.zph with age continuous: > cox.zph(coxph(Surv(TTP,CENSOTTP)~eta)) rho chisq p eta -0.0154 0.0225 0.88 and categorical (eta<60): >

Hi, I used to patch my squid with ZPH patch on http://www.it-academy.bg/zph/ > The idea behind this patch is to allow classification > of packets generated from the squid cache engine towards > clients. > The classification is based on whether the content is > being served from cache (a cache HIT), or > is being retrieved from a remote server (a cache MISS). Very useful

Hi, I''m trying to share internet on a LAN I''ve a linux router with SQUID (with ZPH support) + FAIRNAT The idea is: - fairness sharing internet - priorize interactive traffic - if a web object is on squid-cache (HIT), user can download it, with a rate = LAN rate I''ve: - Last Fairnat Script: www.metamorpher.de/fairnat/ - SQUID 2.5STABLE5 with ZPH patch

Hi, I am using the Survival package, more precisely the cox.zph function, to plot log(Hazard rate) over time. if I type plot(temp.zph[2]) then I get the plot I want. However, I want to change the label of the y axis that cox.zph prints. plot(temp.zph[2],xlab='Days', ylab='log hazard for Fast recovery cluster') then I get Error in plot.default(range(xx), yr, type =

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