search for: agem40

...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.vb)) } Error...

Thank you for Prof.Ripley's suggestion. I fixed the program by adding a lower scope, and the program ran, but I still got warning messages, and don't know what is going on, would this affect my results? ... Step: AIC= 12337.74 Surv(tlfup, cen) ~ MI[[j]]$trt + MI[[j]]$agem40 + MI[[j]]$agem40sq + mhtypeed1 + mhtypeed2 Df AIC <none> 12338 - MI[[j]]$agem40sq 1 12338 - MI[[j]]$agem40 1 12339 - mhtypeed2 1 12353 - mhtypeed1 1 12365 There were 50 or more warnings (use warnings() to see the first 50) Warnin...

...rs are registered that can debug this exception. Unable to JIT debug. In my program, I first source("sump6d2.q") to get data from splus, and attach(sump6d2), Whenever I run this part of program, the error pops up: library(mix) Y<-cbind(race2,nodes,htypeed1,htypeed2,ngraded2,blood2,agem40,agem40sq,e r,pr) #<--Y has to be a matrix for mix s<-prelim.mix(Y,6) #----------------------------run multiple imputation------------------------------# MI<-vector("list",10) #<--vector of complete data after MI fit.model.mi<-vector("list",10) rngseed(1234567...