【发布时间】:2018-01-11 12:11:43
【问题描述】:
x<-rnorm(100)
u<-runif(100)
t1<-(1/0.0105)*log(1-((0.0105/(exp(-8+0.03*x)))*(log(1-u))))
c<-rexp(100,1/1515)
delta<-1*(t1>c)
datinterval<-data.frame(t1[delta==1],c[delta==1],delta[delta==1])
datinterval$interval <- cut(datinterval$t1, breaks=seq(0, 500, 100))
这是获取 delta=1 数据区间的 r 代码。我想使用以下代码将其插入到包含完整数据的数据框中:
x<-rnorm(100)
u<-runif(100)
t1<-(1/0.0105)*log(1-((0.0105/(exp(-8+0.03*x)))*(log(1-u))))
c<-rexp(100,1/1515)
delta<-1*(t1>c)
delta
dat=data.frame(t1,c,delta)
dat$Interval[delta==1]<-cut(dat$t[delta==1], breaks=seq(0, 500, 100))
代码不显示区间,只显示1,2,3,4...
【问题讨论】: