Displaying 1 result from an estimated 1 matches for "ngra".
Did you mean:
gra
2005 Oct 24
1
Error in step() (or stepAIC) for Cox model
...n 1:nrow(reg.bs)){ #<--bootstrap 10 times
mi<-data.frame(tlfup,cen,complete(imp,j)) #<--completed data sets after
MI
in.star<-sample(1:n,n,replace=T) #<--to sample id number 1-1851.
data.star<-mi[in.star,]
M<-coxph(Surv(tlfup,cen)~mi$trt+mi$nodes+mi$htypeed1+mi$htypeed2+mi$ngra
ded2+mi$agem40
+mi$agem40sq+mi$er+mi$pr,data=data.star)
reg.model<-step(M) #<--do stepwise selection
reg.step[[b]]<-c((reg.model$coef)) #<-store selected variables
chosen.vb[[b]]<-names(reg.step[[b]]) #<--store names of selected
variables
}
tot.vb[[j]]<-c(unlist(chosen.v...