library(GART) ################# example ################## set.seed(3434) ss=100 # sample size z=matrix(rgamma(ss,shape=2, rate=2),ss,1) x1=matrix(rbinom(ss,size=1,prob=0.5)) x2=matrix(runif(ss, min=-0.5, max=0.5)) censor1=matrix(NA,ss,1) censor2=matrix(NA,ss,1) cl=1 cu=12 for(i in 1:ss){ w=rbinom(1,1,0.8) censor1[i]=runif(1,min=0,max=cl)*w censor2[i]=runif(1,min=censor1[i],max=cu) } T1=matrix(0,ss,100) T=matrix(13,ss,100) for (i in 1:ss){ T1[i,1]=rexp(1,rate=1) #### Homogeneous Poinon Procen T[i,1]=exp(min(1,T1[i,1]/z[i]/1.5)*x1[i]+x2[i])*T1[i,1]/z[i] #### Nonhomogeneous Poinon Procen j=1 while (T[i,j]n){ Tfinal=matrix(T[ii,(n+1):m]) ID=matrix(ii,m-n,1) break } } for (i in (ii+1):ss) { m=length(which(T[i,]<=censor2[i])) n=length(which(T[i,]n) { Tfinal=rbind(Tfinal,matrix(T[i,(n+1):m])) ID=rbind(ID,matrix(i,m-n,1)) } } id=matrix(1:ss, ss, 1) t0=proc.time() GART(id, x=matrix(cbind(x1, x2), ss, 2), censor1, censor2, ID, T=Tfinal, link="1", grids=150, U=3, inference=TRUE, method="boot", boots=500) t1=proc.time() GART(id, x=matrix(cbind(x1, x2), ss, 2), censor1, censor2, ID, T=Tfinal, link="1", grids=150, U=3, inference=TRUE, method="samp") t2=proc.time() t1-t0 t2-t1