This is probably a stupid question but I need help.
I want to calculate the null distribution for the gene overlap between 2 lists.
The first list is Chip-seq data and the second RNA-seq. And the background genome is 20000 thousand genes. I have this data for 50 samples.
The first sample has a Chip list with 751 genes and a 590 RNA-seq gene list.
I tried it with r but the result looks odd.
ts = replicate(5000,t.test(rnorm(751),rnorm(590))$statistic) range(ts) pts = seq(-3.5, 3.5,length=100) plot(pts,dt(pts,df=25),col='red',type='l') lines(density(ts))