Displaying 10 results from an estimated 10 matches for "surv2".
Did you mean:
surv
2009 Oct 13
1
Time Dependent Cox Model
...k it is working properly except for calculating
the death.time. This column should be 1s or 0s and as I have it only
produces 0s. Any help will be greatly appreciated.
http://www.nabble.com/file/p25881478/Survival_master2.csv
Survival_master2.csv
Here is my code:
sum(!is.na(surv[,16:726]))
surv2<-matrix(0,12329,19)
colnames(surv2)<-c('start', 'stop', 'death.time',
names(surv)[1:15],'habitat')
row<-0 # set record counter to 0
for (i in 1:nrow(surv)) { # loop over individuals
for (j in 16:726) { # loop over 726 days
if (is.na(sur...
2015 Aug 02
3
ayuda con análisis de supervivencia
...link con los datos (.Rdata) y el script que he utilizado en el análisis.
https://www.dropbox.com/s/d96itird8ms42yx/dataframe.Rdata?dl=0 <https://www.dropbox.com/s/d96itird8ms42yx/dataframe.Rdata?dl=0>
sapply(levels(df0$bmi0),function (x){ #####SURVIVAL CURVE
dfx=filter(df0,bmi0==x)
surv2=Surv(dfx$bmi,dfx$MetS)
km2=survfit(surv2~dfx$apoe4)##start.time=20,type='kaplan')
plot(km2,lty=2:1,xlim=c(20,41),xlab='BMI at onset',main=x,mark.time = F)
legend('bottomleft',c('E4','no-E4'),lty=2:1)
cox=list(coxph(surv2~relevel(dfx$apoe4,ref='no-...
2009 Feb 18
1
Age as time-scale in a cox model-How to calculate x-time risk?
...ld be:
stanf.surv1=survest(stanf.cph1, times=1000)
#Age in this case is a covariate.
#I now want to compare the above estimate to the 1000-day survival probability I get using age at entry and exit as my time-scale:
stanf.cph2=cph(Surv(age,age+time, status) ~ t5+id, data=stanford2, surv=T)
stanf.surv2=survest(stanf.cph2, times=1000)
summary(stanf.surv1$surv)
Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
0.1131 0.3370 0.4669 0.4538 0.5633 0.7480 27.0000
> summary(stanf.surv2$surv)
Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
0.07387 0.23240 0.35...
2011 Sep 05
1
SAS code in R
...als for combinations of risk factors for my outcome.
/* Combinations of Risk Factors */
data test2;
input sex treat;
DATALINES;
0 0
1 0
0 1
1 1
;
run;
/* Survival estimates for the above combinations */
proc phreg data = pudat2;
model withtime*wcens(0) = sex treat /ties = efron;
baseline out = surv2 survival = survival lower = slower upper = supper
covariates = test2 /method = ch nomean cltype=loglog;
run;
/* Survival estimates at 1 year */
proc print data = surv2 noobs;
where withtime = 364;
run;
I now would like to do the same thing but in a competing risks setting (cumulative inciden...
2008 Feb 28
0
surv2sample 0.1-2
Dear useRs,
There is a new version 0.1-2 of the package surv2sample available on CRAN.
Users of the previous versions should update because a bug in the function
cif2.ks has been fixed.
General information about the package:
surv2sample provides various two-sample tests for right-censored survival
data. Three main areas and corresponding methods are:...
2008 Feb 28
0
surv2sample 0.1-2
Dear useRs,
There is a new version 0.1-2 of the package surv2sample available on CRAN.
Users of the previous versions should update because a bug in the function
cif2.ks has been fixed.
General information about the package:
surv2sample provides various two-sample tests for right-censored survival
data. Three main areas and corresponding methods are:...
2010 Jul 28
1
anderson-darling test
...of fit test. But I'm not sure which
package in R to use.
I tried ad.test(...) but it does not recognise the test by Vito Ricci in FITTING
DISTRIBUTIONS WITH R
> ad.test(hist_hume_beec[,1],hist_hume_beec[,2])
Error: could not find function "ad.test"
I also tried stat.ad from surv2.ks {surv2sample} but not really sure how to use
it.
Here are the data:
rel_obs rel_gen
[1,] 0.000000000 0.0000000000
[2,] 0.052576236 0.0282122554
[3,] 0.085173502 0.0752419321
[4,] 0.080967403 0.0986985317
[5,] 0.091482650 0.1060457970
[6,] 0.092534175 0.1035785945
[7,] 0.0...
2009 Feb 19
0
Age as time-scale in a cox model
...3 year survival of 45%. I was surprised when you put "id" in the model, but
it turns out to have p=.03! It seems that patients entered later in the study
have better survival.
Now for age scale:
> fit2 <- coxph(Surv(age, age+ time/365.25, status) ~ t5 + id, stanford2)
> surv2<- survfit(fit2)
> surv2
n events median 0.95LCL 0.95UCL
1.0 102.0 12.2 12.2 28.1
This shows a median age at death of 12.2 years. Puzzling, isn't it.
First, note that your code
cph(Surv(age,age+time, status) ~ t5+id, data=stanford2...
doesn't make sen...
2024 Feb 07
2
Difficult debug
...rs advise as a next step?
>
> Here is the last part of the screen
> > fit2 <- coxph(list(Surv(tstart, tstop, bstat) ~ 1,
> + c(1:4):5 ~ age / common + shared), id= id,
> istate=bili4,
> + data=pbc2, ties='breslow', x=TRUE)
> > surv2 <- survfit(fit2, newdata=list(age=50), p0=c(.4, .3, .2, .1, 0))
> > test2 <- mysurv(fit2, pbc2$bili4, p0= 4:0/10, fit2, x0 =50)
> ==31730== Invalid read of size 8
> ==31730== at 0x298A07: Rf_allocVector3 (memory.c:2861)
> ==31730== by 0x299B2C: Rf_allocVector (Rinlinedfu...
2024 Feb 07
2
Difficult debug
...a for( ) loop.??? What
would others advise as a next step?
Here is the last part of the screen
> fit2 <- coxph(list(Surv(tstart, tstop, bstat) ~ 1,
+??????????????????? c(1:4):5 ~ age / common + shared), id= id, istate=bili4,
+?????????????? data=pbc2, ties='breslow', x=TRUE)
> surv2 <- survfit(fit2, newdata=list(age=50), p0=c(.4, .3, .2, .1, 0))
> test2 <- mysurv(fit2, pbc2$bili4, p0= 4:0/10, fit2, x0 =50)
==31730== Invalid read of size 8
==31730==??? at 0x298A07: Rf_allocVector3 (memory.c:2861)
==31730==??? by 0x299B2C: Rf_allocVector (Rinlinedfuns.h:595)
==31730==?...