Displaying 1 result from an estimated 1 matches for "lrpvalue".
Did you mean:
llvalue
2009 Oct 05
0
Unusual error while using coxph
...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","Other"),lty=1,color=c(2:5),xlab="Time from randomisation
(y...