chat_roulette_surv1.R

library(survival)
library(ggplot2)
library(Design)


### Test data generation ###

# Chi-square test
count<-1000
cs1<-rchisq(count,1)
cs2<-rchisq(count,2)

cs_test<-as.data.frame(cbind(cs1,cs2,1:count))
png("cs_test.png",width=1000,height=450,res=100)
ggplot(cs_test,aes(V3))+geom_density(aes(x=cs2,fill="k=2",alpha=0.25))+
    geom_density(aes(x=cs1,fill="k=1",alpha=0.25))+xlim(-1,20)+opts("Chi-square test")
dev.off()