search for: remtime

Displaying 10 results from an estimated 10 matches for "remtime".

Did you mean: reltime
2009 Apr 03
2
Schoenfeld Residuals
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] * sev...
2009 Jun 09
2
Isolating a single plot from plots produced simultaneously
Dear R-Help, I am using the 'mfp' package. It produces three plots (as I am using the Cox model) simultaneously which can be viewed together using the following code: fit <- mfp(Surv(rem.Remtime,rem.Rcens)~fp(age)+strata(rpa),family=cox,data=nearma,select=0.05,verbose=TRUE) par(mfrow=c(2,2)) plot(fit) They can be viewed separately but the return key must be pressed after each graph appears (Click or hit ENTER for next page). I'd like to isolate the second plot produced (the estimated...
2005 Jun 04
6
New HTB-derived qdisc for accounting?
...ed 385953282 bytes 283998 pkt, overlimits 0 requeues 0) Proposal: ---------------------------------- Another idea would be to create a qdisc HTBQ (HTB with quota) derived from HTB with the following characteristics: htb_rate=min(htbq_rate, (alreadysent=>htbq_squota)?((htbq_quota-alreadysent)/remtime):htbq_rate) htb_ceil=htbq_ceil //this is just passed on htb_burst=htbq_burst htb_cburst=htbq_cburst htb_prio=htbq_prio htb_quantum=htbq_quantum //should be set automatically htbq_interval seconds //time after which quota is reset htbq_starttime seconds //time when the first interval starts htbq_qu...
2009 Oct 05
0
Unusual error while using coxph
...type[i] <- 3 if(data$all.type[i]=="G" && data$all.aa[i]>=1 && data$all.tc[i]>=1) stype[i] <- 3 if(data$all.otc[i]>=1) stype[i] <- 4 if(data$all.o[i]>=1) stype[i] <- 4 } return(stype) } fita <- survdiff(Surv(rem.Remtime,rem.Rcens)~stypea(nmarma)+strata(nmrpa),data=nmarma) fita lrpvalue=1-pchisq(fita$chisq,3) xx <- cuminc(nmarma$rem.Remtime/365,nmarma$rem.Rcens,stypea(nmarma),strata=nmrpa) plot(xx,curvlab=c("Simple/Complex","SC+2gentc or 2gentc","TC or My/Ab or My/Ab+gentc","Ot...
2008 Sep 29
1
Hazard plot
Hi All, This sounds a relatively simple query, and I hope it is! I am looking at a continuous variable, age. I am looking at time to 12-month remission and can calculate the HR and 95% confidence interval are follows: coxfita = coxph(Surv(rem.Remtime,rem.Rcens)~nearma$all.age,data=nearma) exp(coxfita$coefficients) exp(confint(coxfita)) However, because I am looking at age as a continuous variable I cannot draw a Kaplan-Meier curve. Instead I need to draw a plot of hazard against age. How do I do this? I don't think plot(nearma$all.age,co...
2008 Sep 30
0
Hazard curves
-- begin included message ----- I am looking at a continuous variable, age. I am looking at time to 12-month remission and can calculate the HR and 95% confidence interval are follows: coxfita = coxph(Surv(rem.Remtime,rem.Rcens)~nearma$all.age,data=nearma) exp(coxfita$coefficients) exp(confint(coxfita)) However, because I am looking at age as a continuous variable I cannot draw a Kaplan-Meier curve. Instead I need to draw a plot of hazard against age. How do I do this? I don't think plot(nearma$all.age,co...
2009 Mar 26
1
Centring variables in Cox Proportional Hazards Model
...ndicator variables, x1 and x2 which represent age categories (x1 is patients less than 16 while x2 is for patients older than 65). If I use the following Cox model, is there anyway I can centre the variables? Do I have to do it before I fit them into the model and if so, how? fit2=coxph(Surv(rem.Remtime,rem.Rcens)~x1(partial)+x2(partial),data=partial,method="breslow") Thank you, Laura
2009 May 22
1
Forcing a variableinto a model using stepAIC
...is stratified by randomisation period). I know that SAS can be used to do this but my SAS coding is poor and consequently it would be easier for me to use R, especially given the fractional polynomial transformations! Currently the model is as follows (without treatment). coxfita=coxph(Surv(rem.Remtime,rem.Rcens)~sind(nearma)+fsh(nearma)+fdr(nearma)+th1(nearma)+th2(nearma)+fp(cage)+fp(fint)+fp(tsb)+strata(rpa),data=nearma) Thank you for your help, Laura
2012 Mar 27
0
Standard error terms from gfcure
Dear R-help, I am using R 2.14.1 on Windows 7 with the 'gfcure' package (cure rate model). I have included the treatment variable in the cure part of the model as shown below: Ø ref_treat <- gfcure(Surv(rem.Remtime,rem.Rcens)~1,~1+strata(drpa)+factor(treat(delcure)),data=delcure,dist="loglogistic") >From that I can obtain the coefficients, standard errors etc as per alternative models (with covariates only fitted to the survival part of the model say). > summary(ref_treat) However, only one...
2008 Oct 15
0
R-help Digest, Vol 67, Issue 31
...f9e8d at mail.gmail.com> Content-Type: text/plain Hi All, This sounds a relatively simple query, and I hope it is! I am looking at a continuous variable, age. I am looking at time to 12-month remission and can calculate the HR and 95% confidence interval are follows: coxfita = coxph(Surv(rem.Remtime,rem.Rcens)~nearma$all.age,data=nearma) exp(coxfita$coefficients) exp(confint(coxfita)) However, because I am looking at age as a continuous variable I cannot draw a Kaplan-Meier curve. Instead I need to draw a plot of hazard against age. How do I do this? I don't think plot(nearma$all.age,co...