Displaying 2 results from an estimated 2 matches for "nomean".
Did you mean:
noman
2011 Sep 05
1
SAS code in R
...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 incidence rather than Cox model version (crr from cmprsk in R) however there is no pre-set routine for...
2011 Mar 07
0
survest() for cph() in Design package
...-69 or >=70,
used as strata in the model
The SAS code is (if this helps to address my question/problem)
proc phreg data=uminers;
model rstime*cc(0)=cr500 smoke25/ entry=rsentry offset=logw rl;
baseline out=ch covariates=covs cumhaz=cumhaz stdcumhaz=secumhaz lowercumhaz=lci
uppercumhaz=uci/ nomean;
strata rstime2;
run;
My R code is:
library(survival)
library(Design)
> uminers<-read.table("uminers.txt",sep='\t',header=T,row.names=NULL)
> fit <- cph(Surv(rsentry,rstime,cc)~cr500+smoke25+strat(rstime2) +offset(logw),
>uminers, x=T, y=T)
> fit
Cox Propo...