Distribution

From ZhangLabWiki
Revision as of 18:16, 7 July 2015 by >Shicheng (Created page with " table.binomial.test<-function(region,event.time){ region1<-table(region) region2<-unique(region) p=length(region)/(length(region2)*event.time) pvalue<-c() for(i in ...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
table.binomial.test<-function(region,event.time){
 
region1<-table(region)
region2<-unique(region)
p=length(region)/(length(region2)*event.time)
pvalue<-c()
for(i in 1:length(table(region))){
 ptmp<-binom.test(table(region)[i],event.time,p,alternative="two.sided")$p.value
 pvalue<-c(pvalue,ptmp)
}
return(pvalue)
print(paste("there are",sum(pvalue<0.05/length(pvalue),"regions were over-preferred in the samples"),sep=" "))
}