chat_roulette_surv2.R

### Generate our data ###

# Time to seeing a penis
penis<-round(rchisq(30,2))
p_cens<-rep(1,30)
p_type<-rep("penis",30)
p_bind<-cbind(penis,p_cens,p_type)

# Time to seeing a lonely dude
lonely<-round(rchisq(40,1))
l_cens<-rep(1,40)
l_type<-rep("lonely",40)
l_bind<-cbind(lonely,l_cens,l_type)

# Time to seeing a two or more drunk persons
drunk<-sample(1:20,20,replace=T)
d_cens<-rep(1,20)
d_type<-rep("drunk",20)
d_bind<-cbind(drunk,d_cens,d_type)

# Time to seeing one or more women
chick<-sample(70:120,10,replace=T)
c_cens<-rep(1,10)
c_type<-rep("chick",10)
c_bind<-cbind(chick,c_cens,c_type)

# Create data frame
surv_data<-as.data.frame(rbind(p_bind,l_bind,d_bind,c_bind))
colnames(surv_data)<-c("time","censor","type")
surv_data$time<-as.numeric(surv_data$time)
surv_data$censor<-as.numeric(surv_data$censor)